Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,747,025 G→C 9.0% intergenic (+20/‑106) ydhR → / → ydhS putative monooxygenase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,747,0250GC 9.0% 529.4 / 32.7 155intergenic (+20/‑106)ydhR/ydhSputative monooxygenase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (72/69);  new base C (5/9);  total (77/78)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 157 total aligned reads displayed.

AAAAGCGCCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATC  >  NC_000913/1746877‑1747172
                                                                                                                                                    |                                                                                                                                                   
aaaaGCGCCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCgg                                                                                                                                                     >  2:764635/1‑149 (MQ=255)
   aGCGCCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTa                                                                                                                                                  <  4:1374535/149‑1 (MQ=255)
        ccTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGcc                                                                                                                                             <  1:1316862/149‑1 (MQ=255)
            gCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCttt                                                                                                                                         >  2:1432037/1‑149 (MQ=255)
            gCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGCTAACCGCCCttg                                                                                                                                         <  3:191439/149‑2 (MQ=255)
             ccTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGCTAACCGCCCttga                                                                                                                                        <  2:1251994/149‑3 (MQ=255)
                aTCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCAt                                                                                                                                     >  2:1053331/1‑149 (MQ=255)
                 tCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATc                                                                                                                                    >  1:1203888/1‑149 (MQ=255)
                  cTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGCTAACCGCCCttgattaaa                                                                                                                                   <  2:27785/149‑8 (MQ=255)
                    ggAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCt                                                                                                                                 >  1:976305/1‑149 (MQ=255)
                    ggAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCt                                                                                                                                 >  3:455052/1‑149 (MQ=255)
                          aCATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGa                                                                                                                           >  4:1131263/1‑149 (MQ=255)
                           cATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGaa                                                                                                                          >  3:1118360/1‑149 (MQ=255)
                           cATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGaa                                                                                                                          >  4:631545/1‑149 (MQ=255)
                                        tGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGt                                                                                                                               >  1:252258/1‑131 (MQ=255)
                                        tGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGt                                                                                                                               <  2:252258/131‑1 (MQ=255)
                                            aaaTCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTaaa                                                                                                         <  3:817273/149‑1 (MQ=255)
                                                ctcGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAAtc                                                                                                     <  3:1101568/149‑1 (MQ=255)
                                                 tcGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAAtct                                                                                                    >  1:186486/1‑149 (MQ=255)
                                                  cGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAAtctc                                                                                                   <  2:611260/149‑1 (MQ=255)
                                                    gCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCaa                                                                                                 <  3:989528/149‑1 (MQ=255)
                                                    gCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATATCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTGAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCaa                                                                                                 >  4:1003717/1‑149 (MQ=255)
                                                            gAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTa                                                                                                           <  4:88445/131‑1 (MQ=255)
                                                               gTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCACCCAAATTCATc                                                                                      <  4:1110689/149‑1 (MQ=255)
                                                                ttGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCt                                                                                     <  1:378979/149‑1 (MQ=255)
                                                                     gCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATaa                                                                                <  2:857872/149‑1 (MQ=255)
                                                                     gCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGCTAACCGCCCttga                                                                                                                                        >  1:123683/1‑91 (MQ=255)
                                                                     gCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGCTAACCGCCCttg                                                                                                                                         >  1:206986/1‑91 (MQ=255)
                                                                     gCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGCTAACCGCCCttg                                                                                                                                         <  2:206986/92‑2 (MQ=255)
                                                                      ccAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAAttt                                                                                                                      <  3:124759/110‑1 (MQ=255)
                                                                      ccAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAAttt                                                                                                                      >  4:124759/1‑110 (MQ=255)
                                                                      ccAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAac                                                                               >  1:472127/1‑149 (MQ=255)
                                                                      ccAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAac                                                                               >  2:510322/1‑149 (MQ=255)
                                                                      ccAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCGCAACCAAACTTATCTGATAac                                                                               >  4:1137343/1‑149 (MQ=255)
                                                                      ccAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGCTAACCGCCCttgat                                                                                                                                       >  3:913086/1‑90 (MQ=255)
                                                                      ccAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGCTAACCGCCCttgat                                                                                                                                       <  4:913086/93‑4 (MQ=255)
                                                                       cAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGTGGTTAgc                                                                                                                                                <  2:1269217/83‑1 (MQ=255)
                                                                       cAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAaca                                                                              >  3:327607/1‑149 (MQ=255)
                                                                       cAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAaca                                                                              >  1:82320/1‑149 (MQ=255)
                                                                       cAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAaca                                                                              >  3:298447/1‑149 (MQ=255)
                                                                           gTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCACATTCTAAAAATCTCAACCAAACTTATCTGATAACACTaa                                                                          <  2:1360781/149‑1 (MQ=255)
                                                                                 gATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGa                                                                    >  2:467642/1‑149 (MQ=255)
                                                                                   tGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGaaa                                                                  <  4:789207/149‑1 (MQ=255)
                                                                                   tGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGaaa                                                                  <  4:422971/149‑1 (MQ=255)
                                                                                   tGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGaaa                                                                  <  2:1022520/149‑1 (MQ=255)
                                                                                      cAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGaa                                                               <  1:1500736/149‑1 (MQ=255)
                                                                                      cAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGaa                                                               <  2:1229044/149‑1 (MQ=255)
                                                                                       aaTGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCaaa                                                                                                                <  1:1122302/99‑1 (MQ=255)
                                                                                       aaTGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCaaa                                                                                                                >  2:1122302/1‑99 (MQ=255)
                                                                                       aaTGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAAt                                                              >  1:229257/1‑149 (MQ=255)
                                                                                       aaTGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAAt                                                              >  3:1005951/1‑149 (MQ=255)
                                                                                         tGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGc                                                            <  1:467642/149‑1 (MQ=255)
                                                                                          gAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCaa                                                                                                 <  1:1398829/111‑1 (MQ=255)
                                                                                            gCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGccc                                                                                                                                            <  1:89872/66‑1 (MQ=255)
                                                                                            gCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGccc                                                                                                                                            >  2:89872/1‑66 (MQ=255)
                                                                                            gCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAGATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTa                                                         <  3:885148/149‑1 (MQ=255)
                                                                                            gCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTa                                                         <  3:1137343/149‑1 (MQ=255)
                                                                                               aCTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAg                                                      >  2:97704/1‑149 (MQ=255)
                                                                                               aCTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAg                                                      <  2:357139/149‑1 (MQ=255)
                                                                                                cTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGc                                                            >  4:1025392/1‑142 (MQ=255)
                                                                                                 ttAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgt                                                    <  1:10653/149‑1 (MQ=255)
                                                                                                 ttAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgt                                                    <  2:186486/149‑1 (MQ=255)
                                                                                                 ttAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgt                                                    <  1:97704/149‑1 (MQ=255)
                                                                                                    gTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgtaag                                                 >  2:1311266/1‑149 (MQ=255)
                                                                                                     tCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCaaa                                                                                            >  2:66588/1‑105 (MQ=255)
                                                                                                     tCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgtaagt                                                <  4:298447/149‑1 (MQ=255)
                                                                                                     tCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgtaagt                                                >  3:792987/1‑149 (MQ=255)
                                                                                                     tCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGCGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgtaagt                                                >  1:518387/1‑149 (MQ=255)
                                                                                                       aaaTCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgtaagtaa                                              >  3:220220/1‑149 (MQ=255)
                                                                                                        aatcaatcAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAAc                                             >  4:591005/1‑149 (MQ=255)
                                                                                                        aatcaatcAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAAc                                             <  1:1432037/149‑1 (MQ=255)
                                                                                                          tcaatcaGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACaa                                           <  2:1321392/149‑1 (MQ=255)
                                                                                                           caatcaGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAAt                                          <  3:379462/149‑1 (MQ=255)
                                                                                                              tcaGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGaa                                       >  2:763739/1‑149 (MQ=255)
                                                                                                                  gCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGaaaaaa                                   >  3:1358481/1‑149 (MQ=255)
                                                                                                                     aaaCTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAAtt                                <  3:631545/149‑1 (MQ=255)
                                                                                                                     aaaCTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTATGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAAtt                                <  1:764635/149‑1 (MQ=255)
                                                                                                                      aaCTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTg                               >  4:1563093/1‑149 (MQ=255)
                                                                                                                         tCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGaaa                                                                  >  2:1380480/1‑111 (MQ=255)
                                                                                                                         tCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCg                                                           >  4:1567296/1‑118 (MQ=255)
                                                                                                                         tCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTa                            >  3:1436280/1‑149 (MQ=255)
                                                                                                                         tCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTa                            >  2:269521/1‑149 (MQ=255)
                                                                                                                           gCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTAtt                          <  1:510322/149‑1 (MQ=255)
                                                                                                                             cTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATGGCTAtttt                        >  3:617172/1‑147 (MQ=255)
                                                                                                                               gACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTgg                      <  3:1140563/149‑1 (MQ=255)
                                                                                                                                   gAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCt                                                                                     >  3:151935/1‑82 (MQ=255)
                                                                                                                                     aTTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCg                >  1:326224/1‑149 (MQ=255)
                                                                                                                                      tttAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCgg               >  4:319677/1‑149 (MQ=255)
                                                                                                                                      tttAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCgg               <  4:1118360/149‑1 (MQ=255)
                                                                                                                                         aaTCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGcc            >  1:1328587/1‑149 (MQ=255)
                                                                                                                                            aaaGGGCGCTAACCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAAtt                                <  2:1165302/125‑1 (MQ=255)
                                                                                                                                              aGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACgg       >  2:892618/1‑149 (MQ=255)
                                                                                                                                               gggCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACggg      <  4:521505/149‑1 (MQ=255)
                                                                                                                                                ggCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACgggg     <  3:10673/149‑1 (MQ=255)
                                                                                                                                                ggCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACg        >  4:571962/1‑146 (MQ=255)
                                                                                                                                                ggCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACg        <  3:571962/146‑1 (MQ=255)
                                                                                                                                                ggCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTGCTCAAATTCTAAAAATCTCAACCAAACTTAACTGATAACACTAAATTCGAAAGACTGGGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCTTACgggg     >  3:1245354/1‑149 (MQ=255)
                                                                                                                                                 gCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGa    <  1:32723/149‑1 (MQ=255)
                                                                                                                                                  cGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGAt   <  4:1436280/149‑1 (MQ=255)
                                                                                                                                                  cGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGAt   <  2:229257/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AAAAGCGCCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATC  >  NC_000913/1746877‑1747172

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: