Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,177,340 A→C 6.1% intergenic (‑136/+172) gatY ← / ← fbaB D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,177,3400AC 6.1% 322.4 / 8.3 117intergenic (‑136/+172)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base A (55/54);  new base C (3/5);  total (58/59)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

Only 100 of 117 total aligned reads displayed.

CACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGG  >  NC_000913/2177196‑2177488
                                                                                                                                                |                                                                                                                                                    
 aCGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAAAGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAg                                                                                                                                                >  1:816180/1‑150 (MQ=255)
 aCGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAg                                                                                                                                                <  2:1323560/150‑1 (MQ=255)
   ttACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCTATCGGGAGGcc                                                                                                                                              <  1:1037405/149‑1 (MQ=255)
   gTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGcc                                                                                                                                              <  2:4621/150‑1 (MQ=255)
    tACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCg                                                                                                                                             <  1:153701/150‑1 (MQ=255)
    tACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATAAAATGAAATTATTTGTTTTTAAATAACGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGCGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCg                                                                                                                                             <  2:780027/150‑1 (MQ=255)
      cATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                                                                           >  2:1331694/1‑150 (MQ=255)
          ttCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAAGGCCTCCAGATAGGGAAGCCGTAGCaa                                                                                                                                       >  2:155499/1‑150 (MQ=255)
          ttCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATAACAAAAACGACAATATGAAAATTATTAGAGGATTGAAATGCAAAAAAAAGGCCTACCAATAGGGAATCCGTAGCaa                                                                                                                                       >  1:786989/1‑150 (MQ=255)
           tCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAg                                                                                                                                     <  1:533683/150‑1 (MQ=255)
           tCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATTAAAATTATTAGAGGAGTGAAAGGCAAAAAAACCGCCTACCGATAGGGAAGCCGTAGCaaa                                                                                                                                      >  2:1776587/1‑150 (MQ=255)
            cATATCCTGTCGTTTTTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtttttttg                                                                                                                                     <  2:1314963/150‑8 (MQ=255)
            cATATCCTGTCGTTTGTTTTGGATTTGAAAATATAATGAAATTATTTGTTTTTTAATATAGAGATAAAGATCACAAAAAAGACAAAATTAAAATTATTAGAGGAGTGAAAGGCAAAAAAACGGCCACCCGATATGGAAGCCGTAGCaaat                                                                                                                                     >  2:808292/1‑149 (MQ=255)
              tatCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGATGAGTGAAAGGCAAAAAAACGGCCTCCAGATAGGGAAGCCGTAGCAAAGTg                                                                                                                                   >  2:349468/1‑150 (MQ=255)
              tatCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTg                                                                                                                                   >  1:205761/1‑150 (MQ=255)
                tCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgcg                                                                                                                                 >  1:231293/1‑150 (MQ=255)
                  cTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGcgcg                                                                                                                               >  2:476655/1‑150 (MQ=255)
                   tGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTTAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgt                                                                                                                              >  2:702999/1‑150 (MQ=255)
                       gtttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCAGATAGGGAAGCCGTAGCAAAGTGCGCGTGttt                                                                                                                          >  1:978070/1‑150 (MQ=255)
                          tgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa                                                                                                                       >  1:1489057/1‑150 (MQ=255)
                            ttttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGCGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                                                                                                                     <  1:208431/150‑1 (MQ=255)
                                cGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCgg                                                                                                                 <  1:1556929/150‑1 (MQ=255)
                                cGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCgg                                                                                                                 <  2:924879/150‑1 (MQ=255)
                                  aTTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAAAGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAATGTGCGGGTGTTTTTATGCCGGAt                                                                                                               >  1:1737962/1‑150 (MQ=255)
                                     taaaaaTATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACATTATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCg                                                                                                            >  1:985299/3‑150 (MQ=255)
                                     tCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCg                                                                                                            >  2:1333262/1‑150 (MQ=255)
                                     tCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCg                                                                                                            >  2:892135/1‑150 (MQ=255)
                                     tCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCg                                                                                                            >  1:1281676/1‑150 (MQ=255)
                                     tCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCg                                                                                                            >  1:930718/1‑150 (MQ=255)
                                      cAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATACCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg                                                                                                           <  2:35879/150‑1 (MQ=255)
                                      cAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg                                                                                                           <  1:317606/150‑1 (MQ=255)
                                      cAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg                                                                                                           >  2:208714/1‑150 (MQ=255)
                                       aaaaTATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgt                                                                                                          <  1:1109877/150‑1 (MQ=255)
                                             taATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAAc                                                                                                    <  2:1734534/150‑1 (MQ=255)
                                             taATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCCCGTGTTTTTATGCCGGATGCGGTGTAAAc                                                                                                    <  1:1081705/150‑1 (MQ=255)
                                              aaTGAAATTATTTTTTTTTAAATCTCGAGATAACCAGCACACAAACCTCAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGTGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACa                                                                                                   <  1:81886/150‑1 (MQ=255)
                                                 gAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAc                                                                                                <  2:522134/150‑1 (MQ=255)
                                                    attattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTg                                                                                             <  1:929610/150‑1 (MQ=255)
                                                    attattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTg                                                                                             <  2:1453458/150‑1 (MQ=255)
                                                      tattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGc                                                                                           <  2:549750/150‑1 (MQ=255)
                                                          tGTTTTTAAATATCAAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTg                                                                                       >  2:1256177/1‑150 (MQ=255)
                                                           gTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGc                                                                                      >  2:814924/1‑150 (MQ=255)
                                                            tttttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                                                                                     <  2:1005133/150‑1 (MQ=255)
                                                            tttttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                                                                                     <  2:1456479/150‑1 (MQ=255)
                                                              tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTAt                                                                                   >  2:1602607/1‑150 (MQ=255)
                                                               ttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATg                                                                                  >  2:207367/1‑150 (MQ=255)
                                                               ttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTCCTATg                                                                                  >  2:774976/1‑150 (MQ=255)
                                                                tAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGt                                                                                 <  2:134758/150‑1 (MQ=255)
                                                                tAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCACGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGACGGATGCGGTGTAAACATACGTGGCAGTGCTATGt                                                                                 >  2:1797609/1‑150 (MQ=255)
                                                                 aattataGAGATAACTATCACAAAAACTAAAATATGAAAATTATTGGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGGTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTa                                                                                <  2:1169738/147‑1 (MQ=255)
                                                                    tatCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGc                                                                             <  2:419269/150‑1 (MQ=255)
                                                                       cGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATg                                                                          <  2:684862/150‑1 (MQ=255)
                                                                          taTAACGATCACAAAAACGACAATATGAAAAATATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa                                                                       <  2:1330717/149‑1 (MQ=255)
                                                                          gaTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa                                                                       <  2:1165487/150‑1 (MQ=255)
                                                                          gaTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa                                                                       <  2:1070927/150‑1 (MQ=255)
                                                                            tAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAg                                                                     >  2:889047/1‑150 (MQ=255)
                                                                             aaCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGa                                                                    >  1:1679267/1‑150 (MQ=255)
                                                                               cGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGCTAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACg                                                                  >  2:1675248/1‑150 (MQ=255)
                                                                               cGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACg                                                                  >  2:1869438/1‑150 (MQ=255)
                                                                                gATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGc                                                                 >  1:1611105/1‑150 (MQ=255)
                                                                                   cacaAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTc                                                              >  1:1652818/1‑150 (MQ=255)
                                                                                     caAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAg                                                            >  2:12438/1‑150 (MQ=255)
                                                                                      aaaaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGc                                                           >  2:1338469/1‑150 (MQ=255)
                                                                                       aaaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCg                                                          >  1:893776/1‑150 (MQ=255)
                                                                                          aCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCg                                                       >  2:745667/1‑150 (MQ=255)
                                                                                                atatGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAg                                                 <  1:992578/150‑1 (MQ=255)
                                                                                                 aatcaaaaTTATCCGACGAGGGAAAGGCAAAAAACCGGCCCCCCGATCGGGAAGCCGTAGCAACGGGCGCGGGTTTTTAGGCCGGATGCGGTGAAAACATACGGGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAgg                                                <  2:962502/146‑1 (MQ=255)
                                                                                                   tGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCa                                              <  1:1409637/150‑1 (MQ=255)
                                                                                                     aaaaTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATGCGACGCTGAAGCGTCTTATCATGCCTa                                            <  1:451059/150‑1 (MQ=255)
                                                                                                      aaaTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAg                                           <  2:610296/150‑1 (MQ=255)
                                                                                                             tCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCg                                    <  1:1447148/150‑1 (MQ=255)
                                                                                                              cGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCgg                                   <  2:1019982/150‑1 (MQ=255)
                                                                                                               gaggagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTTCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGa                                  <  1:1869007/150‑1 (MQ=255)
                                                                                                                aggagTGAAAGGCAAAAATACGGCCTCCCGATAGGGAAGCGGTAGCTTAGGGCGCGTGTTTTTATGCCGGTTGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGAt                                 <  2:1768355/150‑1 (MQ=39)
                                                                                                                aggagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGAt                                 <  1:476655/150‑1 (MQ=255)
                                                                                                                 ggagTGACAGTCAAAAAAACGTCCTCCCGATAGGTAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATTTAGGCATTATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATg                                <  1:1366371/150‑1 (MQ=255)
                                                                                                                 ggagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGAtt                                >  2:1546221/1‑149 (MQ=255)
                                                                                                                    gTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACTCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAg                             >  1:1616348/1‑150 (MQ=255)
                                                                                                                    gTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAg                             <  1:1256177/150‑1 (MQ=255)
                                                                                                                      gAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCCGGCATAGGTTGCCGGATGCAGCg                           >  2:1424475/1‑150 (MQ=255)
                                                                                                                      gAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGc                                               <  2:387526/130‑1 (MQ=255)
                                                                                                                      gAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGc                                               >  1:387526/1‑130 (MQ=255)
                                                                                                                      gAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAGGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGAGAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCg                           >  2:886298/1‑150 (MQ=255)
                                                                                                                      gAAAGGCAAAAAAACGGCCTCACGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGAGTACACATCCGGGGCAGTGCTATGTAGGCATGATACGACGCTGCAGCGTCGCATCAGGCTTAGGTTGCCGGATGCAGCg                           >  1:116714/1‑150 (MQ=255)
                                                                                                                        aaGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTa                         <  2:231293/150‑1 (MQ=255)
                                                                                                                         aGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTaa                        >  2:520699/1‑150 (MQ=255)
                                                                                                                            cAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACg                     <  2:816180/150‑1 (MQ=255)
                                                                                                                            cAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACg                     >  1:1369205/1‑150 (MQ=255)
                                                                                                                                 ctaCGGCTTCCCTATCGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTa                >  1:747496/3‑150 (MQ=255)
                                                                                                                                  aaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTAt               <  2:1281676/150‑1 (MQ=255)
                                                                                                                                    cGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATcc             <  1:23404/150‑1 (MQ=255)
                                                                                                                                     ggCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCg            <  1:666785/150‑1 (MQ=255)
                                                                                                                                     ccccacccGATCCGGAAGCCGTAGCAACGTGCGCGCGTCTTTCTGCCGGATGCGGTGTAAACACACCTGGCAGTGCTATGCAGGCATGATACGACGCATCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCg            <  2:197702/145‑1 (MQ=37)
                                                                                                                                      gCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGt           >  2:1469529/1‑150 (MQ=255)
                                                                                                                                      gCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCGTATCCGt           >  1:767222/1‑150 (MQ=255)
                                                                                                                                       ccTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTc          <  2:1042167/150‑1 (MQ=255)
                                                                                                                                       ccTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTc          <  2:205761/150‑1 (MQ=255)
                                                                                                                                           ccGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTAc      <  2:1407589/150‑1 (MQ=255)
                                                                                                                                            cGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACg     >  1:79044/1‑150 (MQ=255)
                                                                                                                                               tAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACgggg  <  1:950718/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGG  >  NC_000913/2177196‑2177488

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: