Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,081,058 C→T 5.8% V69M (GTG→ATG) ‡ sbmC ← DNA gyrase inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,081,0580CT 5.8% 277.0 / 4.8 89V69M (GTG→ATG) ‡sbmCDNA gyrase inhibitor
Reads supporting (aligned to +/- strand):  ref base C (41/41);  new base T (2/3);  total (44/45)
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.14e-01

AGACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCACGGGCCAA  >  NC_000913/2080913‑2081207
                                                                                                                                                 |                                                                                                                                                     
agaCTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCacgacg                                                                                                                                                   <  1:508055/150‑1 (MQ=255)
  aCTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGgt                                                                                                                                                 >  2:1307248/1‑150 (MQ=255)
  aCTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGgt                                                                                                                                                 <  2:537696/150‑1 (MQ=255)
   cTATTAAAGAACTGATACCAGGGTTTAGAAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGCGATGCGTAATGTAACCCGGAACCGTCACGGCGgtg                                                                                                                                                >  2:61609/1‑150 (MQ=255)
     aatAAATAACAGATACAAGGGTTTAGCAAAATCATCAACGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTc                                                                                                                                              <  1:36837/148‑1 (MQ=255)
      ttAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCg                                                                                                                                             <  2:737931/150‑1 (MQ=255)
      ttAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGACACCGTCGTGACGGTGCCg                                                                                                                                             <  2:293083/150‑1 (MQ=255)
       tAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCGTGACGGTGCCgg                                                                                                                                            <  1:998924/150‑2 (MQ=255)
       tAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACAACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGc                                                                                                                                            >  1:253941/1‑150 (MQ=255)
         aaGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAg                                                                                                                                          <  1:717513/150‑1 (MQ=255)
         aaGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGACAGAATGCCGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAg                                                                                                                                          >  1:927005/1‑150 (MQ=255)
         aaGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTGCTGTCAGACTGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCGTGACGGTGCCgggt                                                                                                                                          <  1:78671/150‑4 (MQ=255)
               tGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCTTACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAAt                                                                                                                                    >  2:903654/1‑150 (MQ=255)
               tGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAAt                                                                                                                                    <  1:270707/150‑1 (MQ=255)
               tGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAAt                                                                                                                                    <  2:565171/150‑1 (MQ=255)
                    ccAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGTGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTc                                                                                                                               <  1:537010/150‑1 (MQ=255)
                    ccAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGCCGGTGTCGCAGCGTAATTTTTc                                                                                                                               >  1:840318/1‑150 (MQ=255)
                     cAGGGTTTAGCAAAATCATCACCGACTACACGCGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTATTTTTTcg                                                                                                                              >  2:1248126/1‑150 (MQ=255)
                     cAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTcg                                                                                                                              >  1:453244/1‑150 (MQ=255)
                      cgggTTTAGCAAAATCATCACCCACTACACGCGCAACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTcgg                                                                                                                             <  2:819840/149‑1 (MQ=255)
                               cAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGttt                                                                                                                    >  1:370738/1‑150 (MQ=255)
                                 aaaTCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCa                                                                                                                  >  1:850876/1‑150 (MQ=255)
                                  aaTCATCAGCGACTGCACGAGCTACCGCCACCGCATACTGACCTGCTGTAATTTCTGTGAGAATGACGCCCGGACTGTTTTCGGGAAGCGTAAAGTACCCCGGCGCCGTCAGGGCGGTGTCGCAGGGTAATTTTTCGGGGGgtggtgcag                                                                                                                 >  2:1163869/1‑144 (MQ=255)
                                   atcatcACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATc                                                                                                                <  1:1217582/150‑1 (MQ=255)
                                     catcaCCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTg                                                                                                              >  2:368244/1‑150 (MQ=255)
                                        caCCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGAt                                                                                                           <  1:1111801/150‑1 (MQ=255)
                                           cGACTACACGAGCTACCGCCACCGCATACTGACCCCCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGt                                                                                                        <  1:858904/150‑1 (MQ=255)
                                             aCTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCg                                                                                                      <  2:1003991/150‑1 (MQ=255)
                                             aCTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCg                                                                                                      <  2:977209/150‑1 (MQ=255)
                                              cTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGt                                                                                                     <  2:357184/150‑1 (MQ=255)
                                              cTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGt                                                                                                     <  1:279383/150‑1 (MQ=255)
                                                 cacGAGCTACCGCAACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGATGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAAt                                                                                                  >  2:552753/1‑150 (MQ=255)
                                                  acGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATa                                                                                                 >  1:1000883/1‑150 (MQ=255)
                                                   cGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAg                                                                                                >  1:812969/1‑150 (MQ=255)
                                                   cGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAg                                                                                                >  1:1248218/1‑150 (MQ=255)
                                                     aGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGTCGGGTGTTTCATCTGGATTGTCGTAATAGAc                                                                                              <  2:890560/150‑1 (MQ=255)
                                                     aGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGAc                                                                                              <  1:739343/150‑1 (MQ=255)
                                                     aGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGAc                                                                                              <  1:748500/150‑1 (MQ=255)
                                                     aGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGAc                                                                                              <  2:748479/150‑1 (MQ=255)
                                                      gCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACa                                                                                             >  2:416888/1‑150 (MQ=255)
                                                         aCCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGcc                                                                                          >  1:868127/1‑149 (MQ=255)
                                                            gCCCCCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGGCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCCAcc                                                                                       >  2:703330/1‑150 (MQ=255)
                                                                ccGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACt                                                                                   <  2:38911/150‑1 (MQ=255)
                                                                  gCATACTGACCACCTGTACTTTCTGTCAGAAAGACGCCCTCACTGTTTTCGGGAAGCGTAACGTACCCGGGAACCGTCCCGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTAATCTGGATTGTCGTAATAGACAGCAACCCACTcc                                                                                 <  2:70925/150‑1 (MQ=255)
                                                                     tACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTc                                                                              <  1:472368/150‑1 (MQ=255)
                                                                     aaCTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTc                                                                              <  1:541153/149‑1 (MQ=255)
                                                                       cTGACCACCTGAAATTTCTGTCAGAATGACGCCCTCACTGTTTTAGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCgg                                                                            <  2:834751/150‑1 (MQ=255)
                                                                         gcccaccTGTAATTTCTGTCAGAATGACGCCCGCAGTGGTGTCGGGAAGCGTAAAGTAACCCGGCCCCGTCACGAATGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGca                                                                          <  2:206051/148‑1 (MQ=255)
                                                                         gcacaccTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGCTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCCCTCCTTCGGca                                                                          >  1:820294/4‑150 (MQ=255)
                                                                            caccTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGGAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACaa                                                                       <  2:542214/150‑1 (MQ=255)
                                                                              ccTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAata                                                                     >  1:631772/1‑150 (MQ=255)
                                                                              ccTGTAATTTCTGTCAGAATGACGCCCTCACTGTTGTCGGGCAGCGGATGGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTGATCTGGATTGTCGGAATAGACCGCAACCCAATCCTTCGGCAGAatt                                                                     >  2:665659/1‑149 (MQ=255)
                                                                               cTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCATCCCACTCGTTCGGCACAatat                                                                    >  1:786397/1‑150 (MQ=255)
                                                                                tGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATAtt                                                                   >  2:804427/1‑150 (MQ=255)
                                                                                 gTAATTTCTGTCAGAATGACGCCGTCAGTGTTTTCGGTAAGAGTAAAGTAACTAGGCACAGTCAAGAAGGTGTCGCAGAGTAATTTTTAGGCGTGTGTTTAATCTGGTTTTTCGTAATAGTCAGCAACCCACTCCTTCGGCACAATAttt                                                                  <  1:61609/150‑1 (MQ=255)
                                                                                    ttttCTGTCAGAATGAGGCCCTCACTGTTTTGGGTATGCTTAATGTAACCCTGCACCGTCTCGACGGTGTCGAAGCGTGATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGTCAGCAACCCACTCCTTCGGCACAATATTTTTg                                                               <  2:249837/149‑1 (MQ=255)
                                                                                    aTTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTg                                                               >  2:232522/1‑150 (MQ=255)
                                                                                      ttCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGCCGGGGTCGCCGCGTAATTTTTCGGCTGGTGTTTCATCTGTGTTGTGGTATTAGCCAGCAACCCAATACTTCGGCACAATTTTTTTGCt                                                             >  2:823039/1‑150 (MQ=255)
                                                                                        cTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCAACTGGATTGTCGTAATAGACAGCAAACCACTCCTACGGCACAATATTTTTGCTAt                                                           >  1:26416/1‑150 (MQ=255)
                                                                                         tGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATc                                                          <  2:338102/150‑1 (MQ=255)
                                                                                          gTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTACCCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCACACCACTCCGTCGGCACAATATTTTTGCTATcg                                                         >  2:186609/1‑149 (MQ=255)
                                                                                          gTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCt                                                         >  2:530339/1‑150 (MQ=255)
                                                                                           tCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTTATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTa                                                        >  2:322876/1‑150 (MQ=255)
                                                                                               aaTGACGCCCGCACTGTTTGCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTATTTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCca                                                    <  2:507147/150‑1 (MQ=255)
                                                                                                 tGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCAACGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCGCTCCTTCGGCACACTATTTTTGCTATCTACCctca                                                  >  1:820262/1‑150 (MQ=255)
                                                                                                  gACGCCCTCACTGTTTTCGGGAAGCGTATAGTAACCCTGCACAGTCACGACGGTGTCGCAGCGTAATTGTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATTTTTTTGCTATATAACCCcat                                                 >  2:1192977/1‑150 (MQ=255)
                                                                                                    cGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCAcatca                                               <  1:1137576/150‑1 (MQ=255)
                                                                                                     gCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCCCTCCTTCGGCACAATATTTTTGCTATCTACCCAcatcat                                              >  1:1199450/1‑150 (MQ=255)
                                                                                                      cccTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCCCGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCCCTCCTTCGGCACAATATTTTTGCTATCTACCCACatcatc                                             >  1:68822/1‑150 (MQ=255)
                                                                                                       ccTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACAtcatca                                            <  2:511060/150‑1 (MQ=255)
                                                                                                         tCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAAc                                          <  1:1112654/150‑1 (MQ=255)
                                                                                                          cACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACt                                         >  1:1230790/1‑150 (MQ=255)
                                                                                                            cTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGc                                       >  2:417968/1‑150 (MQ=255)
                                                                                                             tGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCt                                      >  1:1160548/1‑150 (MQ=255)
                                                                                                                    gggAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGcc                               <  1:1307248/150‑1 (MQ=255)
                                                                                                                    gggAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGGTTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGcc                               <  1:903654/150‑1 (MQ=255)
                                                                                                                      gAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCctt                             <  1:1073467/150‑1 (MQ=255)
                                                                                                                      gAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCctt                             <  2:894405/150‑1 (MQ=255)
                                                                                                                      gAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCctt                             >  2:656978/1‑150 (MQ=255)
                                                                                                                         gCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCctttct                          <  1:368244/150‑1 (MQ=255)
                                                                                                                              aaGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTCCCCACATCAGCACCTGCTCAAAGCCGTGCTTTccc                     >  2:1231118/1‑150 (MQ=255)
                                                                                                                               aGTAACCCGGCACCGTCCCGGCGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTCTCTACCCACATCATCAGCTTCTCAAAGCCTTTCTTTTCCg                    >  2:213443/1‑150 (MQ=255)
                                                                                                                               aGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCg                    <  1:909833/150‑1 (MQ=255)
                                                                                                                                 tAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTc                  >  1:608253/1‑150 (MQ=255)
                                                                                                                                   ctgcgacACCGTCGTGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTg                >  1:560241/7‑150 (MQ=255)
                                                                                                                                   ctgcgacACCGTCGTGACGGTGTCGAAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCAACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATAAACTGCTCAAAGCCTTTCGTTACCGTATg                >  1:1175239/7‑150 (MQ=255)
                                                                                                                                    cccGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGt               >  1:352466/1‑150 (MQ=255)
                                                                                                                                      cGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTc             <  1:163199/150‑1 (MQ=255)
                                                                                                                                             gTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCACggg      <  2:840318/150‑1 (MQ=255)
                                                                                                                                               cACGACGGTGTCGCAGCGTAATTTTTAGGAGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCACCTGCTCAAAGTCTTTCTTTACCGTCTGTTCCCACGGGcc    <  1:186609/150‑1 (MQ=255)
                                                                                                                                                 cgacgGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGATATATACCCACATCATCAACTGCTCAAAGCCTTTCTTGACCGTCTGTTCCCACGGGCCaa  >  2:374455/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
AGACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCACGGGCCAA  >  NC_000913/2080913‑2081207

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: ‑