Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 715,347 C→A 6.0% intergenic (+149/‑65) pgm → / → ybfP phosphoglucomutase/lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913715,3470CA 6.0% 314.6 / 8.2 100intergenic (+149/‑65)pgm/ybfPphosphoglucomutase/lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (46/48);  new base A (2/4);  total (48/52)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

Only 100 of 105 total aligned reads displayed.

AACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTG  >  NC_000913/715198‑715490
                                                                                                                                                     |                                                                                                                                               
aaCACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                 <  1:1009045/150‑1 (MQ=255)
aaCACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                 <  2:1043765/150‑1 (MQ=255)
aaCACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                 <  2:235527/150‑1 (MQ=255)
aaCACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                 <  1:390214/150‑1 (MQ=255)
aaCACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                 >  1:1096971/1‑150 (MQ=255)
 acacATTTAATAAAAAAATGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCa                                                                                                                                                <  1:156708/150‑1 (MQ=255)
 acacATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCa                                                                                                                                                <  2:1085923/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)
             aaaaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCGGATGGTGCAATATAGTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGAATCCgg                                                                                                                                    >  1:1711759/1‑150 (MQ=255)
             aaaaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                    <  1:750358/150‑1 (MQ=25)
             aaaaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                    <  2:236217/150‑1 (MQ=255)
               aaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                  <  2:1386473/150‑1 (MQ=255)
               aaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                  <  1:723959/150‑1 (MQ=255)
                aaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                                                                 >  1:799902/1‑150 (MQ=255)
                 aaGGGCGGTCGCAAGATCGCACTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTAAAGCGATTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGccta                                                                                                                                <  1:1561386/150‑4 (MQ=21)
                  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                                                                                                                         <  2:726376/150‑1 (MQ=255)
                        gTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACg                                                                                                                         >  1:438835/1‑150 (MQ=255)
                         tCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGACTACGAGGATGGTGCAATATATTGAATTTTAGCGATTTTGTCGGCCGGATAAGGCGTTCACGCCGCCTCCGGCAAAAACTACGa                                                                                                                        >  2:321165/1‑150 (MQ=25)
                           gCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAac                                                                                                                      <  2:1152586/150‑1 (MQ=255)
                           gCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAac                                                                                                                      <  1:849539/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                                                                                                  >  1:1822991/1‑150 (MQ=255)
                                               gTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGt                                                                                                  <  2:1096971/150‑1 (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)
                                                                         gcgcTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGAcc                                                                        <  1:47557/150‑1 (MQ=255)
                                                                           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)
                                                                                                                      ttAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTa                           >  1:908948/1‑150 (MQ=255)
                                                                                                                        aaGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTAct                         >  2:1790629/1‑150 (MQ=255)
                                                                                                                          gCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTActgc                       <  1:1179969/150‑1 (MQ=255)
                                                                                                                          gCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTActgc                       >  1:15/1‑150 (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:124564/150‑1 (MQ=255)
                                                                                                                                    aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATg             >  2:556429/1‑150 (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   >  2:1028063/1‑150 (MQ=255)
                                                                                                                                               ggCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGGCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTg  >  2:1569614/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
AACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTG  >  NC_000913/715198‑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: