Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,928,3430GA33.7% 101.9 / 75.9 98L39L (CTG→CTAsdaCputative serine transporter
Reads supporting (aligned to +/- strand):  ref base G (35/30);  new base A (16/17);  total (51/47)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTATTTGCATTCCTCCAGGAGAAATAGATGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACTGTCTGGTAAAAACCCGGGCGAAGACATCACCGAGGTTG  >  NC_000913/2928200‑2928491
                                                                                                                                               |                                                                                                                                                    
ttattTGCATTCCTCCAGGAGAAATAGATGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCt                                                                                                                                                 >  1:499693/1‑149 (MQ=255)
 tattTGCATTCCTCCAGGAGAAATAGATGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTg                                                                                                                                                <  2:70730/149‑1 (MQ=255)
  attTGCATTCCTCCAGGAGAAATAGATGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCGAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGc                                                                                                                                               <  2:1334225/149‑1 (MQ=255)
   tttGCATTCCTCCAGGTGAAATAGAGGGAAACGAGTTAAACCAGCACGATTGCGTCTGAAGGCTCTCGGAGCGACTGGCGCAGGACAGACAACATGTGGATACTGGCCCTTTAAGGCACGGCGATCGGCGCGGGCGTGCTGTTCCTGcc                                                                                                                                              >  2:365094/1‑149 (MQ=255)
   tttGCATTCCTCCAGGAGAAATAGATGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGcc                                                                                                                                              <  1:1228635/149‑1 (MQ=255)
      gCATTCCTCCAGGAGAAATAGATGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAAt                                                                                                                                           <  2:615816/149‑1 (MQ=255)
      gCATTCCTCCAGGAGAAATAGATGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAAt                                                                                                                                           <  2:1333616/149‑1 (MQ=255)
        aTTCCTCCAGGAGAAATAGATGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCa                                                                                                                                         >  2:1317070/1‑149 (MQ=255)
         atcctgcCGGCGAAATACAGGGAAACCACTCAAACCACCACGATTGCGTCCAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCaa                                                                                                                                        <  2:243159/148‑1 (MQ=255)
            ctccAGGAGAAATAGATGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCATAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGc                                                                                                                                     <  2:920947/149‑1 (MQ=255)
                  gagaAATAGATGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCACTCGGCGCGGGCGTGCTATTCCTGCCAATCCACGCCGgggt                                                                                                                               >  2:122861/1‑149 (MQ=255)
                  gagaAATAGATGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGgtgt                                                                                                                               >  1:1546350/1‑149 (MQ=255)
                    gaAATAGATGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTg                                                                                                                             >  1:447207/1‑149 (MQ=255)
                    gaAATAGATGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTg                                                                                                                             <  1:274091/149‑1 (MQ=255)
                    gaAATAGATGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTg                                                                                                                             >  1:185838/1‑149 (MQ=255)
                        tAGATGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCgg                                                                                                                         <  1:707351/149‑1 (MQ=255)
                         aGATGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGt                                                                                                                        <  2:1421503/149‑1 (MQ=255)
                            tGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCAGTATg                                                                                                                     >  1:1516418/1‑149 (MQ=255)
                              gaacccaCTCAAACCAGCACGATGGCGTCGAAAGACTCTCGCAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGAt                                                                                                                   <  2:75941/143‑1 (MQ=255)
                              gAAACGACTCAAACCAGCACGATTGCGTCGCTTGCCGCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGGGTTGGCGGTATGAt                                                                                                                   >  1:759600/1‑149 (MQ=255)
                              gAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCGGTATGAt                                                                                                                   <  1:1426919/149‑1 (MQ=255)
                                aaCGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATcc                                                                                                                 >  1:1439073/1‑149 (MQ=255)
                                  cGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCCATCAACGCCGGTGTTGGCGGTATGATCCCg                                                                                                               >  1:619218/1‑149 (MQ=255)
                                  cGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCg                                                                                                               <  2:1183051/149‑1 (MQ=255)
                                    gcTCAACCCAGCACGATTGCGTCGAAAGACTCTCGGAGTGCCTGGCGCAAGACTGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGTGTGTGCTATTCCTGCCAATCCACACCGGTGTTGGCGGTATGATCCCGCt                                                                                                             <  1:873842/148‑1 (MQ=255)
                                      tCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTc                                                                                                                                                   >  1:22216/1‑109 (MQ=255)
                                      tCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTc                                                                                                                                                   <  2:22216/109‑1 (MQ=255)
                                      tCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCACCCCCGGTGTTGGCGGTATGTCCCCgtcga                                                                                                           >  2:1315955/1‑145 (MQ=255)
                                      tCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGGGTTGGCGGTATGATCCCGCTGa                                                                                                           >  1:1418969/1‑149 (MQ=255)
                                      tCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGa                                                                                                           <  1:1478258/149‑1 (MQ=255)
                                      tCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCAGGCGTGGGGTTCCGGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGa                                                                                                           >  2:1164666/1‑149 (MQ=255)
                                         aaCCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACCGACACCGTGTGGATGCTGGGCCTTTACGGCCCGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGAGCCCGCTGatca                                                                                                        >  2:240773/1‑149 (MQ=255)
                                           ccAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGGCAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGatcatc                                                                                                      >  1:1027268/1‑149 (MQ=255)
                                            cAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGAtcatca                                                                                                     >  2:1197799/1‑149 (MQ=255)
                                            cAGCACGATGGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGAtcatca                                                                                                     <  2:1401914/149‑1 (MQ=255)
                                             aGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACGCCCTATGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATcatcat                                                                                                    <  1:44267/149‑1 (MQ=255)
                                              gCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCAACTACACGTCTTTGTTGGCGGTATGATCCCGCTGATCATCATg                                                                                                   <  1:1539398/149‑1 (MQ=255)
                                                  gATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTa                                                                                               >  2:1087032/1‑149 (MQ=255)
                                                  gATTGCCTCGAAAGACTCTCGTAGTGCCTGGCGCCCGCCAGACACCAAGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGAGTTGCCGGCAACTGACCGCGGATCATCATTGCTa                                                                                               >  2:1088448/1‑149 (MQ=255)
                                                       cGTCGAAAGACTCTCGTTGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACTCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCtt                                                                                          <  1:610622/149‑1 (MQ=255)
                                                         gccaaaGGCTCTTGTAGTGCCTGGCGCAAGACGGACACCATGTGGAGGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGc                                                                                        <  1:257738/146‑1 (MQ=255)
                                                           gAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGt                                                                                      >  2:245830/1‑149 (MQ=255)
                                                           gAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTGGCGt                                                                                      >  2:1294077/1‑149 (MQ=255)
                                                            aaaGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGtt                                                                                     >  1:1062327/1‑149 (MQ=255)
                                                            aaaGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGAGGTATGATCCCGCTGATCAGCATGGCTATCCTTGCGtt                                                                                     >  1:323873/1‑149 (MQ=255)
                                                               gcctctCGGAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTccc                                                                                  <  2:476484/147‑1 (MQ=255)
                                                               gACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTccc                                                                                  >  1:670815/1‑149 (MQ=255)
                                                                  tctcGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGTCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGAt                                                                               >  1:46434/1‑149 (MQ=255)
                                                                  tctcGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATg                                                                                                   >  1:358745/1‑129 (MQ=255)
                                                                  tctcGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATg                                                                                                   <  2:358745/129‑1 (MQ=255)
                                                                        aGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGtt                                                                         <  2:764163/149‑1 (MQ=255)
                                                                        aGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGtt                                                                         >  1:79794/1‑149 (MQ=255)
                                                                          tGCCTGGCGCAAGACAGAAACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCTCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGtttt                                                                       <  2:1388541/149‑1 (MQ=255)
                                                                            ccTGGCGCAAGGCAGACACCCTGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGtttttt                                                                     <  2:424163/149‑1 (MQ=255)
                                                                            ccTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCg                                                                                                                          >  1:889016/1‑96 (MQ=255)
                                                                            ccTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCg                                                                                                                          <  2:889016/96‑1 (MQ=255)
                                                                              tGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTTCGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGATGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTGTTTTGc                                                                   >  1:1005789/1‑149 (MQ=255)
                                                                               ggCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCt                                                                  <  1:1197799/149‑1 (MQ=255)
                                                                               ggCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCt                                                                  <  1:474956/149‑1 (MQ=255)
                                                                                gcgcAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTc                                                                 <  1:1087032/149‑1 (MQ=255)
                                                                                gcgcAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTc                                                                 >  1:1088864/1‑149 (MQ=255)
                                                                                 cccAACACAGACACCATCTGGATTCTGGGCCCTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCa                                                                <  1:540000/147‑1 (MQ=255)
                                                                                  gcAAGACAGAGACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCAc                                                               <  2:1514317/149‑1 (MQ=255)
                                                                                   cAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCACCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCAcc                                                              >  1:142759/1‑149 (MQ=255)
                                                                                   cAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTAGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCAcc                                                              >  1:1180428/1‑149 (MQ=255)
                                                                                    aaGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACcg                                                             >  2:1513312/1‑149 (MQ=255)
                                                                                       acagacaCCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCgg                                                          >  1:1492903/1‑149 (MQ=255)
                                                                                       acagacaCCATGTGGATACTGGGCCTTTAAGGCAACGAAATCGGAGCGGGCGAGTTATGCTTGCCATTCATCTCCGTTTTTGGCGGTGTGATCTCGCTGTTCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCgg                                                          >  2:311454/1‑149 (MQ=255)
                                                                                         agacaCCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGcc                                                        >  2:210243/1‑149 (MQ=255)
                                                                                            cacCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGa                                                     <  2:499693/149‑1 (MQ=255)
                                                                                              ccATGGGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACt                                                   <  2:1512494/149‑1 (MQ=255)
                                                                                               cATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGTTAATCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTc                                                  >  1:752657/1‑149 (MQ=255)
                                                                                               cATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCACGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTCACTc                                                  >  2:1096754/1‑149 (MQ=255)
                                                                                                aTGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGAGCCCGCTGATCATCTTGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGGGGCCTGACTCg                                                 >  2:861582/1‑149 (MQ=255)
                                                                                                 ggtgGAGGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACAGAGGCCTTACgcgc                                                <  2:1179595/148‑1 (MQ=255)
                                                                                                  gtgGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCt                                               >  1:999504/1‑149 (MQ=255)
                                                                                                  gtgGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCGTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCt                                               <  1:331283/149‑1 (MQ=255)
                                                                                                    ggATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTTCTCTCCTCGCCCTGCCGCGCTTc                                             >  2:60286/1‑149 (MQ=255)
                                                                                                         cTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACt                                        <  1:1317070/149‑1 (MQ=255)
                                                                                                          tGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCAGCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACTg                                       >  1:570950/1‑149 (MQ=255)
                                                                                                           gggCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACTGt                                      >  1:191018/1‑149 (MQ=255)
                                                                                                           gggCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCGTGCGTTCCCGGTGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACTGt                                      >  1:1166813/1‑149 (MQ=255)
                                                                                                            ggCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTCTCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACTGTc                                     <  2:619218/149‑1 (MQ=255)
                                                                                                              ccTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTCTGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACTGTCTg                                   <  2:424877/149‑1 (MQ=255)
                                                                                                               cTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGTTCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACTGTCTgg                                  <  2:447207/149‑1 (MQ=255)
                                                                                                                tttACGGCACGGCAATCGGCGCGGGCGTGCTATTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACGGCGGCCTGACTCGCGTCGTTCTGTCTGGt                                 >  1:20890/1‑149 (MQ=255)
                                                                                                                tttACGGCACGGCAATCGGCGCGGGAGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTGCTGTCTGGt                                 >  1:775453/1‑149 (MQ=255)
                                                                                                                    cggcacggcaATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACTGTCTGGTaaaa                             <  1:1337821/149‑1 (MQ=255)
                                                                                                                        acggcaATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGGTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACTGTCTGGTAAAAcccc                         >  1:587181/1‑149 (MQ=255)
                                                                                                                        acggcaAGCGGCGCGGGCGTGCTATTCCTGCCAATCAGCGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACTGTCTGGTAAAAAccc                         <  2:302022/149‑1 (MQ=255)
                                                                                                                          ggcaATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACTGTCTGGTAAAAACCCgg                       >  1:34672/1‑149 (MQ=255)
                                                                                                                             aaTCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGATTTTCGGACTGCTCGGCCGGACTCGCTTCGTACTGTCTGGTAAAAACCCGGGCg                    >  1:120393/1‑149 (MQ=255)
                                                                                                                              aTCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACTGTCTGGTAAAAACCCGGGCGa                   <  2:1546350/149‑1 (MQ=255)
                                                                                                                               tCGGCGCGGGCGTGCTGTTCCTGCCAATCATCGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACTGTCTGGTAAAAACCCGGGCGaa                  <  1:1347079/149‑1 (MQ=255)
                                                                                                                               tCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACTGTCTGGTAAAAACCCGGGCGaa                  <  1:565224/149‑1 (MQ=255)
                                                                                                                                            gCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGCTCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTGTGCTCACCGCGGCCTGACTCGCTTCGTACTGTCTGGTAAAAACCCGGGCGAAGACATCACCGAgg     <  1:176557/149‑1 (MQ=255)
                                                                                                                                            gCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACTGTCTGGTAAAAACCCGGGCGAAGACATCACCGAgg     <  1:568941/149‑1 (MQ=255)
                                                                                                                                            gCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGACTATCCTTGCGTTCCCGATGACGTTTTTTGCTCTCCGCGGCCTGACTCGCTTCGTACTGTCTGGTAAAAACCCGGGCGAAGACTTCACCgggg     >  2:995300/1‑149 (MQ=255)
                                                                                                                                               gTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATTCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACTGTCTGGTAAAAACCCGGGCGAAGACATCACCGAGGTTg  <  2:1027268/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TTATTTGCATTCCTCCAGGAGAAATAGATGGAAACGACTCAAACCAGCACGATTGCGTCGAAAGACTCTCGTAGTGCCTGGCGCAAGACAGACACCATGTGGATGCTGGGCCTTTACGGCACGGCAATCGGCGCGGGCGTGCTGTTCCTGCCAATCAACGCCGGTGTTGGCGGTATGATCCCGCTGATCATCATGGCTATCCTTGCGTTCCCGATGACGTTTTTTGCTCACCGCGGCCTGACTCGCTTCGTACTGTCTGGTAAAAACCCGGGCGAAGACATCACCGAGGTTG  >  NC_000913/2928200‑2928491

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