Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,575,222 A→G 5.7% intergenic (‑20/+25) pqqL ← / ← yddB putative zinc peptidase/putative TonB‑dependent receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,575,2220AG 5.7% 381.9 / 9.2 123intergenic (‑20/+25)pqqL/yddBputative zinc peptidase/putative TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base A (65/51);  new base G (2/5);  total (67/56)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 123 total aligned reads displayed.

GCCATTGTCCAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAATATCGGCACTGATCGTC  >  NC_000913/1575074‑1575370
                                                                                                                                                    |                                                                                                                                                    
gCCATTGTCCAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTAc                                                                                                                                                     <  1:1179958/150‑1 (MQ=255)
  cATTGTCCAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTAccc                                                                                                                                                   >  2:1356569/1‑150 (MQ=255)
  cATTGTCCAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTAccc                                                                                                                                                   <  1:1218694/150‑1 (MQ=255)
   aTTGTCCAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCg                                                                                                                                                  >  1:543006/1‑150 (MQ=255)
   aTTGTCCAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCg                                                                                                                                                  <  2:564783/150‑1 (MQ=255)
     tGTCCAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCTGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGTAAATTCGCCCGGt                                                                                                                                                <  2:141303/150‑1 (MQ=255)
      gTCCAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGtt                                                                                                                                               >  1:513021/1‑150 (MQ=255)
         cAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAg                                                                                                                                            <  2:633869/150‑1 (MQ=255)
          aGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGTAAATTCGCCCGGttttcg                                                                                                                                           <  2:1454230/150‑5 (MQ=255)
           gTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGTAAATTCGCCCGGttttcgc                                                                                                                                          <  2:662133/150‑6 (MQ=255)
              gCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTaaaa                                                                                                                                       >  2:727200/1‑150 (MQ=255)
               cccGGTAATAAACTTTTCATACTGAGGTAAGGCGGCTACAATAAGCAGCACGGTGATCAAAAAAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAAt                                                                                                                                      >  2:1598946/1‑150 (MQ=255)
                  ggTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCg                                                                                                                                   >  1:337150/1‑150 (MQ=255)
                     ttcaggtttttCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAAc                                                                                                                                <  1:1207657/143‑1 (MQ=255)
                     aaTTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAAc                                                                                                                                <  1:1115773/150‑1 (MQ=255)
                          aCTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTa                                                                                                                           >  1:865028/1‑150 (MQ=255)
                           cTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTaa                                                                                                                          >  1:221224/1‑150 (MQ=255)
                           cTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTaa                                                                                                                          >  1:737365/1‑150 (MQ=255)
                            ttttCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTaaa                                                                                                                         >  2:26735/1‑150 (MQ=255)
                            ttttCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTaaa                                                                                                                         <  2:438709/150‑1 (MQ=255)
                            ttttCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTaaa                                                                                                                         >  1:1031364/1‑150 (MQ=255)
                             tttCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTaaaa                                                                                                                        >  2:1012433/1‑150 (MQ=255)
                                cATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGCCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAAttt                                                                                                                     >  2:1162307/1‑150 (MQ=255)
                                  tCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCa                                                                                                                   <  1:258421/150‑1 (MQ=255)
                                  tCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCa                                                                                                                   >  2:1258457/1‑150 (MQ=255)
                                    cTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCAGGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATg                                                                                                                 >  1:747450/1‑150 (MQ=255)
                                         gTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGAc                                                                                                            <  2:1623314/150‑1 (MQ=255)
                                           aaGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACAt                                                                                                          <  1:392291/150‑1 (MQ=255)
                                                  cggcAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCa                                                                                                   >  1:915178/1‑150 (MQ=255)
                                                   ggcAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAg                                                                                                  <  2:678036/150‑1 (MQ=255)
                                                   ggcAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAg                                                                                                  <  1:1012433/150‑1 (MQ=255)
                                                      aaTCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGaaa                                                                                               >  2:1087494/1‑150 (MQ=255)
                                                       aTCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGATTTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAg                                                                                              >  2:417743/1‑150 (MQ=255)
                                                       aTCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGGACTATTAAAATTTCATGCTGACATCAAGCCAGAAAg                                                                                              >  1:1312265/1‑150 (MQ=255)
                                                          aGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTAc                                                                                           >  2:175449/1‑150 (MQ=255)
                                                            ccgccCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACgg                                                                                         <  1:1605571/150‑1 (MQ=255)
                                                             cgccCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGc                                                                                        <  1:1162307/150‑1 (MQ=255)
                                                             cgccCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGc                                                                                        >  2:1234519/1‑150 (MQ=255)
                                                              gccCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCAAAATTATCTCCATGCGAAAACCGGTCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGcc                                                                                       <  2:1584932/150‑1 (MQ=255)
                                                                ccGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCAc                                                                                     <  1:1350910/150‑1 (MQ=255)
                                                                  ggggAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTg                                                                                   <  1:263813/150‑1 (MQ=255)
                                                                     gAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCg                                                                                >  1:588410/1‑150 (MQ=255)
                                                                     gAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCg                                                                                <  1:943222/150‑1 (MQ=255)
                                                                      aGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGt                                                                               <  1:728278/150‑1 (MQ=255)
                                                                         aaCAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAgg                                                                            >  2:1604420/1‑150 (MQ=255)
                                                                         aaCAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAgg                                                                            >  2:179235/1‑150 (MQ=255)
                                                                            agagTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATCAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCg                                                                         <  1:1506464/150‑1 (MQ=255)
                                                                              agTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGcc                                                                       <  1:859191/150‑1 (MQ=255)
                                                                                  gCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCacac                                                                   >  2:376207/1‑150 (MQ=255)
                                                                                    cACTAACGACAGTAAGAAACAGAGGTTTCTCATAATTATCTCCAAGCGAAAACCCGCCGAATTTAACAGGTTAACAAAAATCCGAACTATTAAAAATTCATACTGACATCAACCCAGAAAGTACGGCAACTGGCGTAGGTCGCCacacac                                                                 >  2:1479613/1‑148 (MQ=255)
                                                                                     aCTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCgg                                                                >  2:1519345/1‑150 (MQ=255)
                                                                                        aaCGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGgtgt                                                             >  2:1133423/1‑150 (MQ=255)
                                                                                        aaCGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGgtgt                                                             >  1:685483/1‑150 (MQ=255)
                                                                                           gTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCCCACCGGTGTTCg                                                          >  2:275103/1‑150 (MQ=255)
                                                                                           gTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCg                                                          >  2:118761/1‑150 (MQ=255)
                                                                                            tCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGt                                                         >  1:1051224/1‑150 (MQ=255)
                                                                                              aGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTg                                                       <  1:1408722/150‑1 (MQ=255)
                                                                                               gTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGt                                                      >  1:1423996/1‑150 (MQ=255)
                                                                                               gTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGt                                                      >  1:1580042/1‑150 (MQ=255)
                                                                                               gTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGt                                                      <  1:848412/150‑1 (MQ=255)
                                                                                                tAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTa                                                     <  2:337150/150‑1 (MQ=255)
                                                                                                tAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTa                                                     <  2:132230/150‑1 (MQ=255)
                                                                                                  aGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATc                                                   <  2:1600786/150‑1 (MQ=255)
                                                                                                  aGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATc                                                   <  2:737365/150‑1 (MQ=255)
                                                                                                  aGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATc                                                   <  1:1258457/150‑1 (MQ=255)
                                                                                                   gAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCa                                                  >  2:686264/1‑150 (MQ=255)
                                                                                                     aaCAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAAc                                                >  1:1222498/1‑150 (MQ=255)
                                                                                                         gagGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCg                                            >  2:1390761/1‑150 (MQ=255)
                                                                                                         gagGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCg                                            >  1:511140/1‑150 (MQ=255)
                                                                                                          agGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCgg                                           <  1:995260/150‑1 (MQ=255)
                                                                                                               tctcATAATTATCTCCATGCGAAAACCGGGCGAATTTAACCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTg                                      <  1:1479613/150‑1 (MQ=255)
                                                                                                                ctcATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCTAACTATTAAAATTTCATGCTGACATCAAGCCAGACAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGc                                     <  2:1619763/150‑1 (MQ=255)
                                                                                                                ctcATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTTTTCGTTGTATCAACATCGGTTTTGc                                     <  2:747450/150‑1 (MQ=255)
                                                                                                                ctcATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTGGTTGTATCAACAGCGGTTTTGc                                     <  2:887952/150‑1 (MQ=255)
                                                                                                                 tcATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCt                                    <  1:289628/150‑1 (MQ=255)
                                                                                                                  cATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTa                                   >  1:136175/1‑150 (MQ=255)
                                                                                                                  cATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTa                                   >  2:1300962/1‑150 (MQ=255)
                                                                                                                  cATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTa                                   >  2:554120/1‑150 (MQ=255)
                                                                                                                  cATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCAACACCGGTGTTCGTTGTATCAACAGCGATTTTGCTa                                   >  1:668386/1‑150 (MQ=255)
                                                                                                                    tAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATc                                 >  2:1068449/1‑150 (MQ=255)
                                                                                                                     aaTTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATcc                                <  1:641117/150‑1 (MQ=255)
                                                                                                                        tATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGt                             <  1:96966/150‑1 (MQ=255)
                                                                                                                            tcCATGCGAAATCCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCAACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGGTATCCAGTACAt                         <  2:543006/150‑1 (MQ=255)
                                                                                                                             ccATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACAtt                        >  1:1301875/1‑150 (MQ=255)
                                                                                                                               aTGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTga                      <  1:1598946/150‑1 (MQ=255)
                                                                                                                                  cGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGaa                   >  2:254504/1‑150 (MQ=255)
                                                                                                                                   gAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAat                  >  1:734355/1‑150 (MQ=255)
                                                                                                                                   cttaaCCGGGTAAATTCGCCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTg                                                       >  2:216331/4‑113 (MQ=255)
                                                                                                                                   cttaaCCGGGTAAATTCGCCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTg                                                       <  1:216331/110‑1 (MQ=255)
                                                                                                                                       aCCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAATATCg              <  2:583036/150‑1 (MQ=255)
                                                                                                                                       aCCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAATATCg              <  1:211290/150‑1 (MQ=255)
                                                                                                                                        ccGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAATATCgg             >  1:760206/1‑150 (MQ=255)
                                                                                                                                        ccGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTGTTGGTATCCAGTACATTGAGAATATCgg             >  2:1612715/1‑150 (MQ=255)
                                                                                                                                           ggCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAATATCGGCAc          >  1:364661/1‑150 (MQ=255)
                                                                                                                                            gCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAATATCGGCACt         >  1:605091/1‑150 (MQ=255)
                                                                                                                                             cGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAATATCGGCACTg        >  2:573434/1‑150 (MQ=255)
                                                                                                                                               aaTTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAATATCGGCACTGAt      <  2:1375016/150‑1 (MQ=255)
                                                                                                                                                aTTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGACACTGGCGTAGGTCGCCACAACGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAATATCGGCACTGATc     >  2:1322340/1‑150 (MQ=255)
                                                                                                                                                 tttACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAATATCGGCACTGATCg    <  1:1558379/150‑1 (MQ=255)
                                                                                                                                                   tACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAATATCGGCACTGATCGTc  >  1:477268/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCATTGTCCAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAATATCGGCACTGATCGTC  >  NC_000913/1575074‑1575370

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: