Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913595,1650CT70.9% 62.1 / 36.3 55L93L (TTG→TTA) cusRDNA‑binding transcriptional activator CusR
Reads supporting (aligned to +/- strand):  ref base C (8/8);  new base T (17/22);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTT  >  NC_000913/594920‑595411
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
aTAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTc                                                                                                                                                                                                                                                   <  1:451082/251‑1 (MQ=255)
          ggggCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTAAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCAAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTACCCCCAACTCCAGCCCCTTg                                                                                                                                                                                                                                         <  1:162030/251‑1 (MQ=255)
               agcaCTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGAcgcg                                                                                                                                                                                                                                    >  1:223753/1‑251 (MQ=255)
                     tCGCTCTGAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACGTTGCGGCTGACGAGATCGACCATCAAACCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCCCCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGtt                                                                                                                                                                                                                               <  1:121981/250‑1 (MQ=255)
                     tCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTc                                                                                                                                                                                                                              >  2:129427/1‑251 (MQ=255)
                          cTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATgg                                                                                                                                                                                                                         >  1:30399/1‑251 (MQ=255)
                               ggCGAAGGAAGAACTCCAGCAGAGTAAACTCTGTGCTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCGCCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGcc                                                                                                                                                                                                                     <  2:30399/250‑1 (MQ=255)
                                  gaaggaagAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAgc                                                                                                                                                                                                                 <  1:3573/251‑1 (MQ=255)
                                                agcagAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAAc                                                                                                                                                                                                     >  1:552886/1‑249 (MQ=255)
                                                              tttACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATAACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCATTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCtt                                                                                                                                                                                     >  1:317792/1‑251 (MQ=255)
                                                                 aCTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTa                                                                                                                                                                                   <  1:473092/250‑1 (MQ=255)
                                                                    ggtaaaaGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAAGGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTg                                                                                                                                                                                <  2:575166/246‑1 (MQ=255)
                                                                    ggTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCGGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAAGGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGTTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTgg                                                                                                                                                                               <  2:626114/251‑1 (MQ=255)
                                                                    ggTCAAAGGGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTAtt                                                                                                                                                                                 <  2:223512/249‑1 (MQ=255)
                                                                     gTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGCCCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGc                                                                                                                                                                              >  1:369356/1‑251 (MQ=255)
                                                                     gTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGc                                                                                                                                                                              <  2:447161/251‑1 (MQ=255)
                                                                          aGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGCCGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGCTAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGc                                                                                                                                                                         >  2:97361/1‑251 (MQ=255)
                                                                                    gTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAACCTGACTTTCGCCAATCACCGCCGCCCCGCGACGCAGTAATGTGCGCACCCGCGCCCGCAGTTCAGCAAAAGCGTATGGCTTCACCAGGCAGCCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGGGCCAAGCGCAGTACGCACCAGAATCGGCATCCCTTTaagatcggacgatcgtagtgg                                                                                                                                                               >  2:103110/1‑231 (MQ=255)
                                                                                    gTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTa                                                                                                                                                                                   <  1:103110/231‑1 (MQ=255)
                                                                                      gCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCACGCGCCGCAGTCATGTGCGAACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTGACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCTTAAAATgcgc                                                                                                                                                              >  2:233873/1‑250 (MQ=255)
                                                                                        cGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACGTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACg                                                                                                                                                           <  1:381753/251‑1 (MQ=255)
                                                                                               ggTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGCTGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAACCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATccc                                                                                                                                                    >  1:451303/1‑251 (MQ=255)
                                                                                                  gACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGc                                                                                                                                                 <  1:308503/251‑1 (MQ=255)
                                                                                                          ggCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACg                                                                                                                                         <  1:382295/251‑1 (MQ=255)
                                                                                                                             cAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCg                                                                                                                       >  1:201800/1‑250 (MQ=255)
                                                                                                                             cAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCg                                                                                                                       <  1:571582/250‑1 (MQ=255)
                                                                                                                             cAAATCGGCAACCGGAAACTGACTTTCGATAATCGCCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGa                                                                                                                      <  2:375279/251‑1 (MQ=255)
                                                                                                                                 tCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGAt                                                                                                                  >  1:88939/1‑251 (MQ=255)
                                                                                                                                  cGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGAtt                                                                                                                 >  2:536622/1‑251 (MQ=255)
                                                                                                                                           gAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGAGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGTGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATc                                                                                                         >  1:347793/1‑250 (MQ=255)
                                                                                                                                                  aCTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCAc                                                                                                 >  2:137056/1‑251 (MQ=255)
                                                                                                                                                     tcctaTAAACACCGCCGCCCCGCGCCGCAGTAATGGGCGCACCCTCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTATTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATGGGCGGAGCGGAACTTGCTCACGATATCCTAGCCTTTCACGTCCGGCAGCATAATATCTATGATTAGCAGATCAGAATctcc                                                                                                <  2:372500/245‑4 (MQ=255)
                                                                                                                                                       cGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGt                                                                                             <  1:26541/250‑1 (MQ=255)
                                                                                                                                                       cGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTc                                                                                            >  1:374758/1‑251 (MQ=255)
                                                                                                                                                            aTCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGc                                                                                       >  1:52391/1‑251 (MQ=255)
                                                                                                                                                                     gccCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCACCAGATGGTa                                                                              <  1:145451/251‑1 (MQ=255)
                                                                                                                                                                           cgccgcAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAAGGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCAtt                                                                        <  2:201800/251‑1 (MQ=255)
                                                                                                                                                                                cagtatggtacgcACACGCGCATGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGATCCTTGACGCGATGGTCAATGGTGTCTAGCGCGGTAAGCAACAGAATCGGCTTCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGCATTATCAGATCATAATCGCCGTTCATCGCCAGCTGGTAGCCATTAAGct                                                                   <  2:161397/241‑2 (MQ=255)
                                                                                                                                                                                 aGTAATGTGCGCACCCTCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCACCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGccc                                                                  <  2:494949/251‑1 (MQ=255)
                                                                                                                                                                                    aaTGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGtt                                                               <  2:374758/251‑1 (MQ=255)
                                                                                                                                                                                    aaTGTCCGCACCCGCGCCAGCTGTTCCGCATAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCCGCCCCTTGCCGCTATGTTCAATGGTGCCAAGCGCGGTAAGCTACCGAATCGGCATCCCTTTATTTGCGGTGCGTACCATTCGCTCGATATCCCTGCCGTTCACGTCCGGCAGCATAATATCGATGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGccc                                                                  <  1:422144/248‑1 (MQ=255)
                                                                                                                                                                                      tgtgCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGt                                                             >  1:113979/1‑251 (MQ=255)
                                                                                                                                                                                       gtgCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTc                                                            <  1:457420/251‑1 (MQ=255)
                                                                                                                                                                                              cccGCGCCAGCAGTTCATCAAAAGAGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCaaa                                                     <  2:314372/251‑1 (MQ=255)
                                                                                                                                                                                               ccGCGCCAGCATTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGAGCATAATCACCGGTCATCGCCAGAGGGTAGCCATTCAGCCCGTTGTCGGCCAAAt                                                    <  2:463293/251‑1 (MQ=255)
                                                                                                                                                                                                  cgcCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCg                                                  >  1:496819/1‑250 (MQ=255)
                                                                                                                                                                                                             cAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCg                                      <  1:137056/251‑1 (MQ=255)
                                                                                                                                                                                                                 aaaaGCGAATGGCTTCACCAGGAAGTCATCTGCCCCCAACTCCAGCCCCTTGCCGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATAACGAGGATTATCAGATCATCATCACCGGTCATCGCCCGATGGTAGCCATTCAGCCCGTTGTCGGCCCACTCGACCACAAAACCGGCtt                                  >  2:243739/1‑251 (MQ=255)
                                                                                                                                                                                                                       gAATGGCTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTa                            <  1:97361/251‑1 (MQ=255)
                                                                                                                                                                                                                        aaTGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTaa                           <  1:73629/251‑1 (MQ=255)
                                                                                                                                                                                                                             cTTCACCAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAACTCGACCACAAAACCGGCTTCGGTTAACCCtt                      >  2:378010/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ccAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCTATGTTCAAAGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAAACGGCTTCGGTTAACCCTTTGGTc                 >  2:59496/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cAGGTAGTCATCTGCCCCTAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGCTTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTc                 >  1:610378/1‑250 (MQ=255)
                                                                                                                                                                                                                                      gTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGt             >  1:589778/1‑251 (MQ=255)
                                                                                                                                                                                                                                       tAGTCATCTGCCCCTACCTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGt             >  2:428774/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 cccccAACTCCAGCCCCTTGACGCTATGTTCAATGGTTCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGATCGTAACATGCGCACGATATCCCATCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGtt  <  1:165723/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
ATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTT  >  NC_000913/594920‑595411

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