Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,687,043 G→T 10.7% intergenic (+18/+27) hmp → / ← glnB nitric oxide dioxygenase/nitrogen regulatory protein PII‑1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,687,0430GT10.7% 293.4 / 17.4 93intergenic (+18/+27)hmp/glnBnitric oxide dioxygenase/nitrogen regulatory protein PII‑1
Reads supporting (aligned to +/- strand):  ref base G (41/42);  new base T (5/5);  total (46/47)
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.19e-01

CCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACA  >  NC_000913/2686897‑2687192
                                                                                                                                                  |                                                                                                                                                     
ccGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGaaa                                                                                                                                                    <  1:18955/150‑1 (MQ=255)
 cGACAATGCAGTTCTATTTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAAc                                                                                                                                                   >  1:1469451/1‑150 (MQ=255)
  gACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACa                                                                                                                                                  >  2:528998/1‑150 (MQ=255)
   aCAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAAAATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAAAAt                                                                                                                                                 >  2:731908/1‑150 (MQ=255)
     aaTGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATcc                                                                                                                                               <  2:1462927/150‑1 (MQ=255)
        gCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGc                                                                                                                                            <  2:1417590/150‑1 (MQ=255)
        gCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGc                                                                                                                                            <  1:320424/150‑1 (MQ=255)
         cAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGTTTCCATCCGGCa                                                                                                                                           <  1:1406801/150‑1 (MQ=255)
         cAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGGGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGTTTCCATCCGGCa                                                                                                                                           <  2:130036/150‑1 (MQ=255)
           gTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGGGCTGTAATTTGATGTTGCCGGATGTTTCCATCCGGCAAc                                                                                                                                         <  1:842393/150‑1 (MQ=255)
             tCTATCTCTGCGGCCCCGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAAccc                                                                                                                                       <  2:1458800/150‑1 (MQ=255)
                aTCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTg                                                                                                                                    >  2:940863/1‑150 (MQ=255)
                     tGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGAcgcg                                                                                                                               >  1:1274958/1‑150 (MQ=255)
                     tGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGAcgcg                                                                                                                               >  2:319653/1‑150 (MQ=255)
                      gCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCgg                                                                                                                              <  1:678536/150‑1 (MQ=255)
                        ggCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGtt                                                                                                                            >  1:215217/1‑150 (MQ=255)
                         gCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGttt                                                                                                                           >  1:1357963/1‑150 (MQ=255)
                           ccGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTaa                                                                                                                         >  2:430596/1‑150 (MQ=255)
                             ggTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAAt                                                                                                                       <  1:319653/150‑1 (MQ=255)
                                  gCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCcg                                                                                                                  >  1:798782/1‑150 (MQ=255)
                                  gCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCAGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCcg                                                                                                                  >  1:1485593/1‑150 (MQ=255)
                                   cTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCcgc                                                                                                                 >  1:1169223/1‑150 (MQ=255)
                                      cATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGcgtc                                                                                                              <  1:110161/150‑1 (MQ=255)
                                          cAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTcc                                                                                                          <  1:1416442/150‑1 (MQ=255)
                                          cAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTcc                                                                                                          <  2:1480230/150‑1 (MQ=255)
                                             tttACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCctct                                                                                                       <  2:82988/150‑1 (MQ=255)
                                               tACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTc                                                                                                     >  2:300686/1‑150 (MQ=255)
                                                 ccGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTa                                                                                                                          >  2:1167445/1‑127 (MQ=255)
                                                 ccGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTa                                                                                                                          <  1:1167445/127‑1 (MQ=255)
                                                  cgcgAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCAcc                                                                                                  <  1:228407/150‑1 (MQ=255)
                                                  cgcgAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCAcc                                                                                                  <  1:528998/150‑1 (MQ=255)
                                                  cgcgAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCAcc                                                                                                  <  1:940863/150‑1 (MQ=255)
                                                   gcgAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGTTTCCATCCGGc                                                                                                                                            <  2:1622152/107‑1 (MQ=255)
                                                   gcgAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGTTTCCATCCGGc                                                                                                                                            >  1:1622152/1‑107 (MQ=255)
                                                    cgAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCgg                                                                                                <  1:1219236/150‑1 (MQ=255)
                                                          aGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGa                                                                                          <  2:1448477/150‑1 (MQ=255)
                                                             tAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGc                                                                                       <  1:1046650/150‑1 (MQ=255)
                                                              aGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCg                                                                                      <  2:1437347/150‑1 (MQ=255)
                                                               gTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGa                                                                                     >  2:1112090/1‑150 (MQ=255)
                                                                tGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGaa                                                                                    <  1:883935/150‑1 (MQ=255)
                                                                 ggATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAAt                                                                                   <  2:1169223/150‑1 (MQ=255)
                                                                  gATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATg                                                                                  >  1:818660/1‑150 (MQ=255)
                                                                      tGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGAccc                                                                              >  2:1331026/1‑150 (MQ=255)
                                                                       gggCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCg                                                                             >  2:27458/1‑150 (MQ=255)
                                                                         gCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTg                                                                           <  2:798782/150‑1 (MQ=255)
                                                                              aagaagGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACg                                                                      <  1:363913/150‑1 (MQ=255)
                                                                               aGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTg                                                                           <  2:1619800/144‑1 (MQ=255)
                                                                               aGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTg                                                                           >  1:1619800/1‑144 (MQ=255)
                                                                               aGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGt                                                                     >  1:1029381/1‑150 (MQ=255)
                                                                                 ggccaaagCATTCGTAATGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCa                                                                   >  1:1516170/5‑150 (MQ=255)
                                                                                 cAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCa                                                                   >  1:1337817/1‑150 (MQ=255)
                                                                                  gccaaagCATTCGTAATGAATGTTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCaa                                                                  >  2:252162/4‑150 (MQ=255)
                                                                                      aaaCATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGAc                                                              <  1:1030939/150‑1 (MQ=255)
                                                                                       aaCATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACg                                                             <  2:567264/150‑1 (MQ=255)
                                                                                        aCATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGa                                                            <  2:1274958/150‑1 (MQ=255)
                                                                                         cattcattACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTTGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGaa                                                           >  2:453058/1‑150 (MQ=255)
                                                                                              attACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGAttt                                                      >  1:356026/1‑150 (MQ=255)
                                                                                                tACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTa                                                    >  2:666454/1‑150 (MQ=255)
                                                                                                  cGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTAcc                                                  >  1:106051/1‑150 (MQ=255)
                                                                                                   gAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCg                                                 >  1:1595398/1‑150 (MQ=255)
                                                                                                   gAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCg                                                 >  1:351812/1‑150 (MQ=255)
                                                                                                   gAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCg                                                 <  2:809823/150‑1 (MQ=255)
                                                                                                      tGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCa                                              <  1:430596/150‑1 (MQ=255)
                                                                                                       gCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCAc                                             >  1:471354/1‑150 (MQ=255)
                                                                                                        cTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCAcc                                            >  1:1343064/1‑150 (MQ=255)
                                                                                                             gCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGAttt                                       <  1:707491/150‑1 (MQ=255)
                                                                                                                cGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGc                                    <  2:772368/150‑1 (MQ=255)
                                                                                                                cGCATAAGGTGCTGTAATTTGATGTTACCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGc                                    >  1:726042/1‑150 (MQ=255)
                                                                                                                 gCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGcc                                   <  1:1331026/150‑1 (MQ=255)
                                                                                                                  cATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCg                                  <  1:1112090/150‑1 (MQ=255)
                                                                                                                    tAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGt                                >  2:143134/1‑150 (MQ=255)
                                                                                                                      aGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGttt                              >  1:1575068/1‑150 (MQ=255)
                                                                                                                        gTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTgc                            >  2:637205/1‑150 (MQ=255)
                                                                                                                           cTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGcc                         >  2:430863/1‑150 (MQ=255)
                                                                                                                               aaTTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGc                     <  2:818660/150‑1 (MQ=255)
                                                                                                                                  ttGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGaa                  >  2:432500/1‑150 (MQ=255)
                                                                                                                                   tGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGaat                 <  1:1329562/150‑1 (MQ=255)
                                                                                                                                      tGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGaataat              >  1:663338/1‑150 (MQ=255)
                                                                                                                                      ggTTGCCGGATGTTTCCATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGc                                    <  2:1420758/127‑1 (MQ=255)
                                                                                                                                      ggTTGCCGGATGTTTCCATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGc                                    >  1:1420758/2‑128 (MQ=255)
                                                                                                                                       gTTGCCGGATGTTTCCATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATg             >  2:992101/1‑150 (MQ=255)
                                                                                                                                       gTTGCCGGATGTTTCCATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATg             >  1:204831/1‑150 (MQ=255)
                                                                                                                                         tGCCGGATGTTTCCATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGt           >  1:1222453/1‑150 (MQ=255)
                                                                                                                                          gCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGACTGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTa          >  1:595583/1‑150 (MQ=255)
                                                                                                                                          gCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGa                                          <  2:1523617/118‑1 (MQ=255)
                                                                                                                                          gCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGa                                          >  1:1523617/1‑118 (MQ=255)
                                                                                                                                          gCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTa          <  2:1078654/150‑1 (MQ=255)
                                                                                                                                          gCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTa          <  1:432500/150‑1 (MQ=255)
                                                                                                                                            cGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATc        <  1:696584/150‑1 (MQ=255)
                                                                                                                                            cGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATc        <  1:430863/150‑1 (MQ=255)
                                                                                                                                             ggATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCg       <  1:27504/150‑1 (MQ=255)
                                                                                                                                               aTGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGac     >  1:64597/1‑150 (MQ=255)
                                                                                                                                               aTGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGac     <  1:143134/150‑1 (MQ=255)
                                                                                                                                                  gAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATg             <  2:1256276/139‑1 (MQ=255)
                                                                                                                                                  gAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGAcaca  <  2:1357963/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCATTACGAATGCTTTGGCCCGCATAAGGTGCTGTAATTTGATGTTGCCGGATGGAAACATCCGGCAACCCTTGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACA  >  NC_000913/2686897‑2687192

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: