Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,133,952 C→T 6.1% Q40* (CAA→TAA)  katG → catalase‑peroxidase HPI, heme b‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,133,9520CT 6.1% 208.9 / 2.5 65Q40* (CAA→TAA) katGcatalase‑peroxidase HPI, heme b‑containing
Reads supporting (aligned to +/- strand):  ref base C (27/34);  new base T (2/2);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

AACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGCCTGAAAAAAGATCTGAAAGCCCTGTTGACAGAATCTCAACCGTG  >  NC_000913/4133805‑4134097
                                                                                                                                                   |                                                                                                                                                 
aaCGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATc                                                                                                                                                   >  6:164338/1‑148 (MQ=255)
     tAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTc                                                                                                                                             >  4:293236/1‑149 (MQ=255)
       acacTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCg                                                                                                                                            <  6:118206/148‑1 (MQ=255)
        cacTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCgtg                                                                                                                                          <  2:159762/149‑1 (MQ=255)
             tAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGAc                                                                                                                                     <  1:147667/149‑1 (MQ=255)
             tAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGAc                                                                                                                                     >  1:341365/1‑149 (MQ=255)
              agagGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATTAACTTCGTGTTGAcc                                                                                                                                    >  4:180033/1‑149 (MQ=255)
              agagGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATTAACTTCGTGTTGAcc                                                                                                                                    >  6:200021/1‑149 (MQ=255)
                         caTTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCaa                                                                                                                         >  2:15088/1‑149 (MQ=255)
                          aTTGATGAGCACGTCAGACGATATCCATAACCCCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCaa                                                                                                                         <  6:165554/148‑1 (MQ=255)
                           ttGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATTAACTTCGTGTTGACCTGTTAAACCAACa                                                                                                                       <  3:180033/149‑1 (MQ=255)
                           ttGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATTAACTTCGTGTTGACCTGTTAAACCAACa                                                                                                                       <  5:200021/149‑1 (MQ=255)
                               tgaGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCt                                                                                                                   <  3:293236/149‑1 (MQ=255)
                                   cACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATc                                                                                                               <  1:160067/149‑1 (MQ=255)
                                   cACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTCGACCTGTTAACCCAACATTCTAATc                                                                                                               >  8:157141/1‑149 (MQ=255)
                                         aGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCt                                                                                                          <  7:291460/148‑1 (MQ=255)
                                         aGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTa                                                                                                         <  4:130339/149‑1 (MQ=255)
                                         aGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGG‑‑CAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTaa                                                                                                        <  3:42663/148‑1 (MQ=255)
                                          gACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAGCCAACATTCTAATCGTTCTaa                                                                                                        >  2:347578/1‑149 (MQ=255)
                                          gACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTaa                                                                                                        >  3:328819/1‑149 (MQ=255)
                                           aCGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAAc                                                                                                       <  2:282883/149‑1 (MQ=255)
                                              atatCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCa                                                                                                    <  8:15649/149‑1 (MQ=255)
                                                  ccATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTgg                                                                                                >  8:15519/1‑149 (MQ=255)
                                                    aTAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGTCTGGTGGCCAAATCAACTTCGTGTTGACCTGTTTAACCAACATTCTAATCGTTCTAACCCACTGGGt                                                                                              <  1:288446/149‑1 (MQ=255)
                                                     tAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGTGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGt                                                                                              >  4:91617/1‑148 (MQ=255)
                                                     tAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGt                                                                                              >  4:170212/1‑148 (MQ=255)
                                                       acacCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAg                                                                                           <  5:23518/149‑1 (MQ=255)
                                                         accacAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGa                                                                                         <  3:317861/149‑1 (MQ=255)
                                                                ccACTGGCAAATGCCCGTTCCATCCGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGAc                                                                                  <  6:279693/149‑1 (MQ=255)
                                                                  aCTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTa                                                                                <  5:305998/149‑1 (MQ=255)
                                                                   cTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTAc                                                                               >  2:7850/1‑149 (MQ=255)
                                                                    tGGCAAATGCCCGTTCCATCAGGGCGGTCACGCCCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTATAACCAACTGGTGGAGGACTTTGACTAcc                                                                              >  3:68084/1‑149 (MQ=255)
                                                                    tGGCAAATGCCCGTTCCATAAGGGCGGTCCCGACCAGAGTGCGGGGGCGGGCACAACCCCTCGCGACGGGTGGCCAAATCAACTTCGTGGTGACCTGTTAACCCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGCCTAcc                                                                              >  3:63584/1‑149 (MQ=255)
                                                                         aaTGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCaaa                                                                         <  1:15088/149‑1 (MQ=255)
                                                                             cccGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGaa                                                                      <  6:214015/148‑1 (MQ=255)
                                                                              ccGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGCCCTGTTAAACCAACAGTCTAATCGTTCTAACCCACTGGGGGAGTACTTTGACTACCGCATAGAAtt                                                                    >  8:90976/1‑149 (MQ=255)
                                                                              ccGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAAtt                                                                    >  3:155421/1‑149 (MQ=255)
                                                                                   ccATCAGGGCGGTCACGACCAGAGGGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGc                                                                <  1:299493/148‑1 (MQ=255)
                                                                                    cATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCa                                                               <  6:157988/148‑1 (MQ=255)
                                                                                    cATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACAATCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCaa                                                              <  5:164338/149‑1 (MQ=255)
                                                                                     aTCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTACCCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCaaa                                                             >  4:93887/1‑149 (MQ=255)
                                                                                                 aCGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATtactac                                                 >  7:342456/1‑149 (MQ=255)
                                                                                                    aCCAGAGTGCGGGGGCGGGCACAACCACTCTCGACTTGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGTGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGc                                              <  6:107447/149‑1 (MQ=255)
                                                                                                             cGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGCGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGCCTGaaaaaa                                     <  4:328819/149‑1 (MQ=255)
                                                                                                             cGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGCCTGaaaa                                       >  5:97510/1‑147 (MQ=255)
                                                                                                              gggggCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAAAATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGCCTGAAAAAAg                                    <  2:341365/149‑1 (MQ=255)
                                                                                                                     ggcACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGCCTGAAAAAAGATCTGaa                             <  7:15519/149‑1 (MQ=255)
                                                                                                                      gcACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGCCTGAAAAAAGa                                   <  5:267763/142‑1 (MQ=255)
                                                                                                                      gcACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGCCTGAAAAAAGa                                   >  6:267763/1‑142 (MQ=255)
                                                                                                                      gcACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGCCTGAAAAAAGATCTGaaa                            <  3:170212/149‑1 (MQ=255)
                                                                                                                            ccACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGAt                                                       >  8:77662/1‑116 (MQ=255)
                                                                                                                            ccACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGAt                                                       <  7:77662/116‑1 (MQ=255)
                                                                                                                            ccACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGCCTGAAAAAAGATCTGAAAGCCCTg                      <  8:342456/149‑1 (MQ=255)
                                                                                                                            ccACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGCCTGAAAAAAGATCTGAAAGCCCTg                      <  1:7850/149‑1 (MQ=255)
                                                                                                                             cACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGGGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGCCTGAAAAAAGATCTGAAAGCCCTGt                     >  7:256238/1‑149 (MQ=255)
                                                                                                                                tCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCAATGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCATATTAGATTACTACGGCATGAAAAAAGATCTGAACGCCATGTTGa                  >  1:279113/1‑149 (MQ=255)
                                                                                                                                  gcgACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGAt                                                       <  6:183617/110‑1 (MQ=255)
                                                                                                                                  gcgACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGAt                                                       >  5:183617/1‑110 (MQ=255)
                                                                                                                                  gcgACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGCCTGAAAAAAGATCTGAAAGCCCTGTTGACa                <  8:251531/149‑1 (MQ=255)
                                                                                                                                       tggtggCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGCCTGAAAAAAGATc                                 <  3:44945/127‑1 (MQ=255)
                                                                                                                                       tggtggCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGCCTGAAAAAAGATc                                 >  4:44945/1‑127 (MQ=255)
                                                                                                                                            gCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGCCTGAAAAAAGATCTGAAAGCCCTGTTGACAGAATCTCAAc      >  4:51503/1‑149 (MQ=255)
                                                                                                                                              cAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTATTCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACg                                                >  3:291083/1‑105 (MQ=255)
                                                                                                                                              cAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTATTCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACg                                                <  4:291083/105‑1 (MQ=255)
                                                                                                                                               aaaTCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGCCTGAAAAAAGATCTGAAAGCCCTGTTGACAGAATCTCAACCgt   >  4:292833/1‑149 (MQ=255)
                                                                                                                                                aaTCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACACACTGGGTGAGGACTTTGACTACCGCAAAGAATTAAGAAAATTAGATTACTACTGCCTGAAAAAAGATCTGAAAGCACTGTTGACAGAATCTCAACCgtg  >  3:241508/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAGACGATATCCATAACACCACAGCCACTGGCAAATGCCCGTTCCATCAGGGCGGTCACGACCAGAGTGCGGGGGCGGGCACAACCACTCGCGACTGGTGGCCAAATCAACTTCGTGTTGACCTGTTAAACCAACATTCTAATCGTTCTAACCCACTGGGTGAGGACTTTGACTACCGCAAAGAATTCAGCAAATTAGATTACTACGGCCTGAAAAAAGATCTGAAAGCCCTGTTGACAGAATCTCAACCGTG  >  NC_000913/4133805‑4134097

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑