Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,852,038 A→C 7.5% Q51P (CAG→CCG) ‡ pncA → nicotinamidase/pyrazinamidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,852,0380AC 7.5% 220.2 / 8.2 80Q51P (CAG→CCG) ‡pncAnicotinamidase/pyrazinamidase
Reads supporting (aligned to +/- strand):  ref base A (36/38);  new base C (2/4);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.34e-01

GCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTG  >  NC_000913/1851889‑1852186
                                                                                                                                                     |                                                                                                                                                    
gCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAACCGGTTATCGCCAGTCa                                                                                                                                                      >  2:593492/1‑150 (MQ=255)
       cgcgCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTgg                                                                                                                                               <  2:1728334/150‑1 (MQ=255)
        gcgcCCTGTTACTGGTCGCTTTACAAAATGAGATCTGGGCTGGTGGCGCGCTCGCCGGGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCCTGACTGGc                                                                                                                                              <  2:1081174/150‑1 (MQ=255)
         cgcCATGTCACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCCTGACTGGcg                                                                                                                                             <  1:852653/150‑2 (MQ=255)
           cccTGTTACTGGTCGATTTACAAAATGATTTCTGTGCGGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCCTGACTGGcgat                                                                                                                                           <  1:383850/150‑4 (MQ=255)
              tGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCCTGACTGGcgataac                                                                                                                                        <  1:1216546/150‑7 (MQ=255)
              tGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCgg                                                                                                                                        >  1:1695581/1‑150 (MQ=255)
                ttACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCg                                                                                                                                      <  1:953813/150‑1 (MQ=255)
                     ggTCGATTAACAAAAAGAATTCTGTCCTGGTAAAGCGCTCGCAGTGCCGGAAGGTGACAGTACGGTGGATGTAGCTAACCGCCTGAATGAGTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGTGTGGCCCCCGGCGaagcc                                                                                                                                 >  2:280973/1‑147 (MQ=255)
                          agttACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACTGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCa                                                                                                                            <  2:841554/148‑1 (MQ=255)
                           tttACAACATGATTGCTGTGCGGGTGGCGCGCTCGCCGTGCCGGAACGTGACAGTACGGGGGCTGGCGCTCCCCGCCGGCTTGACTGGTGCCAGGCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAg                                                                                                                           <  1:112509/150‑1 (MQ=255)
                             tACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGACGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGtt                                                                                                                         >  2:173309/1‑150 (MQ=255)
                                  aaTGATTTCTTTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGcc                                                                                                                    <  2:1569786/150‑1 (MQ=255)
                                  aaTGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGcc                                                                                                                    <  2:1576757/150‑1 (MQ=255)
                                    tGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCCGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCCCCCGGCGAATCCCGGCAGTTTTGCGAg                                                                                                                  >  2:118139/1‑150 (MQ=255)
                                       tttCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCa                                                                                                               >  1:23646/1‑150 (MQ=255)
                                         tCTGTGCGGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGGCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGGTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGc                                                                                                             <  2:348009/150‑1 (MQ=255)
                                          cTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTACCCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTCTCGCCAGTCAGGACTGGCACCCGGCGAATCCCGGCAGTTTTGCCAGTCAGcc                                                                                                            >  1:102514/1‑149 (MQ=255)
                                          atgtgCTGGTGGCGCGATCGCCGTGCCGGAAGGAGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCGAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCa                                                                                                            >  1:102510/2‑150 (MQ=255)
                                             tgCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACgg                                                                                                         <  1:118914/150‑1 (MQ=255)
                                               cTGGTGGCGCGCTCGCCGGGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGgtg                                                                                                       <  2:1420639/150‑1 (MQ=255)
                                                    ggCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTACCCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGGCAGCACGGTGTagag                                                                                                  >  1:95717/1‑150 (MQ=255)
                                                    ggCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTagag                                                                                                  <  1:1202204/150‑1 (MQ=255)
                                                           tCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTtata                                                                                           <  1:1025838/150‑1 (MQ=255)
                                                            cGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATAc                                                                                          >  1:691460/1‑150 (MQ=255)
                                                                gTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCa                                                                                      >  1:1373520/1‑150 (MQ=255)
                                                                 ggCCGGAAGGGGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGGTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAg                                                                                     <  1:327505/149‑1 (MQ=255)
                                                                    cGGAAGGTGACAGTACGGTGGATGGCGCTAACCGCCTGATTGACTGGTGCCAGGCGCGCGGTGAAGCGGGTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGcc                                                                                  <  1:603435/150‑1 (MQ=255)
                                                                     ggAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCa                                                                                 >  1:409842/1‑150 (MQ=255)
                                                                        aGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACt                                                                              <  1:871732/150‑1 (MQ=255)
                                                                        aGGTGACAGTACGGGGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACt                                                                              <  2:1443387/150‑1 (MQ=255)
                                                                            gACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGAc                                                                          <  2:1670649/150‑1 (MQ=255)
                                                                            gACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGGGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGAc                                                                          <  1:1241571/150‑1 (MQ=255)
                                                                             aCAGTACGGTGGACGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACg                                                                         <  2:36302/150‑1 (MQ=255)
                                                                              cAGTACGGTGTATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACt                                                                              <  2:160485/144‑1 (MQ=255)
                                                                              cAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACt                                                                              >  1:160485/1‑144 (MQ=255)
                                                                              cAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACgg                                                                        >  1:514995/1‑150 (MQ=255)
                                                                               aGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGt                                                                       >  1:1624023/1‑150 (MQ=255)
                                                                                gTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGTGAATCACGGCAGTTTTGCCTGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGtt                                                                      >  2:970390/1‑150 (MQ=255)
                                                                                  aCGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTg                                                                    >  1:28344/1‑150 (MQ=255)
                                                                                  aCGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCAAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTg                                                                    >  1:1591025/1‑150 (MQ=255)
                                                                                   cGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGc                                                                   <  1:781069/150‑1 (MQ=255)
                                                                                      gggATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCcac                                                                <  2:1353019/149‑1 (MQ=255)
                                                                                       ggATGTCGCTAACCGCCTGATTGACTGGTGCCAGACGAGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGTCAGTTTTGCCAGTCAGCTCGGTGTAGAGCCTTATACGCCAGGCCAAATAGACGGTTTGCcaca                                                               >  1:180094/1‑150 (MQ=255)
                                                                                           gTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAAcc                                                           >  2:1409710/1‑150 (MQ=255)
                                                                                                 aaCCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTgg                                                     >  2:726317/1‑150 (MQ=255)
                                                                                                 aaCCGCCTGAGTGAGTGGTGCCAGGCGGGCGGTGAAGCGGTTATCGTCAGTCAGGACTGGCACCCGGCGAATCACGGGAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACTATCCTTCTgg                                                     <  2:682831/150‑1 (MQ=255)
                                                                                                  cccGCCAGATTGACTGGTGCCAGCCGCGCGGTGAACCGGGTATCTCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGTGGCAAACCGTcgagt                                                    <  2:1153747/149‑5 (MQ=255)
                                                                                                  aCCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGc                                                    <  1:143161/150‑1 (MQ=255)
                                                                                                   ccGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGcc                                                   <  2:1343593/150‑1 (MQ=255)
                                                                                                     gcatGATTGTCTGGTGCCAGTCGCGAGGTGATGCTGTTATCGCCAGTGAGGACTGTCACCCGGCGAATAACGGCAGTTTTGTCAGTCAGCACGGTGTAGAGCCTTATACTCCAGGCCAACTCGACGGTTTACCACAAACCTTCTGGCCAg                                                 <  1:292438/147‑1 (MQ=255)
                                                                                                     gCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAg                                                 >  1:1334921/1‑150 (MQ=255)
                                                                                                      ccTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGa                                                >  2:18258/1‑150 (MQ=255)
                                                                                                      ccTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATAGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGa                                                <  1:1584151/150‑1 (MQ=255)
                                                                                                         gATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCa                                             >  2:741166/1‑150 (MQ=255)
                                                                                                          aTTGACTGCTGCCCGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCAc                                            <  2:291920/150‑1 (MQ=255)
                                                                                                             gACTGGTGCCTGTCGCGCGGTGATGCGGTTATCGGCAGTCAGGACTGGCTCCCGGAGAATCACGGCAGTCTTTCCAGTCAGCACGGTGTAGAGCCTTAGACGCCAGGCCATCTCGAGGGGTTGCCACAAACCTTCTGGCCAGATCACtgt                                         <  1:280973/150‑1 (MQ=255)
                                                                                                             gACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGCATCACGGCAGTTTTGCCAGTCAGCACGGTGTTGCGCCTTATCCGCCAGGCCAACTCGACGGTTTACCACAAACCTTCTGGCCAGATCACtgt                                         >  1:274379/1‑150 (MQ=255)
                                                                                                                    gCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGCGCCTTATACGCCAGGCCAACTCGACGGTTTGCCGCAAGCCTTCTGGCCAGATCACTGTGTGCAGa                                  >  1:1058281/1‑150 (MQ=255)
                                                                                                                    gCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGa                                  >  1:728682/1‑150 (MQ=255)
                                                                                                                      cAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAAc                                <  2:766764/150‑1 (MQ=255)
                                                                                                                      cAGTCGCGCGGGGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAAc                                <  1:914956/150‑1 (MQ=255)
                                                                                                                        gTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAg                              <  1:726317/150‑1 (MQ=255)
                                                                                                                         tCGCGCGGTGAAGCGGGTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGt                             <  2:615587/150‑1 (MQ=255)
                                                                                                                          cgcgcgGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACACACCTTCTGGCCAGATCACTGTGTGCAGAACAGTg                            >  1:693815/1‑150 (MQ=255)
                                                                                                                          cgcgcgGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACATTg                            >  2:1317629/1‑150 (MQ=255)
                                                                                                                          cgcgcgGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTg                            >  2:403070/1‑150 (MQ=255)
                                                                                                                            cgcgGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGaa                          <  1:542725/150‑1 (MQ=255)
                                                                                                                              cgGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAgg                        <  1:689842/150‑1 (MQ=255)
                                                                                                                                gTGAAACGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCCGCACGGTGTAGAGCCTTATCCGCCAGGCCAACTCGACGGTTTCCCACACACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGgcg                      >  1:262405/1‑150 (MQ=255)
                                                                                                                                   aaGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTGGCCAGTCAGCACGGGGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTGCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGcac                   <  1:173309/150‑1 (MQ=255)
                                                                                                                                    aGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCTACTGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGcaca                  >  1:1307616/1‑150 (MQ=255)
                                                                                                                                       ggTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCCCAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAAtt               >  1:122333/1‑150 (MQ=255)
                                                                                                                                       ggTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAAtt               <  2:1639905/150‑1 (MQ=255)
                                                                                                                                       ggTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAAtt               <  1:700792/150‑1 (MQ=255)
                                                                                                                                         cgggtgCCAGTCCTGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGCGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTAc             >  2:849308/6‑150 (MQ=255)
                                                                                                                                          tATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCGTCTGGCCAGATCACTGTGTGCAGCACAGTGAAGGCGCAGAATTACa            >  1:1272620/1‑150 (MQ=255)
                                                                                                                                           ggtgCCAGTCCTGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAATCCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACAt           >  1:1011057/4‑150 (MQ=255)
                                                                                                                                           aTCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACAt           >  2:841307/1‑150 (MQ=255)
                                                                                                                                            tCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGCTCACTGTGTGCTGAACAGTGAAGGCGCACAAttccgtc          >  1:946244/1‑145 (MQ=255)
                                                                                                                                                    cAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTg  >  2:931364/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTG  >  NC_000913/1851889‑1852186

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