Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,223,932 G→A 6.6% intergenic (+25/+347) ycgI → / ← minE pseudogene/cell division topological specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,223,9320GA 6.6% 404.1 / 11.9 136intergenic (+25/+347)ycgI/minEpseudogene/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base G (65/62);  new base A (6/3);  total (71/65)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

Only 100 of 137 total aligned reads displayed.

CAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAA  >  NC_000913/1223785‑1224081
                                                                                                                                                   |                                                                                                                                                     
gaGGCAGCTATTAAAGCCTGATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGgc                                                                                                                                                     <  2:1203593/149‑1 (MQ=255)
ccggCCGCTATTAAAGCCTCATGCAGCTGTACGTTTAGCGGACATCTCGGGGTGGTGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTGCGgc                                                                                                                                                     <  2:1404072/148‑1 (MQ=17)
cAGGCAGCTATTAAAGCCTCATTCAGCATTACGTTTAGCGGACATCTCGGGGGTGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGgc                                                                                                                                                     <  2:1755269/150‑1 (MQ=255)
cAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGgc                                                                                                                                                     >  1:978224/1‑150 (MQ=255)
 aGGCAGCTATTAAAGCCTCAATCAGCAGTACGTGTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGgcg                                                                                                                                                    <  2:1531886/150‑1 (MQ=255)
    cAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAg                                                                                                                                                 <  2:639459/150‑1 (MQ=255)
    cAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAg                                                                                                                                                 >  1:1848124/1‑150 (MQ=255)
    cAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAg                                                                                                                                                 >  2:30553/1‑150 (MQ=255)
    cAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAg                                                                                                                                                 >  1:1690036/1‑150 (MQ=255)
      gCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGtt                                                                                                                                               >  2:764629/1‑150 (MQ=255)
       cTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCCGAAGCGCAGttt                                                                                                                                              <  1:634146/150‑1 (MQ=255)
       cTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGGGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCCGAAGCGCAGttt                                                                                                                                              <  1:1395384/150‑1 (MQ=17)
        tATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGtttt                                                                                                                                             >  1:610065/1‑150 (MQ=255)
        tATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGtttt                                                                                                                                             >  2:778024/1‑150 (MQ=255)
         aTTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTc                                                                                                                                            >  2:1714037/1‑150 (MQ=255)
         aTTAAAGCCTCATTCAGCAGTACGTTAAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTc                                                                                                                                            <  1:225468/150‑1 (MQ=255)
           tAAAGCCTCATTCCGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCCGAAGCGCAGtttaatc                                                                                                                                          <  1:388686/150‑5 (MQ=17)
           tAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGt                                                                                                                                          >  1:830713/1‑150 (MQ=255)
           tAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGt                                                                                                                                          <  1:711068/150‑1 (MQ=255)
            aaaGCCTCATTCACCAGTACGTTTAGCGGACATCGCGGGGTGGGGTCTTGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAATCTGCGCTTCGTCGCAGTTTTCGtt                                                                                                                                         <  2:1477821/150‑1 (MQ=18)
             aaGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGttt                                                                                                                                        <  2:1451515/150‑1 (MQ=255)
              aGCCTCTTTCCGCAGTTCGTTTATCGGAAATCTCGGGGTTGGGTATAGCCATGGTGGCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGGGTGATCTGGTCGGTCTGACCTTCATCTATTAAACTGCGCTTCGGCGCAGTTTTCGTTTa                                                                                                                                       <  1:893927/150‑1 (MQ=18)
              aGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTa                                                                                                                                       <  1:1278375/150‑1 (MQ=255)
               gCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTAc                                                                                                                                      >  1:552864/1‑150 (MQ=255)
                ccTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACa                                                                                                                                     <  2:829393/150‑1 (MQ=255)
                  tCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAgg                                                                                                                                   >  2:877843/1‑150 (MQ=255)
                   cATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGa                                                                                                                                  >  2:1484567/1‑150 (MQ=255)
                      tCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGt                                                                                                                               <  1:877843/150‑1 (MQ=255)
                      tCAGCAGTACGTGTAGCGGACAGCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCTTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGt                                                                                                                               <  1:148445/150‑1 (MQ=255)
                       cagcagTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGtt                                                                                                                              >  2:1806046/1‑150 (MQ=255)
                           agTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGaaa                                                                                                                          >  2:1874466/1‑150 (MQ=255)
                                gTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGaa                                                                                                                     >  1:1635948/1‑150 (MQ=255)
                                 tttAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGACCGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGaaa                                                                                                                    >  1:1057928/1‑150 (MQ=255)
                                 tttAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGaaa                                                                                                                    >  2:807001/1‑150 (MQ=255)
                                 tttAGCGGACATCTCGGGGTGGGGTAATGCAATGGTGCCTGTGTGTAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCAGAGCATCAACGATTAAACTGCGCTGCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGaaa                                                                                                                    <  1:100101/150‑1 (MQ=25)
                                        gACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTgg                                                                                                             <  1:777913/150‑1 (MQ=255)
                                         aCATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTggg                                                                                                            <  1:1732065/150‑1 (MQ=255)
                                           aTCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGc                                                                                                          <  2:804950/150‑1 (MQ=255)
                                           aTCTCGGGGGGGGGTCTAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGc                                                                                                          <  2:353986/150‑1 (MQ=255)
                                            tctcGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCa                                                                                                         <  1:142194/150‑1 (MQ=255)
                                             ctcGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAATGGAAAATTCTGGGGCaa                                                                                                        >  1:158373/1‑150 (MQ=255)
                                             ctcGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCaa                                                                                                        >  1:36075/1‑150 (MQ=255)
                                              tcGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCaaa                                                                                                       >  1:1547967/1‑150 (MQ=255)
                                                   gtggggtATAGCCTTGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTTGCTGGTGTGTACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAg                                                                                                  <  2:627933/150‑1 (MQ=255)
                                                      gtgtATAGCCATGGTGCCGGTGTGGAATCCCCGTGTAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGccc                                                                                               <  2:795886/148‑1 (MQ=255)
                                                          ataGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCa                                                                                           >  1:9806/1‑150 (MQ=255)
                                                          ataGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCa                                                                                           >  1:2767/1‑150 (MQ=255)
                                                            aGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGt                                                                                         >  2:1776591/1‑150 (MQ=255)
                                                                   gTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCAAGTTAcggcg                                                                                  >  1:1750205/1‑150 (MQ=255)
                                                                     gCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCggg                                                                                >  1:1060195/1‑150 (MQ=255)
                                                                       cGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGaa                                                                              >  1:1612360/1‑150 (MQ=255)
                                                                        ggTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGTTGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGaaa                                                                             >  2:1451030/1‑150 (MQ=255)
                                                                         gtgtGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAc                                                                            <  1:1445279/150‑1 (MQ=255)
                                                                          ggtgGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGAGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAcc                                                                           <  1:1825505/149‑1 (MQ=255)
                                                                           gtgGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCt                                                                          >  2:614192/1‑150 (MQ=255)
                                                                            gggAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTTGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTc                                                                         <  1:1308702/149‑1 (MQ=255)
                                                                             ggAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAACCCtaa                                                                        >  1:1345720/1‑148 (MQ=255)
                                                                                  ccccGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTa                                                                   >  1:923463/1‑150 (MQ=255)
                                                                                   cccGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCCGTTATCGTTTACAGGATGTTGAAAGGGATAATTCTGGGGCATAAAAAGCCAGCCAGTTACGGCGGGAAACCTCATCCTAt                                                                  >  2:1147677/1‑150 (MQ=38)
                                                                                    ccGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATg                                                                 <  1:387785/150‑1 (MQ=255)
                                                                                    ccGGGGAACGCGGGGGCTGGGGTGAACTGGGCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATg                                                                 <  2:123515/150‑1 (MQ=255)
                                                                                      gTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTACACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGGTACGGCGGGAAACCTCATCCTATggg                                                               <  1:1368642/150‑1 (MQ=255)
                                                                                       tGGAACGCGGTGGCGGGTGTGAACTGGTGGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGAAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATggta                                                              <  2:839337/150‑3 (MQ=255)
                                                                                        ggAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGgag                                                             <  1:1001666/150‑1 (MQ=255)
                                                                                        ggAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGgag                                                             <  1:778024/150‑1 (MQ=255)
                                                                                         gAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGgaga                                                            >  2:8444/1‑150 (MQ=255)
                                                                                            cgcgGTGGCTGGTGTGAACTGGTCCTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACa                                                         <  2:566846/150‑1 (MQ=255)
                                                                                             gcgGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACaa                                                        >  1:182701/1‑150 (MQ=255)
                                                                                               ggTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATg                                                      <  1:1635025/150‑1 (MQ=255)
                                                                                                 tGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGaa                                                    <  2:543595/150‑1 (MQ=255)
                                                                                                     tGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGaataat                                                <  2:487408/150‑1 (MQ=255)
                                                                                                       gtgtGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGa                                              <  2:1084815/150‑1 (MQ=255)
                                                                                                        tgtgAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGaa                                             >  2:1668475/1‑150 (MQ=255)
                                                                                                           gAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAAtt                                          >  1:260599/1‑150 (MQ=255)
                                                                                                                ggTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCggg                                     >  2:1116781/1‑150 (MQ=255)
                                                                                                                ggTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGATTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTAATCCTATGGGAGAACAATGAATAATGAAATTGCggg                                     >  1:477231/1‑150 (MQ=255)
                                                                                                                 gTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCgggg                                    <  1:1556684/150‑1 (MQ=255)
                                                                                                                      tCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTatc                               <  1:1030054/150‑1 (MQ=255)
                                                                                                                             aTCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTccc                        <  2:1690036/150‑1 (MQ=255)
                                                                                                                             aTCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTccc                        <  2:1547967/150‑1 (MQ=255)
                                                                                                                             aTCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTccc                        >  1:483691/1‑150 (MQ=255)
                                                                                                                             aTCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTccc                        >  1:1546903/1‑150 (MQ=255)
                                                                                                                               cAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCTTCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAg                      >  2:1476142/1‑150 (MQ=255)
                                                                                                                               cAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAg                      >  2:969893/1‑150 (MQ=255)
                                                                                                                               cAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAg                      <  2:125738/150‑1 (MQ=255)
                                                                                                                               cAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGCATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAg                      >  2:1437565/1‑150 (MQ=255)
                                                                                                                                 aCGATTAAACTGCGCTTAGGCGCAGTTTTCGTTTACAGGAGGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGTCGGGCTACCTCATCCTATGGGAGAACAATGAATAAAGAACTGGCGGTGTTATCATAGCCCAGta                    >  2:947122/1‑150 (MQ=255)
                                                                                                                                  cGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGtat                   >  2:868722/1‑150 (MQ=255)
                                                                                                                                  aacgaaaaCTGCGCCGAAGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACTGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGtat                   >  2:881313/6‑150 (MQ=255)
                                                                                                                                   acgaaaaCTGCGCCGAAGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATTATCTCCCAGtata                  >  2:1598699/5‑150 (MQ=255)
                                                                                                                                   acgaaaaCTGCGCCGAAGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGtata                  >  1:723876/5‑150 (MQ=255)
                                                                                                                                    aTTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTatat                 <  1:1874466/150‑1 (MQ=255)
                                                                                                                                    aTTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTatat                 >  1:589718/1‑150 (MQ=255)
                                                                                                                                     ttAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATc                <  1:1806046/150‑1 (MQ=255)
                                                                                                                                       aaaCTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCa              <  2:1389177/150‑1 (MQ=255)
                                                                                                                                        aaCTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCAt             >  2:1150117/1‑150 (MQ=255)
                                                                                                                                          cTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATAc           >  1:163448/1‑150 (MQ=255)
                                                                                                                                              gcTTCGTCGAAGTTTTCGTATACAGGATGTTGATAGGGAAAAGTCTGGGTCATAAAAATCCCGCCAGTTACGGCGGGAAACCTCATCCTATCGGAGAACAAAGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAAc       <  1:109605/150‑1 (MQ=255)
                                                                                                                                               cTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACa      <  2:1645486/150‑1 (MQ=255)
                                                                                                                                               cTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACa      <  2:1060195/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAA  >  NC_000913/1223785‑1224081

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: