Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913827,5070TG32.3% 150.3 / 12.7 65D492A (GAC→GCC) ybhFputative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (3/41);  new base G (21/0);  total (24/41)
Fisher's exact test for biased strand distribution p-value = 3.33e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACTGAATAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGA  >  NC_000913/827283‑827744
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
aCTGAATAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGGGACGCCgcacgac                                                                                                                                                                                                                                         >  2:342542/1‑233 (MQ=255)
aCTGAATAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGGGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGTGGGGGGGGACGCCCGAAgga                                                                                                                                                                                                                                         >  2:710039/1‑237 (MQ=255)
     aTAAAAGCTTGATCCGTGGTGGTTTCCGGTTGCTCATCGTTAGCCATCTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACGCCAGGCCGATGCGGTCGCAATATTCCGCGTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGtt                                                                                                                                                                                                                                     <  2:2118939/238‑1 (MQ=255)
           gCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtcgtc                                                                                                                                                                                                                                 <  1:1190373/236‑1 (MQ=255)
                     tggtggGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGGACGCCCGGAGGCGGTGTGTCGGGGAAAAGaaac                                                                                                                                                                                                                   >  2:1190373/1‑237 (MQ=255)
                        tggGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTc                                                                                                                                                                                                                 <  1:342542/238‑1 (MQ=255)
                          ggATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTAGGGTCGTCGGGAGAAAGAATgg                                                                                                                                                                                                                  >  1:917440/1‑234 (MQ=255)
                                            ttAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAATGTCCGGGTCATGGATCAGCGaaaa                                                                                                                                                                                            >  1:1811865/1‑237 (MQ=255)
                                             tAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTc                                                                                                                                                                                                                 <  1:1040855/217‑1 (MQ=255)
                                                 cGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgcca                                                                                                                                                                                       <  1:834980/239‑1 (MQ=255)
                                                    cTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGTGGGGGGGGACGCCCGGACGCGGGTCCGCGTGAAAACGAAAGTCCGGGTGATGGATCAGAGAACaaagccag                                                                                                                                                                                     >  2:2363959/1‑231 (MQ=255)
                                                           ttCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCg                                                                                                                                                                             <  1:1965380/239‑1 (MQ=255)
                                                                                         aaTTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGt                                                                                                                                             <  1:424701/239‑1 (MQ=255)
                                                                                                           aGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGa                                                                                                                               <  1:201124/235‑1 (MQ=255)
                                                                                                               cGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTGCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCAGGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACtttt                                                                                                                      <  1:1084364/238‑1 (MQ=255)
                                                                                                                  tGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGGAGTCGGGTCGTCGAGGAAAAGAAAGTCCGGTTCATGCATCAGCGAACAAGCCAGCGCCCGCCCCCGGTTAAAAACTAAAGGGAGGGTGTCGGGGGCGGGGGGGGGGGGACtttttcg                                                                                                                   >  2:1900316/1‑235 (MQ=255)
                                                                                                                  tGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCCACGCCGGCAGGTGGGTCGTCGAGAAAAAGAATGTCCGGGTCTATCATCAGCGAAAAAGGCAGCGCCAGCCGCCGTTTAAAAACTCATGGCAGTTCATCAGGGGCGTGGGGGGGGGTACTTTTcca                                                                                                                   >  2:679426/1‑236 (MQ=255)
                                                                                                                         gCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACg                                                                                                         <  1:679426/239‑1 (MQ=255)
                                                                                                                           aaTATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCAAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGc                                                                                                        <  1:574864/238‑1 (MQ=255)
                                                                                                                                     ttCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATg                                                                                              <  1:1491656/238‑1 (MQ=255)
                                                                                                                                       cATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGGCACGCCGGGAGGCGGGTTGTCGGGGAAAAAAAAGGGCCGGTCATGCATTCGCGAACAGGCCCACGCCCGCCGCCGTTTTAAAACTAATAGCAGTTCATCGGTGGCGGGGGGGGCGGGACATTTTCAGACCAAACGCCCCCTCATTCgg                                                                                          >  2:1820082/1‑239 (MQ=255)
                                                                                                                                          ccATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaa                                                                                                                                                                    <  1:699035/168‑1 (MQ=255)
                                                                                                                                                 gtgggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtcg                                                                                                                                                                                                                                   <  1:334486/100‑1 (MQ=255)
                                                                                                                                                    ggtggtGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCGGGAGTCGGGTCGTCGAGAAAAAGAAAGTCCGGGTTATTCATCAGCGAACAGGCCCCGGCCAGCCCCCGTTTTAAAACCAAttg                                                                                                                                                          >  2:699035/1‑166 (MQ=255)
                                                                                                                                                    ggtggtGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCCACGCCCGAAGTCGGGTCGTCGAGAAACAGAATGTCCGGGTCATGCATTAGCGAACAGGCCAGCGCCCGCCGCCGTTTTAAAACTAATGGCCGTTTATCGGCGGCGGGGGGGGGGGTACTTTTTCTGCCGAAACGCCCGCCCCTGCGGGGGAGTTTTTTGt                                                                             >  2:1351180/1‑239 (MQ=255)
                                                                                                                                                    ggtggtGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTTCTCGAGAAACAGAAAGGCAGGTTCATGCATCCCCGCACAGGCCCGCGCCAGCCCCCGTTTAGAAAACAAATGCCGGGTAACCGTGGCGTGGGGGGGGAGAACTTTTACGCCGCAAGCCCCGCTCCTACGGGGGAGTGTTTTGt                                                                             >  2:2067130/1‑239 (MQ=255)
                                                                                                                                                      tggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaa                                                                                                                                                                    <  1:1535285/156‑1 (MQ=255)
                                                                                                                                                           aCCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGaa                                                                                                                                                                                              <  1:840904/126‑1 (MQ=255)
                                                                                                                                                           aCCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGGAGGCGGGTCGTCGAGaaaaa                                                                                                                                                                                                                        >  2:334486/1‑98 (MQ=255)
                                                                                                                                                             cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATg                                                                                                                                                                                                        <  1:1659119/114‑1 (MQ=255)
                                                                                                                                                                     tGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGcc                                                                                                                                                                                                                                                >  1:432602/1‑67 (MQ=255)
                                                                                                                                                                     tGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGcc                                                                                                                                                                                                                                                <  2:432602/67‑1 (MQ=255)
                                                                                                                                                                     tGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGGGACGCCGGGAGGCGGGGCGTCGGGAAAAAGAAAGTCCGGTGCTTGCATCCGCGAAAAGggcaggg                                                                                                                                                                                    >  2:840904/1‑124 (MQ=255)
                                                                                                                                                                      gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGt                                                            <  1:1351180/238‑1 (MQ=255)
                                                                                                                                                                      gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGt                                                            <  1:707569/238‑1 (MQ=255)
                                                                                                                                                                      gACGCCTTTCTCTACCATGCTGTGGATGGGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGt                                                            <  1:901076/238‑1 (MQ=255)
                                                                                                                                                                       aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGTCGGGTTGTCGAGAAACAGAATGTCCGGTTCATTCCTCAGCGaa                                                                                                                                                                                              >  2:1659119/1‑114 (MQ=255)
                                                                                                                                                                        cGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAg                                                         <  1:1900316/239‑1 (MQ=255)
                                                                                                                                                                             ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTc                                                                                                                                                                                                                 <  1:538791/89‑1 (MQ=255)
                                                                                                                                                                             ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCa                                                                                                                                                                                                          <  2:1689256/96‑1 (MQ=255)
                                                                                                                                                                             ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATg                                                                                                                                                                                                        <  2:1111672/98‑1 (MQ=255)
                                                                                                                                                                             ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATa                                                    <  1:1016467/239‑1 (MQ=255)
                                                                                                                                                                             ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATa                                                    <  1:740528/239‑1 (MQ=255)
                                                                                                                                                                             ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATa                                                    <  1:2067130/239‑1 (MQ=255)
                                                                                                                                                                             ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGTGGGGGGCCACCCCCGGAGTCGGGGTCTCGAGAAACAGAATGGTCGGGTCCTGCATTAACGCAACGGCCCACGCCACCCCCCGTTTAAAAACTACAGGCCGGTTATCGGTGGGGGGGGGGGGGGGACTATTTTAGGCGCAAGCCCCCCCCCATGGGGGGGGTTTTTTTGTCTGCGCCCGACCCCGCCAGcccca                                                    >  2:1969818/1‑236 (MQ=255)
                                                                                                                                                                               ctctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGaa                                                                                                           <  1:2013338/183‑1 (MQ=255)
                                                                                                                                                                                  tACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtcgtcg                                                                                                                                                                                                                               <  1:2214034/70‑1 (MQ=255)
                                                                                                                                                                                  tACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGgaga                                                                                      <  1:1099830/200‑1 (MQ=255)
                                                                                                                                                                                       tGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGAGAAACAGAAAGTGCGGGTCATGc                                                                                                                                                                                                       >  2:538791/1‑89 (MQ=255)
                                                                                                                                                                                       tGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGTCGGTTCGTCGGGAAAAAGAATGGCCGGGTTATGCATCAGCgca                                                                                                                                                                                              >  1:1111672/1‑96 (MQ=255)
                                                                                                                                                                                       tGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGGAGGCGGTTTGGCGAGAAAAAGAATGTCCGGTTTATGCATCAgtg                                                                                                                                                                                                >  1:1689256/1‑94 (MQ=255)
                                                                                                                                                                                        gCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAgaga                                         <  1:1227651/239‑1 (MQ=255)
                                                                                                                                                                                         cTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGTCAACGCCCGAAGTCGGGTTGTCGAGAAAAAGAATGTCCGGTTCATGCATCAGCGGAAAGGCCAGCGCCAGCCGCTGTTTAAAACCCAATTGCAGGTTCTTCGTGGGGGGGGGGGGGGTACTTTTTCCGCCGCACGCCCCGccc                                                                                                 >  2:2013338/1‑181 (MQ=255)
                                                                                                                                                                                         cTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTc                                                                                                                                                                                                                 >  1:2175040/1‑77 (MQ=255)
                                                                                                                                                                                         cTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTc                                                                                                                                                                                                                 <  2:2175040/77‑1 (MQ=255)
                                                                                                                                                                                            ttGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAATCCTAATGGCAGTTCGTCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGtt                                                                            <  1:2282048/200‑1 (MQ=255)
                                                                                                                                                                                            ttGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAgg                                                                                                                                  <  1:494860/151‑1 (MQ=255)
                                                                                                                                                                                            ttGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGgaga                                                                                      <  1:286696/190‑1 (MQ=255)
                                                                                                                                                                                            ttGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAATCGGGTTTTCGAGAAAAAGaa                                                                                                                                                                                                                     >  2:2214034/1‑70 (MQ=255)
                                                                                                                                                                                            ttGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAAAAAAATGTCCGGGTCCTGCATCAGCGGACAGGCCAGCGCCCGCCGCCGGTTTAAAACTAATGGGAGGTCCTTGGGGGGGGGGGGGGGGATACTTTTTCGGCCGGAACGCCCGCCCCTGCCGGGGGATTTTTTGTTTTTCGCCCCCACCCGCAACCCCAACACACCACAGaaaaa                                     >  2:2055537/1‑239 (MQ=255)
                                                                                                                                                                                                 gtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGtt                                                                            <  1:413459/195‑1 (MQ=255)
                                                                                                                                                                                                          aaaTTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGtt                       <  1:2112408/239‑1 (MQ=255)
                                                                                                                                                                                                                     gggTGGTGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCGGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGt                <  1:935781/235‑1 (MQ=255)
                                                                                                                                                                                                                     gggTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGt                <  1:362816/235‑1 (MQ=255)
                                                                                                                                                                                                                     gggTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGt                <  1:1037406/235‑1 (MQ=255)
                                                                                                                                                                                                                                tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGAt                                                                                                                              <  2:526451/119‑1 (MQ=255)
                                                                                                                                                                                                                                tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGAt                                                                                                                              >  1:526451/1‑119 (MQ=255)
                                                                                                                                                                                                                                tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTaga  >  2:632590/1‑238 (MQ=255)
                                                                                                                                                                                                                                tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCAGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTaga  <  1:946457/238‑1 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
ACTGAATAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGA  >  NC_000913/827283‑827744

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