Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,621,1180AC38.0% 75.7 / 12.3 51V7G (GTG→GGG) eamAcysteine/O‑acetylserine exporter EamA
Reads supporting (aligned to +/- strand):  ref base A (21/10);  new base C (0/19);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 7.40e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATCTCTAAAACAAAACATAACGAAACGCACTGCCGGACAGACAAATGAACTTATCCCTACGACGCTCTACCAGCGCC  >  NC_000913/1620875‑1621364
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
aGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATc                                                                                                                                                                                                                                                 >  1:8798/1‑251 (MQ=255)
   ccAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCtt                                                                                                                                                                                                                                               <  1:575665/250‑1 (MQ=255)
       cagGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTcgc                                                                                                                                                                                                                                          >  1:62383/1‑251 (MQ=255)
             tACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCAGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACAtt                                                                                                                                                                                                                                    >  1:20895/1‑251 (MQ=255)
                   aGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGTGACATTCTACCg                                                                                                                                                                                                                              <  1:436315/251‑1 (MQ=255)
                           gCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCAAGCAGCAAATGCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGaa                                                                                                                                                                                                                      <  1:38812/251‑1 (MQ=255)
                              caAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAAttt                                                                                                                                                                                                                   >  2:191557/1‑251 (MQ=255)
                                 aaaagaaaagCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCAt                                                                                                                                                                                                                >  2:605366/1‑251 (MQ=255)
                                  aaagaaaagCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATc                                                                                                                                                                                                               >  1:379541/1‑251 (MQ=255)
                                      aaaagCAAACTGCGCAAAAATGATGTTTAATCCATACCTCAGCAGCAAATTCAGTGGAACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATCATGCCCCACTTTGATGACCACAACATTTAGCCCCCATACGACCACTCCCAGTCGCGCCAACACCCCATCTTTTCCCGACACTCTACCTCCTCTGAATTTCATCtttt                                                                                                                                                                                                           >  2:332719/1‑251 (MQ=255)
                                        aagaaaaCGGCGAAAAACTGATGGTTAATCCATACCCCAGCAGCAATTTCAGTGGTACTTGCGGTCGTGCGAAAAAAAAGAAAGCCGGAAAAGCGACCAGCAAAAAGCGCAAACCGGCCAGCATCAGCGGGGGCATGTTATAAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCAGACGACCACTACCAGGAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGt                                                                                                                                                                                                         <  1:492680/247‑1 (MQ=255)
                                                 gcgcAAAACTGATGGTTAATCCATACCCCAGCAGAAAATTCAGTGGTACTTCCGGTCGTGCGAAAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCAATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAt                                                                                                                                                                                                  <  1:266896/249‑1 (MQ=255)
                                                        actagttgttattccaccccccAGCAGCAAATTCAGTGTTCCTTCCGCTCGGGCAAAAAAAAAGATAGCGGAAAAAGGCACCAGAATAAAGCGAAACCCGGCCGGCATCAGCGGTGGCTTTTTAGGAGCCCCCACTTTTACGCCCACAAAATTTACCCCCCATACACCCACTCCCGGGAGCGCCAACCCCCCATCTTTCCGCAACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGc                                                                                                                                                                                         <  2:66678/233‑1 (MQ=255)
                                                            aTGGTTAATCCATCCCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATCCGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGaa                                                                                                                                                                                     <  1:605366/251‑1 (MQ=255)
                                                            aTGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTACCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCACCTTAGCTGa                                                                                                                                                                                      >  1:161702/1‑250 (MQ=255)
                                                                  aaccactACCCCACCAGCAAATGCAGGGGTACTGTCGGTCGTGACACAAAAAAGATAGCCGGAAAAGCCCCCAGCAAAAAGCCAAACCCGCCCAGCTTCAGCGTGGGCATTTTATGATTCCCCACTTTGATGCCCACAAAATTTACCCCCCATACGACCCCCCCCAGTAGCTCCAACACCCCATTTTTCCGCGACATCTTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTAc                                                                                                                                                                                <  2:520998/244‑1 (MQ=255)
                                                                  aaccaaTACCCCAGCAGCAAATCCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATTTTATGAACCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTAc                                                                                                                                                                                <  1:386410/245‑1 (MQ=255)
                                                                            ccAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAAc                                                                                                                                                                     >  1:391335/1‑251 (MQ=255)
                                                                            ccAGCAGCAAAGTCAGTGTTACTTTCGGTCGTGCAAAAAAAAAGATACCCGGAAAAGAGCCCAACATAAAGCGCAACCCGGCCAGCATCAGCGGTGGAATGTTATAAAGCCCCTCTTTGATGCCCAAAAAATTGCGCCCCCATACGACCACTACCATTAGCCCCAACTCCCCATCTTTTCGCCACGTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAAc                                                                                                                                                                     <  2:108906/251‑1 (MQ=255)
                                                                              agcagcAAATTCAGTGGTCCTTCCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAg                                                                                                                                                                   <  1:191557/251‑1 (MQ=255)
                                                                                cagaaaaTTCAGTGTTCCTTCCGGTCGGGCGAAAAAAAAAAGAGCCGGAAAAGCGACCAGCATAAGGCGCAACCCGGCCCGCATCAGCGGTGGCATGTTAGAAAGCCCCACTTTGATGCCCAAAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCTTCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGtt                                                                                                                                                                 <  2:462068/247‑1 (MQ=255)
                                                                                aagcaAATACAGTGGTCCTTACGGGCGTGCGAAAAAAAAACAAGCCGGAAAAGCGACGAGCAAAAAGCCCAAACCGCCCAGCATCAGCGGTGGCTTTTTATGACGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCCCTACCTGTAGCGCCACCCCCCCATTTTTCCCCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGtt                                                                                                                                                                 <  2:173098/250‑1 (MQ=255)
                                                                                     atttaaGTGGTACTTTCGGTCGTGCAAAAAAAAATATAGCGGGAAACCGGCCCAGCATAAAGCGCAACCCGCCCATATTCAGCGGTGGCATGTTATGAACCCCCCCTTTGATGACCAAAAAATTTACCCCCCATCCCCCCCCTACCAGTAGCGCCACCCCCCCTTTTTTCCGTGAATTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTc                                                                                                                                                            <  2:245484/246‑1 (MQ=255)
                                                                                     aaTTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCtt                                                                                                                                                                                                                                               <  2:354424/168‑1 (MQ=255)
                                                                                     aaTTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCtt                                                                                                                                                                                                                                               >  1:354424/1‑168 (MQ=255)
                                                                                                aCTTTCGGTCGTGCGAAAAAAAAGATAGCCGGAAAAGCACCCACAATAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGCCCAAAAAATTTACCCCCCATACGCCCACTACCAGTAGCCCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCg                                                                                                                                                 <  1:14608/251‑1 (MQ=255)
                                                                                                   ttCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCAAAAAGCGAAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCATCCGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGt                                                                                                                                               <  1:590934/250‑1 (MQ=255)
                                                                                                     cGGTCGGGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCAAAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATTTTATGAACCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTac                                                                                                                                            <  1:386552/251‑1 (MQ=255)
                                                                                                       gTCGTGCAACAAAAAAGAGAGCCGGAAAAGCGACCAGCATAAAGCGAAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTAATGCCCAAAAAATTTAGCCCCCATACGACCCCTACCGGGAGCGCCAACCCCCCATCTTTCCGCGCCATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacga                                                                                                                                          <  1:153046/251‑1 (MQ=255)
                                                                                                                    aaaGATAGCCAAAAAAGCGACAACCAAAAGGCGCAACCGGCCAAGCATCAGCGGTGGCATTTTATGAAGCCCCACTTTGATGACCACAAAATTTACCCCCCATCCGACCACTCCCAGGAGCGCCAACCCCCCATCTTTGCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgcg                                                                                                                             <  2:20895/251‑1 (MQ=255)
                                                                                                                      aGATGGCCGGAAAAGCGCCCAAAAAAAAGCGCAACCCGGCAAGAACCAGGGGTGGCATGTTTGGAAGCCCCCTTTTGATGACCACAAAATTTAGCCCCCATCCGACCACTACCAGTAGCCCCAACCCCCCATCTTTCCGCAACATCATCCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCa                                                                                                                           <  2:306580/251‑1 (MQ=255)
                                                                                                                              ggAAAAGCGACAAGCAAAAAGCGAAAACCGGCCAGCATCAGCGGTGGCATGTTACGAAGCCCCACTTTGAGCACCACAAAATTTAGCCCCCATACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTc                                                                                                                   <  2:432442/251‑1 (MQ=255)
                                                                                                                                 aaaGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCt                                                                                                                >  2:221132/1‑251 (MQ=255)
                                                                                                                                        ccAGCAAAAGGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCATTTTGATCACCACAAAATTTACCCCCCATCCGCCCACTCCCAGTAGCGCCAACCCCCCATCTTTCCGCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGGCGGCATTAATGCGCAAATAAGTCGCTATACTTc                                                                                                         <  2:8798/251‑1 (MQ=255)
                                                                                                                                                          cGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                                                                                                                              >  2:231756/1‑196 (MQ=255)
                                                                                                                                                          cGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                                                                                                                              <  1:231756/196‑1 (MQ=255)
                                                                                                                                                               cagcagcagGGGGGGCATGTTAGGACGCCCCACTTTGAAGACCAAAAAATTTAGCCCCCAAACGACCCCCCCCAGTAGCGCCACCCCCCCATCTTTCCTCACCATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCt                                                                                   <  2:493199/251‑1 (MQ=255)
                                                                                                                                                                cacaaTCAGCGTTGGCATGTTAGAACCCCCCTCTTGGAGAACCACAAAATTTCCCCCCCATCCGACCACTACCAGTGTCGCCAACACCCCATTTTTTCGCCACATTCTCCCGCCTCTGAATTTCATCTTTTGTAAGCAAGCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCt                                                                                   <  2:179260/247‑1 (MQ=255)
                                                                                                                                                                  gCATCAGCGGTGGCATTTTATGAAGCCCCACTTTGATGACCACAAAATTTACCCCCCATCCGACCACCACCAGTAGCGCCAACCCCCCATCTTTTCGCAACATCCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCttt                                                                                 <  1:47578/251‑1 (MQ=255)
                                                                                                                                                                  gCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCttt                                                                                 >  1:169815/1‑251 (MQ=255)
                                                                                                                                                                  gCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCttt                                                                                 >  2:237957/1‑251 (MQ=255)
                                                                                                                                                                    ataagtggTGGCATTTTATGACACCCCCTTTTGATGACCAAAAATTTTAGCCCCCATACGCCCACTCCCAGTAGCCCCAACACCCCATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTAc                                                                               <  2:39545/245‑1 (MQ=255)
                                                                                                                                                                               cATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATCTCTaaaac                                                                    >  1:401002/1‑251 (MQ=255)
                                                                                                                                                                                aTGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATCTCTaaaac                                                                    >  1:454493/1‑250 (MQ=255)
                                                                                                                                                                                     aTGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTACCAGTTGATTCGTTAGTCGCCGGTTACGACTGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATCTCTaaaacaaaac                                                               >  2:244836/1‑250 (MQ=255)
                                                                                                                                                                                              cACTTTGAAGACCAAAAAATTTATCCCCCATCCCCCCACTCCCAGTAGCGCCAACCCCCCTTCTTTTCTCAACATTCAACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATCTCTAAAACAAAACATAACGAAAc                                                     <  2:355633/251‑1 (MQ=255)
                                                                                                                                                                                                            acaAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCACCACCCCATCTTTTCGCGACATTCTACAGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATCTCTAAAACAAAACATAACGAAACGCACTGCCGgac                                         >  1:537730/1‑250 (MQ=255)
                                                                                                                                                                                                                     aGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATCTCTAAAACAAAACATAACGAAACGCACTGCCGGACAGACAAATg                                >  2:451326/1‑250 (MQ=255)
                                                                                                                                                                                                                          ccATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATCTCTAAAACAAAACATAACGAAACGCACTGCCGGACAGACAAATGAACtt                           >  1:7482/1‑250 (MQ=255)
                                                                                                                                                                                                                               cgccccctcccAGTAGCCCCACCCCCCCATTTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATCTCTAAAACAAAACATAACGAAACGCACTGCCGGACAGACAAATGAACTTATCCCt                     <  2:391335/242‑1 (MQ=255)
                                                                                                                                                                                                                                                  cAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTACGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGCTTTTTGCCATGCTATTTCTTTACATCTCTAAAACAAAACATAACGCAACGCACTGCCGGACAGACAAATGAACTTATCCCTACGACGCTCTACCAGCGcc  >  1:334728/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
AGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTTTGCCATGCTATTTCTTTACATCTCTAAAACAAAACATAACGAAACGCACTGCCGGACAGACAAATGAACTTATCCCTACGACGCTCTACCAGCGCC  >  NC_000913/1620875‑1621364

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