Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 715,346 T→G 5.7% intergenic (+148/‑66) pgm → / → ybfP phosphoglucomutase/lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913715,3460TG 5.7% 324.0 / 8.4 105intergenic (+148/‑66)pgm/ybfPphosphoglucomutase/lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (46/52);  new base G (2/4);  total (49/56)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

Only 100 of 106 total aligned reads displayed.

AAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTG  >  NC_000913/715197‑715490
                                                                                                                                                     |                                                                                                                                                
aaaCACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGtt                                                                                                                                                  >  2:1300456/1‑150 (MQ=255)
 aaCACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                 <  1:1009045/150‑1 (MQ=255)
 aaCACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                 >  1:1096971/1‑150 (MQ=255)
 aaCACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                 <  1:390214/150‑1 (MQ=255)
 aaCACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                 <  2:1043765/150‑1 (MQ=255)
  acacATTTAATAAAAAAATGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCa                                                                                                                                                <  1:156708/150‑1 (MQ=255)
  acacATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCa                                                                                                                                                <  2:1085923/150‑1 (MQ=255)
       tttAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                           <  1:1592999/150‑1 (MQ=255)
         taataaAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACCAGGATGGTGCTATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCa                                                                                                                                         >  2:1387828/1‑150 (MQ=255)
           ataaAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGAAAAACACAGAATTGCCTTATGCGCTACGCATATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAATGCGTTCCCGCCGCATc                                                                                                                                       >  2:368742/1‑150 (MQ=255)
            taaAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                      <  1:398556/150‑1 (MQ=255)
             aaaaaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                     >  2:1721136/1‑150 (MQ=255)
              aaaaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCGGATGGTGCAATATAGTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGAATCCgg                                                                                                                                    >  1:1711759/1‑150 (MQ=255)
              aaaaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                    <  1:750358/150‑1 (MQ=25)
              aaaaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                    <  2:236217/150‑1 (MQ=255)
               aaagaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATTACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTTCAATATATTGTATTTAAGCGATTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                   <  2:1153949/150‑1 (MQ=25)
                aaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                  <  1:723959/150‑1 (MQ=255)
                 aaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                                                                 >  1:799902/1‑150 (MQ=255)
                   cgggCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaaaa                                                                                                                               <  2:1490740/149‑1 (MQ=255)
                   aGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctac                                                                                                                               <  2:156968/150‑5 (MQ=25)
                   aGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaaaa                                                                                                                               <  1:1357798/150‑1 (MQ=255)
                      gCGGTCGCAAGATCGCCCTTTATTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGCGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAaaca                                                                                                                            <  1:920164/150‑1 (MQ=255)
                        ggTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAaacaac                                                                                                                          <  2:1405165/150‑1 (MQ=255)
                         gTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACg                                                                                                                         >  1:438835/1‑150 (MQ=255)
                         gTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACg                                                                                                                         <  2:726376/150‑1 (MQ=255)
                          tCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGACTACGAGGATGGTGCAATATATTGAATTTTAGCGATTTTGTCGGCCGGATAAGGCGTTCACGCCGCCTCCGGCAAAAACTACGa                                                                                                                        >  2:321165/1‑150 (MQ=25)
                            gCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAac                                                                                                                      <  1:849539/150‑1 (MQ=255)
                            gCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAac                                                                                                                      <  2:1152586/150‑1 (MQ=255)
                             cAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAaca                                                                                                                     >  2:1110462/1‑150 (MQ=255)
                                 aTCGCCCTTTTTTACGTATGACAAACACATAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGCATGGAGCAATAAATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACttt                                                                                                                 <  1:45207/150‑1 (MQ=255)
                                    gCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTc                                                                                                              >  2:432300/1‑150 (MQ=255)
                                     cccTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTca                                                                                                             <  2:81550/150‑1 (MQ=255)
                                      ccTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTcaa                                                                                                            <  2:1363402/150‑1 (MQ=255)
                                        ttttttACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTcaaca                                                                                                          <  1:1104711/150‑1 (MQ=255)
                                          ttttACGTATGACAAACACAGAATTGCCTGATGCGCAACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACaaa                                                                                                        <  2:142376/150‑1 (MQ=255)
                                            ttACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACt                                                                                                      <  2:1666088/150‑1 (MQ=255)
                                             tACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTg                                                                                                     <  1:657165/150‑1 (MQ=255)
                                                gTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGt                                                                                                  <  2:1096971/150‑1 (MQ=255)
                                                gTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGt                                                                                                  >  1:1822991/1‑150 (MQ=255)
                                                   tGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGc                                                                                               <  1:123266/150‑1 (MQ=255)
                                                    gACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCt                                                                                              >  2:1179969/1‑150 (MQ=255)
                                                      cAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATAAATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCa                                                                                            >  1:1766864/1‑150 (MQ=255)
                                                            caGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGaaa                                                                                      >  1:764716/1‑150 (MQ=255)
                                                             aGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAAt                                                                                     >  1:852237/1‑150 (MQ=255)
                                                               aaTTGCCTTATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATcc                                                                                   <  2:158389/150‑1 (MQ=255)
                                                               aaTTGCCTGATGCGCTACGATTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACTAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATcc                                                                                   <  1:1141690/150‑1 (MQ=255)
                                                                aTTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATccc                                                                                  >  1:464881/1‑150 (MQ=255)
                                                                  ggCCGGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCaa                                                                                <  2:1575640/149‑1 (MQ=255)
                                                                   gCCTGATGCGCTACGCTTATCAGGCCTACGATGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAAt                                                                               <  1:1300456/150‑1 (MQ=255)
                                                                       gATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAg                                                                           <  2:927595/150‑1 (MQ=255)
                                                                        aTGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGa                                                                          <  2:1057821/150‑1 (MQ=255)
                                                                        aTGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTACACGCCGCATCCGGCAAAACCAACGAACACTTTGTCAACATACTGCGTAGCTCATGGAAATCCCAATGAAGa                                                                          >  2:1610246/1‑150 (MQ=25)
                                                                            gcTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCaa                                                                      >  1:618401/1‑150 (MQ=255)
                                                                              tACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACa                                                                    >  2:472154/1‑150 (MQ=255)
                                                                               aCGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAg                                                                   <  2:1212110/150‑1 (MQ=255)
                                                                                cGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAgg                                                                  >  1:157273/1‑150 (MQ=255)
                                                                                     aTCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATACGGCGTTCACGCCGCATCCGGCAAAACCAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCt                                                             >  2:288921/1‑150 (MQ=255)
                                                                                       cAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAATGCGTTCACGCCGCATCCTGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTgg                                                           >  1:152616/1‑150 (MQ=255)
                                                                                        aGGCCTACGAGGATGGTGCAATAAAATGAATTTAAGAGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACCAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGgt                                                          >  1:1458008/1‑150 (MQ=255)
                                                                                          gCCTACGAGGATGGTGCAATATATTTAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGgtcg                                                        <  1:1219008/150‑1 (MQ=255)
                                                                                            cTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGgtcgtc                                                      <  1:765267/150‑1 (MQ=255)
                                                                                              aCGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCAt                                                    >  2:1034822/1‑150 (MQ=255)
                                                                                               cGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATa                                                   >  2:1600817/1‑150 (MQ=255)
                                                                                                  ggATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAAAAAACTGAGTAGATCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAgtg                                                >  1:852126/1‑150 (MQ=255)
                                                                                                       gTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACt                                           >  1:1008300/1‑150 (MQ=255)
                                                                                                        tGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTa                                          <  2:799902/150‑1 (MQ=255)
                                                                                                           aaTATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAAtt                                       <  2:795628/150‑1 (MQ=255)
                                                                                                            atatatTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTa                                      <  1:621826/150‑1 (MQ=255)
                                                                                                             tatatTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTAc                                     >  2:1811152/1‑150 (MQ=255)
                                                                                                               tatTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCg                                   >  1:468278/1‑150 (MQ=255)
                                                                                                                atTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTAcctc                                  >  2:1214509/1‑148 (MQ=255)
                                                                                                                atTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACAAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGc                                  >  2:950983/1‑150 (MQ=255)
                                                                                                                atTGAATTCAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGc                                  >  1:538340/1‑150 (MQ=255)
                                                                                                                 ttGAATTTAATCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGcc                                 <  2:376295/150‑1 (MQ=255)
                                                                                                                  tGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCcac                               <  1:396847/150‑1 (MQ=255)
                                                                                                                   gAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAACCTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCcac                               >  1:1145034/1‑150 (MQ=255)
                                                                                                                     ttttAAGCGAGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGACATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACAt                             <  2:282007/149‑1 (MQ=255)
                                                                                                                       ttAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTa                           >  2:1415872/1‑150 (MQ=255)
                                                                                                                         aaGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTAct                         >  2:1790629/1‑150 (MQ=255)
                                                                                                                           gCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTActgc                       >  1:15/1‑150 (MQ=255)
                                                                                                                           gCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTActgc                       <  1:1179969/150‑1 (MQ=255)
                                                                                                                            cGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTActgct                      >  2:1173603/1‑150 (MQ=255)
                                                                                                                             gATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTActgctg                     <  1:1383767/150‑1 (MQ=255)
                                                                                                                                tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACt                  <  2:157273/150‑1 (MQ=255)
                                                                                                                                  tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCTAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGc                <  1:1146796/150‑1 (MQ=255)
                                                                                                                                  tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGc                <  1:226888/150‑1 (MQ=255)
                                                                                                                                  tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGc                >  2:755948/1‑150 (MQ=255)
                                                                                                                                    tttgCCGGATGCGGCGTGAACGCCTTATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCAt              >  1:1789414/4‑150 (MQ=255)
                                                                                                                                     aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATg             >  1:585376/1‑150 (MQ=255)
                                                                                                                                     aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATg             >  2:556429/1‑150 (MQ=255)
                                                                                                                                     aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATg             <  2:124564/150‑1 (MQ=255)
                                                                                                                                     aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAAAAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTTCTGACTGCATg             <  1:1629817/150‑1 (MQ=255)
                                                                                                                                      ggCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACGAATTACCGCCACATTACTGCGGACTGCATgc            >  2:189196/1‑150 (MQ=255)
                                                                                                                                       gCCGGATGCGGCGTGAACGCCTTATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATgcg           >  1:722336/1‑150 (MQ=255)
                                                                                                                                         cGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGcc         <  1:25788/150‑1 (MQ=255)
                                                                                                                                         cGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACAATACTGCTGACTGCATGCGcc         >  1:745733/1‑150 (MQ=255)
                                                                                                                                              aaGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAAcc    <  1:1106811/150‑1 (MQ=255)
                                                                                                                                               aGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCt   >  1:677236/1‑150 (MQ=255)
                                                                                                                                               aGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCt   >  2:1028063/1‑150 (MQ=255)
                                                                                                                                                ggCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGGCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTg  >  2:1569614/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
AAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTG  >  NC_000913/715197‑715490

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