Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,852,039 G→T 7.5% Q51H (CAG→CAT) ‡ pncA → nicotinamidase/pyrazinamidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,852,0390GT 7.5% 221.6 / 6.6 80Q51H (CAG→CAT) ‡pncAnicotinamidase/pyrazinamidase
Reads supporting (aligned to +/- strand):  ref base G (36/38);  new base T (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 = 5.69e-01

CGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTG  >  NC_000913/1851896‑1852186
                                                                                                                                               |                                                                                                                                                   
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)
                                                                                                                                               |                                                                                                                                                   
CGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTG  >  NC_000913/1851896‑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: