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

TCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACA  >  NC_000913/266213‑266507
                                                                                                                                                    |                                                                                                                                                      
tcctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGG‑AAAAACTGCC‑Agggggg                                                                                                                                                     >  1:1139536/1‑149 (MQ=255)
 cctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAAACCGCC‑AGGGGGGGcc                                                                                                                                                  >  1:223312/1‑151 (MQ=255)
   tGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GGACGT‑AAAACCTTCA‑AAGGGGGGCCgg                                                                                                                                                >  1:693665/1‑151 (MQ=255)
    ggTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCgg                                                                                                                                                <  1:231406/150‑1 (MQ=255)
    ggTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCgg                                                                                                                                                <  1:181296/150‑1 (MQ=255)
     gTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGAt                                                                                                                                              <  1:123339/151‑1 (MQ=255)
     gTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAAACTTCA‑AGGGGGGGCCgg                                                                                                                                                >  1:448007/1‑149 (MQ=255)
       aTTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATa                                                                                                                                             <  2:531716/150‑1 (MQ=255)
       aTTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGAGACCAGATa                                                                                                                                             <  2:608130/150‑1 (MQ=255)
         tCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAag                                                                                                                                           <  1:1146695/150‑1 (MQ=255)
          ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAagt                                                                                                                                          <  1:775188/150‑1 (MQ=255)
          ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGG‑AAAAACCGCA‑AGGGGGGACCgggaaa                                                                                                                                            >  1:257644/1‑144 (MQ=255)
           cAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GGGCGT‑AAAACCCGCA‑AGGTGGGACCggggaagg                                                                                                                                          >  2:544935/1‑143 (MQ=255)
           cAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GAACGG‑AAAACCCGCA‑AGGGGGGGCCGGa                                                                                                                                               >  2:665378/1‑144 (MQ=255)
            aGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑TAAACCTGCC‑AGGGGGGGCCGGag                                                                                                                                              >  1:686505/1‑143 (MQ=255)
             gTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GGACGT‑AAAAACCGCC‑AGGGGGGGCCCGGTAAAGAGTc                                                                                                                                      >  1:93247/1‑151 (MQ=255)
              tttGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAAACTGCA‑AGGGGGGACCGGATAagta                                                                                                                                         >  2:927557/1‑147 (MQ=255)
                 gATGTTATCGGCAATACACAGCGCAGTCATCAGTGGCTGGACGGTTTTTACCAGCGGATCCACGAGGTAGTAGTAACGGCCATCCTCGCGGCCATCGAGCGACTGGCGG‑GTACGT‑AAAACCAGCA‑AGGTGGGACCGGATAAGTAGTCa                                                                                                                                     <  2:375139/148‑1 (MQ=255)
                   tGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAAACTGCA‑AAGGGGGGCCGGGTAAGTAGTaaa                                                                                                                                    >  1:993576/1‑147 (MQ=255)
                   tGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTAAGG‑AAAAACTGCA‑AAGGGGGGCCGGGTaaaatagcaaac                                                                                                                                  >  2:832106/1‑139 (MQ=255)
                         cGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAAACTGCA‑AGGGGGGGCCGGGTAAATAGGCAACCTGTGAcc                                                                                                                           >  1:1194287/1‑150 (MQ=255)
                                 cacaGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGG‑AAAAACCACA‑AAGGGGGGCCGGGGAAGGAGGTAAACTGTGGCCActcct                                                                                                                     >  2:871/1‑146 (MQ=255)
                                      cgcATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGAca                                                                                                             <  2:946322/151‑1 (MQ=255)
                                      cgcATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTCGCGG‑GTACGT‑AAATCCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGAca                                                                                                             <  2:152389/151‑1 (MQ=255)
                                         aTTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGAcatca                                                                                                          <  2:1150175/151‑1 (MQ=255)
                                         aTTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGGGGGGACGGAAAAGTAGTCAACCTGTGACCATTCTTCCCTGagcac                                                                                                           >  1:1026960/1‑146 (MQ=255)
                                             atcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GGACGT‑AAAACCTGCA‑AGGGGGGACCGGGGAAAGAGTCAACCTGGGACCCTTATTCCCCGAAATCATcc                                                                                                       >  2:546758/1‑150 (MQ=255)
                                              tcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GGAACT‑AAAAAACGCA‑AGGAGGGGCCCGGTAAATAGTCAACCTGTGGACACTCCTCGCtcgaataatac                                                                                                       >  2:1343184/1‑139 (MQ=255)
                                               caGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GGAAGT‑AAAAACTGCC‑AAGGGGGGCCGGATAAAAAGTCAACCCGTGACCCTTCTTTTCTGAaatcatacctg                                                                                                    >  1:641744/1‑146 (MQ=255)
                                                aGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTAACT‑AAAAACCTCA‑AGGGGGGGCCGGGGAAAAAGGCAAACTGTGGACATTTTTCGCTTACAACATCCtcg                                                                                                    >  1:1185168/1‑148 (MQ=255)
                                                    gCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAAACTGCC‑AGGGGGGGCCGGATAAATAGTCAACCTGTGACCATTTTTCCCTGACATCATCCtggggggt                                                                                               >  1:317068/1‑149 (MQ=255)
                                                      tGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GGACGG‑AAAAACCTCA‑AGGGGGGAACCGCAAAGTAGTCAAACTGTGGACATTCTTCGCCGACCTCCACCttgtgg                                                                                                 >  2:887530/1‑147 (MQ=255)
                                                         aCGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GGAAGTAAAAACCTCCA‑AGG‑GGGGCCGGATAAGAAGGCCACCTGTGACCATTCTTCTCTGACATAATCCTGGGGGGGGAcc                                                                                           >  2:1157438/1‑150 (MQ=255)
                                                          cGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACAGGCGG‑GTACGT‑AAAAACTGCA‑AGGTGGGACCGGATAAGGAGTCAACCTGTGACCATTCTTCGCTGACAACATTCTGGTGGCTGACca                                                                                          >  2:1347609/1‑150 (MQ=255)
                                                           ggTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAAACTGCA‑AGGGGGGACCGGGAAAATAGACAACCTGTGACCATTCTTTGCTGTCATCCTCCTGGGGGGTGAcc                                                                                           >  1:1309814/1‑148 (MQ=255)
                                                             ttttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACTT‑AAAACCTGCA‑AAGGGGGGCCCGGGAAGTAGTTAAACTGTGACCAATTTTTGCCGACATCATCCTGGGGGCTGACcaca                                                                                        >  1:263235/1‑149 (MQ=255)
                                                                tttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑GAAAACTGCA‑AGGTGGGGCCGGATAAGGAGTCAACCGGTGAACATTCTTCGCTGACATCATCCTGGGGGGTGACca                                                                                          >  2:706355/1‑144 (MQ=18)
                                                                     tGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTAAGG‑AAAAACCGCA‑AGGGGGGACCGGGGAAATAGTCAACCTGTGACCCTTCTTCGCTGACATCATCCTGGTGGGTGACCACCTCAC‑TTaaa                                                                              >  1:703067/1‑150 (MQ=21)
                                                                      gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGaaa                                                                              <  2:883343/149‑1 (MQ=255)
                                                                      gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAAcg                                                                            <  1:1034829/151‑1 (MQ=255)
                                                                      gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAAACTGCA‑AGGGGGGGCCGGATAAATAGTCAACCTGTGACCATTCTTCGCTTACATCATCCTGGTGGGTGACCCCATCAC‑TGaaa                                                                              >  1:1196353/1‑149 (MQ=255)
                                                                       cATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAAcgc                                                                           <  1:202434/151‑1 (MQ=255)
                                                                          atCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAAACTGCC‑AGGGGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCggg                                                                        >  1:176322/1‑150 (MQ=255)
                                                                          atCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGC‑AAAACCTGCA‑AGGGGGGGGCGGGGAAGGAGACCACCTGTGACCATTCTTCGCTGACATTATCCTCGTGGGTGACCCCAACAC‑TGAAACGCgg                                                                         >  1:949397/1‑150 (MQ=18)
                                                                           tCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GGACGT‑AAAACCCCCA‑AAGTGGGGCCGGAGAAGTAGTCAACCTGTGACCATTCTTTCCTGACATCATCCTTGTGGCTGACCACATCCC‑TGAAAAGCgg                                                                         >  2:389471/1‑149 (MQ=255)
                                                                            ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTgg                                                                      <  1:724168/151‑1 (MQ=255)
                                                                             cACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GGACGT‑AAAACCTGCA‑AGGGGGGGCCGGGTAAGTAGTCAAACTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGGCCCCATCAC‑TGAAACGCgggggg                                                                     >  1:1259326/1‑151 (MQ=21)
                                                                                gAGGTAGTAGTAACTGCCATCCGCGCAGCCATCGGGCGACTGGCGGCTTAGGC‑AAAACC‑CCA‑ACGGGGGGCCGGGATAGTAGGCAAACTGTGACCAATCCTCCCCGACAACATCCTGGTGGTGGTCCCCATCCC‑TTAAAAACTGTgggggg                                                                  >  2:365829/1‑149 (MQ=255)
                                                                                        agtaACTGCCATCCTCGCAGCCATCGGGCGACTGACGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCtt                                                          <  2:1202081/151‑1 (MQ=255)
                                                                                           aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCTTaaa                                                       <  1:1183613/151‑1 (MQ=255)
                                                                                           aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCTTaa                                                        <  2:1203849/150‑1 (MQ=255)
                                                                                           aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCTTaa                                                        <  1:992729/150‑1 (MQ=255)
                                                                                           aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGG‑AAAAACTGCA‑AGGGGGGGCCGGATAAAGAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGGGGCTGACCACATCAC‑TGAAACGCGGGGGGGGGGGGGCCttta                                                        >  1:682008/1‑148 (MQ=18)
                                                                                                ccATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCCCAACAGGGGGGGCC‑GCTAAAGAGTCAAACCGTGGCCCCTTTTTCCTGACCTCATCCTGGGGGCTGACCACA‑CACACGTAACGCGGGGGGGGGGGGGGCTTAAAatt                                                    >  2:894691/1‑149 (MQ=255)
                                                                                                   tCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGC‑AAAACCTGCA‑AGGTGGGGCCGGCTAAGGAGTCAACCTGTGGCCATTCCTCCCTGACATCATCCCGGGGGCGGACCACATCAC‑TGGAACGCGGGGGGGGGAGGGCCCTTAATTTTCt                                                 >  2:237363/1‑149 (MQ=255)
                                                                                                     cTCGCAGCCATCGGGCGACTGGCGG‑GGACGT‑AAAAACTGCA‑AGGGGGGGCCGGGTAAGTAGTCAACCCGTGACCATTCTTTGCTGACATCATCCCTGGGGCTGACCACATCAC‑TGAAAAACGGGGGGGGGGGGGCCTTAAATTTACTGAt                                              >  1:478886/1‑150 (MQ=255)
                                                                                                       cGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATg                                             <  1:412626/149‑1 (MQ=255)
                                                                                                       cGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATg                                             <  2:814465/149‑1 (MQ=255)
                                                                                                         cAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCt                                         <  1:392724/151‑1 (MQ=255)
                                                                                                          aGCCATCGGGCGACTGGCGG‑GTACGT‑AAAAACTGCA‑AGGGGGGACCGGGTAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGGTGACCACATCAC‑TGAAACGCGGGGGGGGGGGGGCCTTAAATTTTCTGATGGCCtt                                        >  1:1223653/1‑151 (MQ=18)
                                                                                                            ccATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCGCATCTC‑TGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCa                                      <  2:1039310/151‑1 (MQ=255)
                                                                                                               tCGGGCGACTGGCGG‑GTGCGT‑AAAAACTGCA‑AGGTGGGACCGGATAAGAAGTCAACCTGTGACCCTTCTTCGCTGACATCATCCTGGGGGCTGGACACATCAC‑TGAAAAGCGGGGGGGGGAGGTGCTTTAATTTACTGATGGCCTTCAgg                                    >  2:1274933/1‑150 (MQ=255)
                                                                                                                cGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGt                                   <  1:1166831/150‑1 (MQ=255)
                                                                                                                 gggCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtca                                 <  1:480041/151‑1 (MQ=255)
                                                                                                                                 aCGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGACCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCg                  >  1:522552/1‑151 (MQ=255)
                                                                                                                                   gT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCa                <  1:1116019/151‑1 (MQ=255)
                                                                                                                                    t‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCa                <  1:383094/150‑1 (MQ=255)
                                                                                                                                    t‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCAt               >  1:381638/1‑151 (MQ=255)
                                                                                                                                      aaaaCCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATg              <  1:1203656/151‑1 (MQ=255)
                                                                                                                                             gCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATg         <  1:921877/149‑1 (MQ=255)
                                                                                                                                                 aGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGaa     >  1:587068/1‑150 (MQ=255)
                                                                                                                                                   gTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACa  <  2:785758/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGG‑GTACGT‑AAAACCTGCA‑AGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCAC‑TGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACA  >  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: ‑