Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,911,917 G→A 5.3% Q579* (CAA→TAA) ‡ relA ← GDP/GTP pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,911,9170GA 5.3% 295.2 / 2.6 92Q579* (CAA→TAA) ‡relAGDP/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base G (41/46);  new base A (3/2);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCTGGGTAATGAAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGAT  >  NC_000913/2911768‑2912066
                                                                                                                                                     |                                                                                                                                                     
ccccTGGGTAATGAAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTg                                                                                                                                                       >  2:920403/1‑150 (MQ=255)
ccccTGGGTAATGAAGCCGACAATCTCATCTCCAGGAATAGGCTGGAAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTg                                                                                                                                                       >  1:545423/1‑150 (MQ=255)
ccccTGGGTAAGGAAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTg                                                                                                                                                       <  2:1034629/150‑1 (MQ=255)
 cccTGGGTAATGAAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGc                                                                                                                                                      >  1:160746/1‑150 (MQ=255)
  ccTGGGTAATGAAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCt                                                                                                                                                     <  2:627003/150‑1 (MQ=255)
    tGGGTAATGAAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTc                                                                                                                                                   <  2:410080/150‑1 (MQ=255)
     gggTAATGAAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCCGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCa                                                                                                                                                  >  2:480655/1‑150 (MQ=255)
      ggTAATGAAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAg                                                                                                                                                 >  2:852960/1‑150 (MQ=255)
       gTAATGAAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAgc                                                                                                                                                <  1:521367/150‑1 (MQ=255)
        tAATGAAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAgcg                                                                                                                                               <  2:1418394/150‑1 (MQ=255)
        tAATGAAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAgcg                                                                                                                                               >  1:1459103/1‑150 (MQ=255)
         aaTGAAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTCTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGCAACTTCAgcaa                                                                                                                                              <  1:140090/150‑3 (MQ=255)
           tGAAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGCAACTTCAgcaaaa                                                                                                                                            >  1:1082173/1‑146 (MQ=255)
             aaGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCgcggcg                                                                                                                                          <  2:1221271/150‑1 (MQ=255)
             aaGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCgcggcg                                                                                                                                          >  2:83742/1‑150 (MQ=255)
             aaGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGCAACTTCAgcaaaaaa                                                                                                                                          <  2:1082173/150‑7 (MQ=255)
              aGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCgcggcgg                                                                                                                                         <  2:1083123/150‑1 (MQ=255)
              aGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCgcgg                                                                                                                                            <  1:1388776/147‑1 (MQ=255)
              aGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCgcgg                                                                                                                                            >  2:1388776/1‑147 (MQ=255)
                      atctcatctcCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGc                                                                                                                                 <  1:746209/150‑1 (MQ=255)
                      atctcatctcCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGc                                                                                                                                 <  1:1386947/150‑1 (MQ=255)
                      atctcatctcCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGc                                                                                                                                 <  1:1399453/150‑1 (MQ=255)
                             ctcCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCgg                                                                                                                          <  2:1272115/150‑1 (MQ=255)
                              tcCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGc                                                                                                                         <  2:1083817/150‑1 (MQ=255)
                              tcCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGc                                                                                                                         <  1:802830/150‑1 (MQ=255)
                               ccAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCa                                                                                                                        >  2:81173/1‑150 (MQ=255)
                                cAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCAc                                                                                                                       >  1:1264284/1‑150 (MQ=255)
                                  ggAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTc                                                                                                                     <  1:48326/150‑1 (MQ=255)
                                       cGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCtt                                                                                                                >  2:655366/1‑150 (MQ=255)
                                        ggCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCtta                                                                                                               <  2:1264284/150‑1 (MQ=255)
                                          cTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCTCGTCGGCGTCCTGCTCTTCGGCACTCGGCttatt                                                                                                             <  2:45540/150‑1 (MQ=255)
                                          cTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCGCTTCGGCACTCGGCttatt                                                                                                             <  1:292395/150‑1 (MQ=255)
                                           tGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCttatta                                                                                                            <  2:1053270/150‑1 (MQ=255)
                                              cagcagCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAAt                                                                                                         <  1:994786/150‑1 (MQ=255)
                                                gcagcGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTg                                                                                                       >  1:1191832/1‑150 (MQ=255)
                                                gcagcGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTg                                                                                                       <  1:388330/150‑1 (MQ=255)
                                                  agcGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCg                                                                                                     <  1:560078/150‑1 (MQ=255)
                                                  agcGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGATTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCg                                                                                                     <  2:1236058/150‑1 (MQ=255)
                                                        cgATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCa                                                                                               <  1:497246/150‑1 (MQ=255)
                                                           tgtgGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGa                                                                                            <  1:852960/150‑1 (MQ=255)
                                                                 gCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCa                                                                                      >  2:370008/1‑150 (MQ=255)
                                                                   aTCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGctc                                                                                                                               >  1:379026/1‑107 (MQ=255)
                                                                   aTCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGctc                                                                                                                               <  2:379026/107‑1 (MQ=255)
                                                                   aTCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCAcc                                                                                    <  1:188741/150‑1 (MQ=255)
                                                                     cAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCAt                                                                                  <  1:791895/150‑1 (MQ=255)
                                                                     cAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCAt                                                                                  >  2:108059/1‑150 (MQ=255)
                                                                       ggTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCt                                                                                >  1:970228/1‑150 (MQ=255)
                                                                        gTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTg                                                                               <  1:920403/150‑1 (MQ=255)
                                                                           gCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGAtt                                                                            <  2:1432903/150‑1 (MQ=255)
                                                                           cccAACACCTTCGACTACCACGCGACCGTTATCTTTACTCCGGATTTGCGGCGTGTAGCTTTTTAGCACCAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGAtt                                                                            >  2:1304142/2‑150 (MQ=255)
                                                                            ccAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTg                                                                           <  1:1268969/150‑1 (MQ=255)
                                                                                   cTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGa                                                                    >  2:678895/1‑150 (MQ=255)
                                                                                      cGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGatat                                                                 >  2:1153022/1‑150 (MQ=255)
                                                                                          tACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATcccc                                                             >  2:209626/1‑150 (MQ=255)
                                                                                           aCCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTcc                                                                                                                                    >  2:623907/1‑78 (MQ=255)
                                                                                           aCCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTcc                                                                                                                                    <  1:623907/78‑1 (MQ=255)
                                                                                           aCCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCg                                                            >  1:1143000/1‑150 (MQ=255)
                                                                                           aCCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCg                                                            >  1:547082/1‑150 (MQ=255)
                                                                                             cACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACGATCTGATTGAGACGGATATCCCCGcc                                                          >  1:1143374/1‑150 (MQ=255)
                                                                                             cACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGcc                                                          >  2:1213748/1‑150 (MQ=255)
                                                                                             cACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGcc                                                          >  1:626865/1‑150 (MQ=255)
                                                                                                gcgACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGccacc                                                       <  2:1406331/150‑1 (MQ=255)
                                                                                                gcgACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGccacc                                                       <  2:939057/150‑1 (MQ=255)
                                                                                                      gTTAGCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATcgc                                                 <  1:471453/150‑1 (MQ=255)
                                                                                                       ttATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATcgcc                                                <  1:1192069/150‑1 (MQ=255)
                                                                                                        tATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATcgccg                                               <  2:1402456/150‑1 (MQ=255)
                                                                                                        tATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATcgccg                                               <  1:307290/150‑1 (MQ=255)
                                                                                                        tATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATACCCGCCACCAATcgccc                                               >  2:677015/1‑149 (MQ=255)
                                                                                                          tCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCgccgcc                                             >  1:1168267/1‑150 (MQ=255)
                                                                                                          tCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCgccgcc                                             <  1:516143/150‑1 (MQ=255)
                                                                                                          tCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCgccgcc                                             <  2:626865/150‑1 (MQ=255)
                                                                                                             ttACTGCGGTATTGAGGCGTGTAGCTATTTTGCTGAAGTTGCTTCAGCGCGGCGGCCTCCTGCTCTTCGGCACTAGCATTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGc                                          >  2:85265/1‑150 (MQ=255)
                                                                                                               aCTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCaa                                        >  1:685056/1‑150 (MQ=255)
                                                                                                                 tGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACt                                      <  2:1191832/150‑1 (MQ=255)
                                                                                                                         tGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGTGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGAca                              <  2:259898/150‑1 (MQ=255)
                                                                                                                         tGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGTGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGAca                              >  1:259898/1‑150 (MQ=255)
                                                                                                                          gcggcgTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGAcat                             >  1:1199059/1‑150 (MQ=255)
                                                                                                                                tgtAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGa                       <  1:1088197/150‑1 (MQ=255)
                                                                                                                                 gtAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGaa                      >  1:290536/1‑150 (MQ=255)
                                                                                                                                  tAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAg                     >  2:877485/1‑150 (MQ=255)
                                                                                                                                   aGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGt                    >  2:1401585/1‑150 (MQ=255)
                                                                                                                                    gCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGtt                   >  2:1323034/1‑150 (MQ=255)
                                                                                                                                    gCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGtt                   >  1:1236698/1‑150 (MQ=255)
                                                                                                                                     ccgccgcgCTGAAGCAACTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTg                  >  2:230135/8‑150 (MQ=255)
                                                                                                                                     cTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTg                  <  1:316428/150‑1 (MQ=255)
                                                                                                                                       tttttGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTa                >  2:598377/1‑150 (MQ=255)
                                                                                                                                        ttttGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTaa               >  2:572587/1‑150 (MQ=255)
                                                                                                                                         cgcgCTGAAGAAACTACAGCGCGGCGCCGTCCTCCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCCGCAACTCGTCGACATCATTGAAGTTGTAAc              >  2:904434/4‑150 (MQ=255)
                                                                                                                                           tGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAAcgc            >  2:650181/1‑150 (MQ=255)
                                                                                                                                           tGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAAcgc            >  1:328108/1‑150 (MQ=255)
                                                                                                                                            gCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAAcgcg           <  1:480655/150‑1 (MQ=255)
                                                                                                                                             cTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCgg          >  2:98349/1‑150 (MQ=255)
                                                                                                                                                  gTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAACCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGgcagca     <  1:370008/150‑1 (MQ=255)
                                                                                                                                                    tGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACtcgtcg                                 >  2:82343/1‑120 (MQ=255)
                                                                                                                                                    tGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACtcgtcg                                 <  1:82343/120‑1 (MQ=255)
                                                                                                                                                     gCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGAt  <  2:1143000/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCCTGGGTAATGAAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGAT  >  NC_000913/2911768‑2912066

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