Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,6260GA24.5% 80.5 / 15.0 49F180F (TTC→TTTycfTinner membrane protein YcfT
Reads supporting (aligned to +/- strand):  ref base G (28/9);  new base A (0/12);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 3.19e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACTTAACGCCAGCCACTGCAC  >  NC_000913/1174376‑1174868
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
tACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACg                                                                                                                                                                                                                                                       >  1:372369/1‑251 (MQ=255)
tACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACg                                                                                                                                                                                                                                                       >  1:350796/1‑251 (MQ=255)
   cGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGaa                                                                                                                                                                                                                                                     >  1:391090/1‑250 (MQ=255)
    gAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAAt                                                                                                                                                                                                                                                   <  2:348654/251‑1 (MQ=255)
                   ttGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCAACGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCTACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCACGCGCCAAGGCTCTAATAAGGCAAATTGAGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACGGCCACACTc                                                                                                                                                                                                                                     >  2:370154/1‑250 (MQ=255)
                              gcttcagaatacccaaaaccGATCCCAGCGACAACGGCAGCGAGAAAATAAAAAACCAGCCCCCGCCCACCAGAACGGTAAGAAAAGAAGCCATCAGCAGAGGGGGCCGCACCGGCACCTTTTTAACAAGGGTATTTTTTGTTGCGCCAAACCATGCGCCAGGGTTTAAATAGGGAAAATTGCGGATCACACGTTTCTTCCCCCACCACGGGGGGGAACCAAAATTAACCGCCACACTCAGCAGTacaaac                                                                                                                                                                                                                         <  1:154196/232‑1 (MQ=255)
                                  cATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAACAAACAACCAGGCACCGACCGCCAGAACGGAAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGAGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGACATTAACCGCCACACTCTGCAGTACAAACAAgg                                                                                                                                                                                                                    >  1:103562/1‑251 (MQ=255)
                                               aTCGATACCAGCGAAAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATGGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAAGTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCTGCAGTACAAACAAGGCGAATAGCGGc                                                                                                                                                                                                         >  2:275607/1‑249 (MQ=255)
                                                   aaaCCAGCGACAACAGAAGAGAGATATTAAACACCCGGGCCCCGACCGCCAGAACGGTAAGAAAAAACGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACCCAGGTCATTATTGTTGCCCAAAACCATGCCCCAAGGTTTTAATAAGGCAAATGGCGGGTCCCCCTGTTCATTCCCCCCCCCGGGGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                                                                                                                                                                   <  1:346920/249‑1 (MQ=255)
                                                      ccAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGCACGAAATTAACCGCCACACTCAGCAGTACAAACCAGGCGAATAGTGGCAGCGcc                                                                                                                                                                                                   >  1:369765/1‑248 (MQ=255)
                                                             caacaGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTa                                                                                                                                                                                          >  1:26958/1‑250 (MQ=255)
                                                                cagaagagagaTAGAAAACACCCAGGCCCCGCCCGCAAGAACGGTAAGAAAAGAACCCATCAGCAGATGGCGGCGAAACGGCCCCTCTTTAACACAGGTCATTATTGTGGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACCCTGTTCATCCCCCCCCCCGGCGGGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaaaa                                                                                                                                                                                      <  1:55484/244‑1 (MQ=255)
                                                                 agcagcGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAt                                                                                                                                                                                     >  1:448926/1‑251 (MQ=255)
                                                                   cagcGAGATAGTAAACAACCAGGCACCGAACGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAttt                                                                                                                                                                                   >  1:310326/1‑251 (MQ=255)
                                                                       gagaTAGGAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTAc                                                                                                                                                                                >  1:241590/1‑250 (MQ=255)
                                                                       gaagaaaaaaaCACCCAGGCACCGACCGCCAGAACGGTAAAAAAAAACGCAAACAGCGGTGGGCGGCGCAACGGCCCCTCTTAAACACAGGCAATTATTGTTGCGCAAACACATGCGCCAGGGCGAAAAAAAGGCAATTGGCGGATCACACTTTTCATCCCCCCCCACGGGGTGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACa                                                                                                                                                                               <  2:160536/243‑1 (MQ=255)
                                                                                   aacCAGGCACCGACAGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCCGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTACTAAGGCAACTTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACCCAGCAGTACAACAACGGCGAATAGTGGCAGCGCCAGACGCCTCAAAACTTTACATACCACGAAAt                                                                                                                                                                    >  2:207128/1‑250 (MQ=255)
                                                                                     ccAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGa                                                                                                                                                                                                                                                      >  1:148483/1‑167 (MQ=255)
                                                                                     ccAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCGTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGa                                                                                                                                                                                                                                                      <  2:148483/167‑1 (MQ=255)
                                                                                        ggCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGCAATAGACa                                                                                                                                                               >  1:266586/1‑250 (MQ=255)
                                                                                          cACCGACCGCCAGAACGGTAAGCAAAGAAGCAATCCGCAGATGGCGGCGCAACGGCACCTCTTTAAAACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAAtt                                                                                                                                                            >  1:336048/1‑251 (MQ=255)
                                                                                                   cgaGAACGGAAAGCAAAGAACCTAGCAGAAGATGGCGGCGCACCGGCCCCTTTTTAACGATGGTCATTATTGTGGCGCCAAACCAGGCGCAAGGGGTGTAATAAGGCAAATTGGGGATCACGCTTTTCATCCCCCACCCCGGGGGGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATAc                                                                                                                                                   <  2:449015/249‑1 (MQ=255)
                                                                                                     gtaagggaaggaaaaGAACCCACCAGCGGAGGTCGGCGCAACGGCACCTATTAACCAGGGGTCTTTATTGTTGCGCAAAACAAGTGGCAAAGTGTGTAATAAGGGAAATCGCGGACCCACTTGTTCATCCCCCACCGCGGGGGGGAAACAAAATTAGCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAg                                                                                                                                                 <  1:394360/239‑1 (MQ=255)
                                                                                                              aGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCACTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATCAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAATTTACATACCACGAAATAGACAATTAACGCATACAGATACcaca                                                                                                                                         >  2:298652/1‑250 (MQ=255)
                                                                                                                            ccatcagtgggcggcccagtgcacccTCTATAACACAGGACTTTTTTGTTCCCCCAAACAAGGCCCAAGGGCTGTAAAAAGAAAAATTGGGGTGCCCACTGTTCCTTCCCCACCGCGCTGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAAGTAACGCATACAGATACCACAAGCTGGTGCTGGCg                                                                                                                           <  2:164208/227‑1 (MQ=255)
                                                                                                                            cagcagAATGCGACGCAACGGCACCTCTTTAACACATGTCATTATTGTTGCGCCAAACCCTGCGCCACGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACAGCCACACTCAGCAGTACAACCACGGCGAATAGAGACAGCCCCAAACGGAAAAAAATTTAACCCACCACGCAATAGTAGATGAACGCATACAAATACCAAAAGcatgaccaggatg                                                                                                                          >  2:267163/1‑239 (MQ=255)
                                                                                                                              gcagATGGGGGCGCACGTGCACCTTTTAAACACAGGTAATGATGTTGGCGCTAACACATGCGCCAGGGTTGTAAGAAGGCAAATGGGGGATCACACTCTCTATCCCCCACCACGGGGGGGAAACGAAATAAGCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTg                                                                                                                        <  2:84939/251‑1 (MQ=255)
                                                                                                                                   tGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGatcat                                                                                                                   >  1:102895/1‑251 (MQ=255)
                                                                                                                                    gggcggctcaccgcccccTCTTAAACACAGTTCATTTTGGGTGAGCCAAACCACGCGCCAGGGTTTAAAGAAGGCAAATGGCGGATCACCCTGTTCATCCCCCCCCACGCGGGGGAAACGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAGTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGatcat                                                                                                                   <  2:102863/235‑1 (MQ=255)
                                                                                                                                           cAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGgtcaactcgcgc                                                                                                            >  2:98844/1‑242 (MQ=255)
                                                                                                                                             aCGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGatgctccccctccgg                                                                                                         >  2:112025/1‑238 (MQ=255)
                                                                                                                                             aCGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAgg                                                                                                        >  2:291410/1‑251 (MQ=255)
                                                                                                                                                         ttAACAGAGGTTATTTTTGTTGCCACAAACCAGGCGCAAAGGCTGTAAAAAGGAAAATGGCGGACCACTCTGTTCTTTCCCCCCCCCGGCGGGGGAAGAAAATTAGCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCgg                                                                                              <  2:179088/250‑1 (MQ=255)
                                                                                                                                                          aaccacgggtttttatgGTTGCGCCAACCCGTGCGCAAAGTTTGTAAGAAGGAAAATGGCGGATCACCCTTTCCATCCCCCCCCCGGGCGTGGGAAAAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTg                                                                                            <  1:334182/239‑1 (MQ=255)
                                                                                                                                                                          tgttgCGCCAAACCATGCGCCAAGGCTGTAAAAAGGCAAATTGCGGATCACCCTGTTCATTCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCg                                                                             <  1:291410/250‑1 (MQ=255)
                                                                                                                                                                           gttgCGCCAACCCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGc                                                                           >  1:288363/1‑251 (MQ=255)
                                                                                                                                                                                              cgcCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCtt                                                          >  1:288491/1‑251 (MQ=255)
                                                                                                                                                                                              cgcCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGGGGAATCGGCATAAGCGGCATTGGAGGCATTGCtt                                                          >  1:266472/1‑251 (MQ=255)
                                                                                                                                                                                                 cAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAACCAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATATCACCAGCTGGTGCTGGCGGTGAACCTCCCGTGCAGGAACCCCCCCGGGGAATCGGCCTAAGTGGCAATGGAGGCATTGCTTaa                                                        >  1:111391/1‑250 (MQ=255)
                                                                                                                                                                                                   aGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGGGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATc                                                     >  1:259042/1‑251 (MQ=255)
                                                                                                                                                                                                         gTAATAAGGCAAATTGCGGATCACCCTGTTCATTCCCCCCCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCt                                                <  1:66238/251‑1 (MQ=255)
                                                                                                                                                                                                           aataaGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCCTCCCGTGCAGGAACACACCCGGGGAAACGGCCCAAACGGCAATTGGGGCATTTCCTTAATCGCGCTcc                                              >  2:183761/1‑250 (MQ=255)
                                                                                                                                                                                                                ggCAAATTGCGGATCCCACTGTTCATCCCCCCCCACGGGGTGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTgc                                         <  1:67318/251‑1 (MQ=255)
                                                                                                                                                                                                                     gttGCGAATCACACTTTGCAGCCCCCCCCCCGGGGTGGGCACAAAATAACCCGCCACACGCAGCAGTACAAACAAGGCGAATAGTGGCGGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACATATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATGGGAGGCATTGCTGAAATCGCGCTCATGTgccagcc                                    <  2:338937/250‑1 (MQ=255)
                                                                                                                                                                                                                       tGCGGATCCCTCTGTTCTTTCCCCCCCCCGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTgccagcca                                   <  1:275607/250‑1 (MQ=255)
                                                                                                                                                                                                                        gCGGATCACACTGTTCATTCCCCACCACGGCGTGCGAACGACATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCCGAGCCAGACGGCTAAAAATTTTACATAACACGAAATAGACAATTACAGCCTACAGATACCCCAAGCTGGTGCTGGCGGGGATCATCCCCGGCAAGGAAGAAACGGGGGAATCGGCACAACCGGCATTGGAGGCATTGCTTAACTCGCGCCCAGGGgccacacaac                                 >  2:411338/1‑245 (MQ=255)
                                                                                                                                                                                                                                       ccttcccccACCGCGGCGTGGAAAAGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACt                    <  2:350796/245‑1 (MQ=255)
                                                                                                                                                                                                                                       caatccccccccggggggTGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACtt                   <  2:102895/234‑1 (MQ=255)
                                                                                                                                                                                                                                              ccaccggggggTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACTTAACgcca            <  1:411338/241‑1 (MQ=255)
                                                                                                                                                                                                                                              ccaccaCGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAACGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACTTAACgcc             >  1:167061/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        gTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGGATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGGTCAGCGCACTTAACGCCAGCCACTGcac  >  1:341884/1‑251 (MQ=255)
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
TACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACTTAACGCCAGCCACTGCAC  >  NC_000913/1174376‑1174868

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: