Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,334,4060AC72.7% 64.4 / 27.1 44P395P (CCT→CCGeptAphosphoethanolamine transferase EptA
Reads supporting (aligned to +/- strand):  ref base A (7/5);  new base C (18/14);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATCCGGGGCGATGGCATAAGGCAGACCGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGGTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTTGCAGGTTATTGATGTACTCTTCAAGCCCGTGGAACAGCACTTCGTCATAGCATTCGCCGTTGATGCACTGATCAGGTAGATTCAGCGCGGTGACGTTCTGGTGAGGCACGCGGTCGCAGGCACCTTTACAGCCGCCATCGTTGTCATTCCACAGCACGTTGATGCCCGCTCGCTGAAT  >  NC_000913/4334160‑4334654
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        
aTCCGGGGCGATGGCATAAGGCAGACCGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGGTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGAGGCg                                                                                                                                                                                                                                                      >  2:185364/1‑251 (MQ=255)
     gggCGATGGCATAAGGCAGACCGTGCCGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAACGATCACCCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTACACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAAATGCGGCGGAt                                                                                                                                                                                                                                                   >  2:185524/1‑249 (MQ=255)
           tGGCATAAGGCAGAACGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACAAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCCTTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCgggtg                                                                                                                                                                                                                                           >  2:142924/1‑248 (MQ=255)
                tAAGGCAGACCGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAAt                                                                                                                                                                                                                                       >  1:394712/1‑250 (MQ=255)
                 aaGGCAGACCGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAAt                                                                                                                                                                                                                                       >  1:136593/1‑249 (MQ=255)
                  aGGCAGACCGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACATGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAgg                                                                                                                                                                                                                                    <  2:394712/251‑1 (MQ=255)
                                gatagatGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGATGGTAAATTTATCCTGATGTTATTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGGCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGTTCGCGGTTGTAATAGGTAGGACCGAGGcgg                                                                                                                                                                                                                      >  2:361002/1‑249 (MQ=255)
                                    gatGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGGTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCg                                                                                                                                                                                                                   >  1:298131/1‑250 (MQ=255)
                                             ttCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGgtgtgt                                                                                                                                                                                                          <  2:136593/250‑1 (MQ=255)
                                                aCCTAACGATTCACCCTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTATTTCAGCAGATTACTCGCTTTATCAAAAATATAGTCGCCGTAAACCAGAGTGTTGTCGTAAGTGTTCAACAGATGCTCTTTGGACCAGGCCTGGAGCTCATTGGTGTCGCAGGATTGGGTAAATTTACTGAACTGAGGCGTATATCGGTTAGACTAGGTCGGACCGTGTCTGCAGATGGTGGGTacg                                                                                                                                                                                                       >  2:143607/1‑248 (MQ=255)
                                                       gATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCCGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCAc                                                                                                                                                                                               >  1:388712/1‑251 (MQ=255)
                                                                 ggTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCAACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCAccc                                                                                                                                                                                     >  1:404942/1‑250 (MQ=255)
                                                                     aGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGTCGGATAGCGGTTGTAATAGTTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCa                                                                                                                                                                                 <  1:179506/251‑1 (MQ=255)
                                                                                ccAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATg                                                                                                                                                                      <  1:371896/251‑1 (MQ=255)
                                                                                  aGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGt                                                                                                                                                                     >  1:357041/1‑250 (MQ=255)
                                                                                      tggtggTAAATTTAACCTGATGTTCTTTCAGCAGATTAATCGCTTTATCACCAATATAGTCGACGTAAACCCGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTActct                                                                                                                                                                >  1:304836/1‑251 (MQ=255)
                                                                                          ggtAAATTTATCCGGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGGCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTcc                                                                                                                                                             <  2:335441/250‑2 (MQ=255)
                                                                                             aaaTTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGcc                                                                                                                                                         >  1:300881/1‑251 (MQ=255)
                                                                                              aaTTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGccc                                                                                                                                                        <  2:357041/251‑1 (MQ=255)
                                                                                                 ttATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCt                                                                                                                                                                                                                       <  1:17956/185‑1 (MQ=255)
                                                                                                 ttATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTATTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTACATTTCCTGAACTGCGGCGGATTGCGGTTGTAATAGGTCGGACCGTGGCt                                                                                                                                                                                                                       >  2:17956/1‑185 (MQ=255)
                                                                                                  tattctgatgttagtgtatcagATTGATCGCTTTATCAACAATATAGCAGACGTGAACCAGCTTGTTTTCGAAAGTTTTCACCTGTTGCTCTTTGCTACAGTTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGta                                                                                                                                                     <  2:14387/231‑2 (MQ=255)
                                                                                                                      agATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGATATTGATGTACTCTTCCAGCCCGTGGCACAGTACTTCGTCATAGc                                                                                                                                 >  1:93431/1‑250 (MQ=255)
                                                                                                                      agATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTGGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCTGACCGTGGCTGCCGATGGTTTGTAAGACAATCACGCAATCACCCTGGAGGTTATTGATGTACTCTTCCATCCCGTGGAACATTACTTCGTCATAGCa                                                                                                                                <  2:162046/251‑1 (MQ=255)
                                                                                                                       gATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGGTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTTGCAGGTTATTGATGTACTCTTCAAGCCCGTGGAACAGCACTTCGTCATAGCAt                                                                                                                               <  2:298131/251‑1 (MQ=255)
                                                                                                                          tAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTACGTGTTCACCAGTTGCTCTTTGGTACCGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTTGCCGGTTATTGCTGTACTCTTAAAGCCCGTGGAACAGCACTTCGTCATAGCATTct                                                                                                                            >  2:208256/1‑250 (MQ=255)
                                                                                                                           aaTCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGTGGAACAGTACTTCGTCATAGCATTCGc                                                                                                                           >  1:146699/1‑251 (MQ=255)
                                                                                                                            aTCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGTGGAACAGTACTTCGTCATAGCATTCGcc                                                                                                                          <  1:125520/251‑1 (MQ=255)
                                                                                                                                    aTCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGGTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTTGCAGGTTATTGATGTACTCTTCAAGCCCGGGGAACAGCACTTCGTCATAGCATTCGCCGTTGATg                                                                                                                   <  1:388564/250‑1 (MQ=255)
                                                                                                                                                   gACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGGTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTTGCAGGTTATTGATGTACTCTTCAAGCCCGTGGAACAGCACTTCGTCATAGCATTCGCCGTTGATGCACTGATCAGGt                                                                                                       >  1:192709/1‑247 (MQ=255)
                                                                                                                                                      gTAAAACAGCGTGCTGTCGTAAGCGTTCACCAGTTCCTCTTTGGTACAGGTCTGGATCTCATTGGTGTCGCCGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGCACTAGGTCGGACCGCGGCTGCCGCTGGTGCGTAAGACAATCACGCCATCACCTTGCAGGTTATTGCTGTACTCTTCAAGCCCGTGGAACAGCACTTCGTCATAGCATTCGCCGTTGATGCACTGATCCGGCAGATTc                                                                                                 >  2:382466/1‑250 (MQ=255)
                                                                                                                                                         aaCCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGATACAGGTCTGGATCTAATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGTGGAACAGTACTTCGTCATAGCATTCGCCGTTGATGCACTGACCAGGCAGATTCAgcg                                                                                             >  1:182128/1‑251 (MQ=255)
                                                                                                                                                          aCCAGCGCGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGCTCTCATTGGTGTCGCAGGTTGGGGTAACTTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGTGGAACAGTACTTCGTCCTAGCATTCGCCGTTGCTGCCCTGACCAGGCAGATTCAgcgc                                                                                            >  2:201703/1‑251 (MQ=255)
                                                                                                                                                                       cGTAAGTGTTCACCAGTTGCTCTTTGGTACAGGTCTGGATCTCATTGGTGTTGCAGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGCCCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTGGCATGTTATTGATGTACTCTTCAAGCCCGTGGAACAGCACTTCGTCATAGCATTCGCCGTTGATGCACTGATCAGGTAGATTCAGCGCGGTGACGTTCTgg                                                                               <  2:405631/251‑1 (MQ=255)
                                                                                                                                                                       cGTAAGTGTTCACCAGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGTGGAACAGTACTTCGTCATAGCATTCGCCGTTGATGCACTGACCAGGCAGATTCAGCGCGGTGACGTTCTgg                                                                               >  2:88803/1‑251 (MQ=255)
                                                                                                                                                                                     aGTTGCTCTTTGCTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGGGGAACAGTACTTCGTCATAGCATTCGCCGTTGATGCACTGACCAGGCAGATTCAGCGCGGTGACGTTCTGGTGAGGCACGCGGTc                                                                 <  1:185524/251‑1 (MQ=255)
                                                                                                                                                                                                gctACAGGTCGGGAGCTAATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATTGCGGTTGTAATAGGTCGGACCGTGGTTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGTGGAACAGTACTTCGTCATAGCATTCGCCGTTGATGCACTGACCAGGCAGATGCAGCGCGGTGACGTTCTGGTGAGGCACGCGGTCGCAGGCACCt                                                       <  1:50954/248‑1 (MQ=255)
                                                                                                                                                                                                 gTACAGGTCTGGATCTCATTGGTGTCGCAGGCTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTTGCAGGTTATTGATGTACTCTTCAAGCCCGTGGAACAGCACTTCGTCATAGCATTCGCCGTTGATGCACTGCTCAGGTAGATTCAGCGCGGTGACGTTCTGGTGAGGCACGCGGTCGCAGGCACCttt                                                     >  2:208251/1‑251 (MQ=255)
                                                                                                                                                                                                   aCAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTTGCAGGTTATTGATGTACTCTTCAAGCCCGTGGAACAGCACTTCGTCATAGCATTCGCCGTTGATGCACTGATCAGGTAGATTCAGCGCGGTGACGTTCTGGTGAGGCACGCGGTCGCAGGCACCTTTAc                                                   <  1:185364/251‑1 (MQ=255)
                                                                                                                                                                                                         cTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGTGGAACAGTACTTCGTCATAGCATTCGCCGTTGATGCACTGACCAGGCAGATTCAGCGCGGTGACGTTCTGGTGAGGCACGCGGTCGCAGGCACCTTTACAgccgc                                             >  1:199668/1‑251 (MQ=255)
                                                                                                                                                                                                                gcATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGTGGAACAGTACTTCGTCATAGCATTCGCCGTTGATGCACTGACCAGGCAGATTCAGCGCGGTGACGTTCTGGTGAGGCACGCGGTCGCAGGCACCTTTACAGCCGCCATCGtt                                      <  1:361002/250‑1 (MQ=255)
                                                                                                                                                                                                                             aGGTTGGGGTAAATTTCCTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGTGGAACAGTACTTCGTCATAGCATTCGCCGTTGATGCACTGACCAGGCAGATTCAGCGCGGTGACGTTCTGGTGAGGCACGCGGTCGCAGGCACCTTTACAGCCGCCATCGTTGTCATTCCACAGc                         <  1:201703/251‑1 (MQ=255)
                                                                                                                                                                                                                                             ccTGAACTGCGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCCCGCCATCACCCTGCAGGTTATTGATGTACTCTTCCAGCCCGTGGAACCGTACTTCGTCATAGCATTCGCCGTTGATGCACTGACCCGGCAGATTCAGCGAGGTGACGTTCTGGTGAGGCACGCGGTCGCAGGCACCTTTACAGCCGCCATCGTTGTCCTTCCACAGCACGTTGATGCCcgcac         >  2:410428/1‑249 (MQ=255)
                                                                                                                                                                                                                                                    ggAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTTGCAGGTTATTGATGTACTCTTCAAGCCCGTGGAACAGCACTTCGTCATAGCATTCGCCGTTGATGCACTGATCAGGTAGATTCAGCGCGGTGACGTTCTGGTGAGGCACGCGGTCGCAGGCACCTTTACAGCCGCCATCGTTGTCAGTCCACAGCACGTTGATGCCCGCTCGCTGAAt  <  1:208256/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        
ATCCGGGGCGATGGCATAAGGCAGACCGTGCAGATAGATGCCATTTTCACCTAACGATTCACCGTGGTCAGAAAGATAAACCAGGCTGGTGGTAAATTTATCCTGATGTTCTTTCAGCAGATTAATCGCTTTATCAACAATATAGTCGACGTAAACCAGCGTGTTGTCGTAAGTGTTCACCAGTTGCTCTTTGGTACAGGTCTGGATCTCATTGGTGTCGCAGGTTGGGGTAAATTTCCTGAACTGAGGCGGATAGCGGTTGTAATAGGTCGGACCGTGGCTGCCGATGGTGTGTAAGACAATCACGCCATCACCTTGCAGGTTATTGATGTACTCTTCAAGCCCGTGGAACAGCACTTCGTCATAGCATTCGCCGTTGATGCACTGATCAGGTAGATTCAGCGCGGTGACGTTCTGGTGAGGCACGCGGTCGCAGGCACCTTTACAGCCGCCATCGTTGTCATTCCACAGCACGTTGATGCCCGCTCGCTGAAT  >  NC_000913/4334160‑4334654

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