Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913266,3580TG43.5% 97.9 / 20.6 71T66P (ACC→CCC) yafYlipoprotein, inner membrane; degP regulator; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base T (9/30);  new base G (30/0);  total (41/30)
Fisher's exact test for biased strand distribution p-value = 7.88e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACA  >  NC_000913/266213‑266507
                                                                                                                                                    |                                                                                                                                                      
tcctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGGAAAAACTGCCAgggggg                                                                                                                                                     >  1:1139536/1‑149 (MQ=255)
 cctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACCGCCAGGGGGGGcc                                                                                                                                                  >  1:223312/1‑151 (MQ=255)
   tGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAACCTTCAAAGGGGGGCCgg                                                                                                                                                >  1:693665/1‑151 (MQ=255)
    ggTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCgg                                                                                                                                                <  1:231406/150‑1 (MQ=255)
    ggTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCgg                                                                                                                                                <  1:181296/150‑1 (MQ=255)
     gTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGAt                                                                                                                                              <  1:123339/151‑1 (MQ=255)
     gTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTTCAAGGGGGGGCCgg                                                                                                                                                >  1:448007/1‑149 (MQ=255)
       aTTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATa                                                                                                                                             <  2:531716/150‑1 (MQ=255)
       aTTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGAGACCAGATa                                                                                                                                             <  2:608130/150‑1 (MQ=255)
         tCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAag                                                                                                                                           <  1:1146695/150‑1 (MQ=255)
          ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAagt                                                                                                                                          <  1:775188/150‑1 (MQ=255)
          ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGGAAAAACCGCAAGGGGGGACCgggaaa                                                                                                                                            >  1:257644/1‑144 (MQ=255)
           cAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGGCGTAAAACCCGCAAGGTGGGACCggggaagg                                                                                                                                          >  2:544935/1‑143 (MQ=255)
           cAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGAACGGAAAACCCGCAAGGGGGGGCCGGa                                                                                                                                               >  2:665378/1‑144 (MQ=255)
            aGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTTAAACCTGCCAGGGGGGGCCGGag                                                                                                                                              >  1:686505/1‑143 (MQ=255)
             gTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCCAGGGGGGGCCCGGTAAAGAGTc                                                                                                                                      >  1:93247/1‑151 (MQ=255)
              tttGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGACCGGATAagta                                                                                                                                         >  2:927557/1‑147 (MQ=255)
                 gATGTTATCGGCAATACACAGCGCAGTCATCAGTGGCTGGACGGTTTTTACCAGCGGATCCACGAGGTAGTAGTAACGGCCATCCTCGCGGCCATCGAGCGACTGGCGGGTACGTAAAACCAGCAAGGTGGGACCGGATAAGTAGTCa                                                                                                                                     <  2:375139/148‑1 (MQ=255)
                   tGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAAGGGGGGCCGGGTAAGTAGTaaa                                                                                                                                    >  1:993576/1‑147 (MQ=255)
                   tGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAAGGAAAAACTGCAAAGGGGGGCCGGGTaaaatagcaaac                                                                                                                                  >  2:832106/1‑139 (MQ=255)
                         cGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGGCCGGGTAAATAGGCAACCTGTGAcc                                                                                                                           >  1:1194287/1‑150 (MQ=255)
                                 cacaGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGGAAAAACCACAAAGGGGGGCCGGGGAAGGAGGTAAACTGTGGCCActcct                                                                                                                     >  2:871/1‑146 (MQ=255)
                                      cgcATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGAca                                                                                                             <  2:946322/151‑1 (MQ=255)
                                      cgcATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTCGCGGGTACGTAAATCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGAca                                                                                                             <  2:152389/151‑1 (MQ=255)
                                         aTTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGAcatca                                                                                                          <  2:1150175/151‑1 (MQ=255)
                                         aTTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGGGGGGACGGAAAAGTAGTCAACCTGTGACCATTCTTCCCTGagcac                                                                                                           >  1:1026960/1‑146 (MQ=255)
                                             atcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAACCTGCAAGGGGGGACCGGGGAAAGAGTCAACCTGGGACCCTTATTCCCCGAAATCATcc                                                                                                       >  2:546758/1‑150 (MQ=255)
                                              tcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAACTAAAAAACGCAAGGAGGGGCCCGGTAAATAGTCAACCTGTGGACACTCCTCGCtcgaataatac                                                                                                       >  2:1343184/1‑139 (MQ=255)
                                               caGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAAACTGCCAAGGGGGGCCGGATAAAAAGTCAACCCGTGACCCTTCTTTTCTGAaatcatacctg                                                                                                    >  1:641744/1‑146 (MQ=255)
                                                aGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAACTAAAAACCTCAAGGGGGGGCCGGGGAAAAAGGCAAACTGTGGACATTTTTCGCTTACAACATCCtcg                                                                                                    >  1:1185168/1‑148 (MQ=255)
                                                    gCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCCAGGGGGGGCCGGATAAATAGTCAACCTGTGACCATTTTTCCCTGACATCATCCtggggggt                                                                                               >  1:317068/1‑149 (MQ=255)
                                                      tGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCTCAAGGGGGGAACCGCAAAGTAGTCAAACTGTGGACATTCTTCGCCGACCTCCACCttgtgg                                                                                                 >  2:887530/1‑147 (MQ=255)
                                                         aCGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAAACCTCCAAGGGGGGCCGGATAAGAAGGCCACCTGTGACCATTCTTCTCTGACATAATCCTGGGGGGGGAcc                                                                                           >  2:1157438/1‑150 (MQ=255)
                                                          cGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACAGGCGGGTACGTAAAAACTGCAAGGTGGGACCGGATAAGGAGTCAACCTGTGACCATTCTTCGCTGACAACATTCTGGTGGCTGACca                                                                                          >  2:1347609/1‑150 (MQ=255)
                                                           ggTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGACCGGGAAAATAGACAACCTGTGACCATTCTTTGCTGTCATCCTCCTGGGGGGTGAcc                                                                                           >  1:1309814/1‑148 (MQ=255)
                                                             ttttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACTTAAAACCTGCAAAGGGGGGCCCGGGAAGTAGTTAAACTGTGACCAATTTTTGCCGACATCATCCTGGGGGCTGACcaca                                                                                        >  1:263235/1‑149 (MQ=255)
                                                                tttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTGAAAACTGCAAGGTGGGGCCGGATAAGGAGTCAACCGGTGAACATTCTTCGCTGACATCATCCTGGGGGGTGACca                                                                                          >  2:706355/1‑144 (MQ=18)
                                                                     tGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAAGGAAAAACCGCAAGGGGGGACCGGGGAAATAGTCAACCTGTGACCCTTCTTCGCTGACATCATCCTGGTGGGTGACCACCTCACTTaaa                                                                              >  1:703067/1‑150 (MQ=21)
                                                                      gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGaaa                                                                              <  2:883343/149‑1 (MQ=255)
                                                                      gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAAcg                                                                            <  1:1034829/151‑1 (MQ=255)
                                                                      gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGGCCGGATAAATAGTCAACCTGTGACCATTCTTCGCTTACATCATCCTGGTGGGTGACCCCATCACTGaaa                                                                              >  1:1196353/1‑149 (MQ=255)
                                                                       cATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAAcgc                                                                           <  1:202434/151‑1 (MQ=255)
                                                                          atCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCCAGGGGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCggg                                                                        >  1:176322/1‑150 (MQ=255)
                                                                          atCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGCAAAACCTGCAAGGGGGGGGCGGGGAAGGAGACCACCTGTGACCATTCTTCGCTGACATTATCCTCGTGGGTGACCCCAACACTGAAACGCgg                                                                         >  1:949397/1‑150 (MQ=18)
                                                                           tCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAACCCCCAAAGTGGGGCCGGAGAAGTAGTCAACCTGTGACCATTCTTTCCTGACATCATCCTTGTGGCTGACCACATCCCTGAAAAGCgg                                                                         >  2:389471/1‑149 (MQ=255)
                                                                            ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTgg                                                                      <  1:724168/151‑1 (MQ=255)
                                                                             cACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAACCTGCAAGGGGGGGCCGGGTAAGTAGTCAAACTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGGCCCCATCACTGAAACGCgggggg                                                                     >  1:1259326/1‑151 (MQ=21)
                                                                                gAGGTAGTAGTAACTGCCATCCGCGCAGCCATCGGGCGACTGGCGGCTTAGGCAAAACCCCAACGGGGGGCCGGGATAGTAGGCAAACTGTGACCAATCCTCCCCGACAACATCCTGGTGGTGGTCCCCATCCCTTAAAAACTGTgggggg                                                                  >  2:365829/1‑149 (MQ=255)
                                                                                        agtaACTGCCATCCTCGCAGCCATCGGGCGACTGACGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCtt                                                          <  2:1202081/151‑1 (MQ=255)
                                                                                           aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTaaa                                                       <  1:1183613/151‑1 (MQ=255)
                                                                                           aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTaa                                                        <  2:1203849/150‑1 (MQ=255)
                                                                                           aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTaa                                                        <  1:992729/150‑1 (MQ=255)
                                                                                           aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGGAAAAACTGCAAGGGGGGGCCGGATAAAGAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGGGGCTGACCACATCACTGAAACGCGGGGGGGGGGGGGCCttta                                                        >  1:682008/1‑148 (MQ=18)
                                                                                                ccATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCCCAACAGGGGGGGCCGCTAAAGAGTCAAACCGTGGCCCCTTTTTCCTGACCTCATCCTGGGGGCTGACCACACACACGTAACGCGGGGGGGGGGGGGGCTTAAAatt                                                    >  2:894691/1‑149 (MQ=255)
                                                                                                   tCCTCGCAGCCATCGGGCGACTGGCGGGTACGCAAAACCTGCAAGGTGGGGCCGGCTAAGGAGTCAACCTGTGGCCATTCCTCCCTGACATCATCCCGGGGGCGGACCACATCACTGGAACGCGGGGGGGGGAGGGCCCTTAATTTTCt                                                 >  2:237363/1‑149 (MQ=255)
                                                                                                     cTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAGGGGGGGCCGGGTAAGTAGTCAACCCGTGACCATTCTTTGCTGACATCATCCCTGGGGCTGACCACATCACTGAAAAACGGGGGGGGGGGGGCCTTAAATTTACTGAt                                              >  1:478886/1‑150 (MQ=255)
                                                                                                       cGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATg                                             <  1:412626/149‑1 (MQ=255)
                                                                                                       cGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATg                                             <  2:814465/149‑1 (MQ=255)
                                                                                                         cAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCt                                         <  1:392724/151‑1 (MQ=255)
                                                                                                          aGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGACCGGGTAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGGTGACCACATCACTGAAACGCGGGGGGGGGGGGGCCTTAAATTTTCTGATGGCCtt                                        >  1:1223653/1‑151 (MQ=18)
                                                                                                            ccATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCGCATCTCTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCa                                      <  2:1039310/151‑1 (MQ=255)
                                                                                                               tCGGGCGACTGGCGGGTGCGTAAAAACTGCAAGGTGGGACCGGATAAGAAGTCAACCTGTGACCCTTCTTCGCTGACATCATCCTGGGGGCTGGACACATCACTGAAAAGCGGGGGGGGGAGGTGCTTTAATTTACTGATGGCCTTCAgg                                    >  2:1274933/1‑150 (MQ=255)
                                                                                                                cGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGt                                   <  1:1166831/150‑1 (MQ=255)
                                                                                                                 gggCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtca                                 <  1:480041/151‑1 (MQ=255)
                                                                                                                                 aCGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGACCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCg                  >  1:522552/1‑151 (MQ=255)
                                                                                                                                   gTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCa                <  1:1116019/151‑1 (MQ=255)
                                                                                                                                    tAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCa                <  1:383094/150‑1 (MQ=255)
                                                                                                                                    tAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCAt               >  1:381638/1‑151 (MQ=255)
                                                                                                                                      aaaaCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATg              <  1:1203656/151‑1 (MQ=255)
                                                                                                                                             gCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATg         <  1:921877/149‑1 (MQ=255)
                                                                                                                                                 aGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGaa     >  1:587068/1‑150 (MQ=255)
                                                                                                                                                   gTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACa  <  2:785758/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACA  >  NC_000913/266213‑266507

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