Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,333,6630GA76.9% 60.1 / 16.9 40D96D (GAC→GAT) basRDNA‑binding transcriptional dual regulator BasR
Reads supporting (aligned to +/- strand):  ref base G (5/4);  new base A (15/15);  total (21/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 = 3.15e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCGGCAGGGTATATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGT  >  NC_000913/4333416‑4333906
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
aGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATATTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGtcatcg                                                                                                                                                                                                                                                  <  1:98433/251‑1 (MQ=255)
    gTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCAc                                                                                                                                                                                                                                              <  1:412935/251‑1 (MQ=255)
                 aTTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATATTCAGCGTCAGATTGCCAACAAGCAGCTCACTTTCGCCCTGATTATTATGGCGTCGGAGCAGGGCGCGGAGACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGccc                                                                                                                                                                                                                                 <  1:125318/251‑1 (MQ=255)
                        gATGCACAGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATATTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGTGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCCGATAGCAATCGGCACCGACATCCAGCCCGGCgctt                                                                                                                                                                                                                          >  2:106759/1‑248 (MQ=255)
                                  aCTGCCTGCTTTGAGCATTAACCGTGACAGCAGCGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATATTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTc                                                                                                                                                                                                                 >  1:386227/1‑250 (MQ=255)
                                    tgcctgcTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATATTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCa                                                                                                                                                                                                                                               <  1:125215/218‑1 (MQ=255)
                                    tgcctgcTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATATTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCa                                                                                                                                                                                                                                               >  2:125215/1‑218 (MQ=255)
                                         gcTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATATTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATc                                                                                                                                                                                                         <  1:64143/251‑1 (MQ=255)
                                         gcTTTGAGCATTAACAGAGCCCGCAGAGCCTCTTCTTTGGACGTCAGAATCAACTCTTCAACGCCCATCCATAACTGACGGATACCCATCTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTAGCCCTGATTATTATAGCGTCGTAGCAGAGCGAGGATACTGGCGTATAACTCTTCCAGCGCAACAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCACACCGCTTTTGTCGGTCAGCGTATc                                                                                                                                                                                                         >  2:22250/1‑251 (MQ=255)
                                                  aTTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATATTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGt                                                                                                                                                                                                <  1:177449/251‑1 (MQ=255)
                                                   ttAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATATTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTg                                                                                                                                                                                               >  2:105421/1‑251 (MQ=255)
                                                           gACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCa                                                                                                                                                                                        <  1:97328/250‑1 (MQ=255)
                                                             cagcagAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTTACGGCGACCCATATTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGAtt                                                                                                                                                                                                                          <  2:353547/214‑1 (MQ=255)
                                                             cagcagAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATATTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGAtt                                                                                                                                                                                                                          >  1:353547/1‑214 (MQ=255)
                                                               gcagagcagaTTATTTGGTCGTCATACTCAACTCTCCACCGCCCAGCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAGCCGCTCACTTTCGCCCTGATTATTATGGCGTCTTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCCCCAGATAGTCGTCGGCACCGACATCCAGCCCTGCGATTTTGTCGGTCAGCGTATTGCGAGCTGTGAGGATCAGTAcc                                                                                                                                                                                   <  2:305724/251‑1 (MQ=255)
                                                                              tGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATATTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACTGGCAGGGTGTAtttt                                                                                                                                                                    >  1:160535/1‑251 (MQ=255)
                                                                                     cAGAATCTACTCTTCACCGCCCATCCATACCTGACGGCGACCCATATTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGAACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACTGGCAGGGTGTATTTTTTCTGcg                                                                                                                                                             >  1:366696/1‑250 (MQ=255)
                                                                                               tctTCACCGCCCATCCATACCTGACGGCGACCCATATTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTg                                                                                                                                                                                                                       <  2:153144/183‑1 (MQ=255)
                                                                                               tctTCACCGCCCATCCATACCTGACGGCGACCCATATTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTg                                                                                                                                                                                                                       >  1:153144/1‑183 (MQ=255)
                                                                                                  tCACCGCCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGGGAGGATCAGTACCGGCAGGGTATATTTTTTCTGCCGGATACGGGc                                                                                                                                                   <  1:31119/248‑1 (MQ=255)
                                                                                                           atccatACCTGACGGCGACCCATGTTCAGCGTCAGATTACCAACCATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCAGATACGGGCATGTACCTCTTCCAGCGCAAACGGCTTCACCCGATAGTCGTCGGCACCGACATCCAGCCCGGAGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGCTCAGTACCGGCAGGGTATATTTTTTCTGCCGGATACGGGCGAGAAAATGc                                                                                                                                         >  1:253672/1‑249 (MQ=255)
                                                                                                                 aCCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCGGCAGGGTATATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATc                                                                                                                                 >  2:310700/1‑251 (MQ=255)
                                                                                                                      aCGGCGACCCATATTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCCCCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACTGGCAGGGTGTATTTGTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTtcgt                                                                                                                            <  1:182156/251‑1 (MQ=255)
                                                                                                                      aCGGCGACCCATATTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACTGGCAGGGTGTATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTtcgt                                                                                                                            >  1:304337/1‑251 (MQ=255)
                                                                                                                                  gTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCGGCAGGGTATATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCt                                                                                                                 >  1:7067/1‑250 (MQ=255)
                                                                                                                                    tCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACTGGCAGGGTGTATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTaa                                                                                                               >  1:271513/1‑250 (MQ=255)
                                                                                                                                    tCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGAGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCCAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCCGTACTGGCAGGGTGTATTTTTTCTGCCGGTTACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTaa                                                                                                               >  2:178162/1‑250 (MQ=255)
                                                                                                                                                      caaTCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACTGGCAGGGTGTATTTTTTCTGCCGGATACGGGCGGGAAAATGCAGTCCATCTTCGTCGGGGACCCCTAAATCCAGTaccacccgg                                                                                               >  2:237116/1‑245 (MQ=255)
                                                                                                                                                      caaTCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACTGGCAGGGTGTATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAgg                                                                                               <  1:237116/248‑1 (MQ=255)
                                                                                                                                                                        ccTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACTGGCAGGGTGTATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCGGCCTCAAgg                                                                          >  1:71558/1‑251 (MQ=255)
                                                                                                                                                                        ccTGATTAATAAGGCGAAGTAGCAACGCACGGATATGGGCATGTAACACTTAACGCGCAACAGAATTCACCAACTAGTCGTCGGCACCGCCATACAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGCAAGGATAAGTACCGGCAGGGTATATTTTTTCTGCAGGATACGGCCGAGGACATGCAGTCAATATTCGTCGGGTAAACCTAAAACCAGTACAACCAGGCTGCAATGAACTGCATCaacg                                                                          >  2:207385/1‑249 (MQ=255)
                                                                                                                                                                           gATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCGGCAGGGCATATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGGAACCCTAAATCCAGTACCACCAGGCTGTAATGCCCTGCCTCAAGGCtt                                                                       >  2:49501/1‑251 (MQ=255)
                                                                                                                                                                                 tatGGCGTCGTAGCAGGGCGCGGATACGGGCTTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACTGGCAGGGTGTATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAGCCCTAAATCCAGTACCACCAGGCTGTAATGACCGGCCTCAAGGCTTTGTTc                                                                  <  2:160535/250‑1 (MQ=255)
                                                                                                                                                                                           tAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCCTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACTGGCCGGGTGTATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTtcg                                                                                                                             >  2:193073/1‑181 (MQ=255)
                                                                                                                                                                                           tAGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACTGGCAGGGTGTATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTtcg                                                                                                                             <  1:193073/181‑1 (MQ=255)
                                                                                                                                                                                            aGCAGGGCGCGGATACGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACTGGCAGGGTGTATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCcgcg                                                       <  2:69720/250‑1 (MQ=255)
                                                                                                                                                                                                          acgggcgcgTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGGCAGCGGCACCGACATCCAGCCCGGCGATTTTGGCGGTCATCGTAGCGCGAGCGGTGAGGATCAGTACTGGCATGGTGTATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGTGTATCCCTAAAGCCAGTACCACCAGGCGGGAATGACCGGCCTCAAGGCTTTTTTCCGCCATCCGCGCGGTTGTCACGccat                                        <  2:386227/243‑4 (MQ=255)
                                                                                                                                                                                                          aCGGGCGTGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACTGGCAGGGTGTATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCGGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGccat                                        >  1:290181/1‑248 (MQ=255)
                                                                                                                                                                                                                                       gCTTCACCAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACTGGCAGGGTGTATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCGGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCCATCGCACGCGTAGCCTTCGGTTTGCgccgcc           <  1:105421/251‑1 (MQ=255)
                                                                                                                                                                                                                                             ccAGATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACTGGCAGGGTGTATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCGGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCCATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATc     <  2:71558/251‑1 (MQ=255)
                                                                                                                                                                                                                                                gATAGTCATCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACTGGCAGGGTGTATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAGCCCTAAATCCAGTACCACCAGGCTGTAATGACCGGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCCATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGGATCaat  <  1:114287/251‑3 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
AGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCGGCAGGGTATATTTTTTCTGCCGGATACGGGCGAGAAAATGCAGTCCATCTTCGTCGGGTAACCCTAAATCCAGTACCACCAGGCTGTAATGACCTGCCTCAAGGCTTTGTTCCGCCATCCGCGCGGTTGTCACGCTATCGCACGCGTAGCCTTCGGTTTGCGCCGCCAGAATCAGT  >  NC_000913/4333416‑4333906

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