Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT38.2% 43.1 / 12.0 35S62F (TCT→TTT) bglJDNA‑binding transcriptional regulator BglJ
Reads supporting (aligned to +/- strand):  ref base C (2/19);  new base T (13/0);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 1.13e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGT  >  NC_000913/4604105‑4604592
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                            
gcgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCa                                                                                                                                                                                                                                                   <  2:377877/251‑1 (MQ=255)
                  aacttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGGgcgccggg                                                                                                                                                                                                                                 >  2:224304/1‑247 (MQ=255)
                          ttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGTCATGAGAAGAGAgcgccgggggggtatt                                                                                                                                                                                                                         >  2:372171/1‑239 (MQ=255)
                           tgcattCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTCTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTa                                                                                                                                                                                                                         <  1:7637/246‑1 (MQ=255)
                                tCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTg                                                                                                                                                                                                                    <  1:18890/250‑1 (MQ=255)
                                                   cAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGAGGAGGGAGCGCAGAGAGGGGTTATCTTGCctgactgtacaggtggtt                                                                                                                                                                                                 >  1:332905/1‑239 (MQ=255)
                                                        tGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGAGAGCGCAGGGAGGGGTTTTCTTGCCTGCCTTAACTGGCGAtttaggtt                                                                                                                                                                                           >  2:76225/1‑245 (MQ=255)
                                                             cacaGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGTCGATTAAGTTTCcgcg                                                                                                                                                                                    <  2:123169/251‑1 (MQ=255)
                                                                                       cGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatg                                                                                                                                                          <  2:193082/251‑1 (MQ=255)
                                                                                         cACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcgc                                                                                                                                                                                                                                     <  1:192489/176‑1 (MQ=255)
                                                                                                aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTg                                                                                                                                                 <  1:62006/251‑1 (MQ=255)
                                                                                                   aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGCGAGGGGTTATTTTTCCTGACTGAACTGGCGGTTAAGTTTTCGCGTACCCCGCGTTTATTTTTTGCGGGTTATGATATTGa                                                                                                                                                >  2:10654/1‑249 (MQ=255)
                                                                                                     aaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGGGAGCGCAGAGAGGGATTATTTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCt                                                                                                                                            >  1:201172/1‑251 (MQ=255)
                                                                                                          gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGTGAAGTGAGCGCAGAGGGGGGTTATCTTTCCTGACTGTACTGGCGATTTAGTTTTCGCGGTCCCGGCGGTTTGTTATTTCTGGTGATGGTATTTTAGCTCGGcc                                                                                                                                       >  2:1629/1‑250 (MQ=255)
                                                                                                             aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTTCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTTGTTATTTGGGGTGATGATATTGAAGCTCGGCTGAt                                                                                                                                   >  1:4599/1‑251 (MQ=255)
                                                                                                                agtagtATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTgg                                                                                                                               <  1:315650/251‑1 (MQ=255)
                                                                                                                         ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTc                                                                                                                      <  1:166750/251‑1 (MQ=255)
                                                                                                                              tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCa                                                                                                                  <  1:309944/250‑1 (MQ=255)
                                                                                                                               gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTGTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTGCATTGTCGCCAt                                                                                                                 <  2:301335/250‑1 (MQ=255)
                                                                                                                                agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCCGAGAGGGGTTATCTTGCCTGACTGAACTGGCGGTTAAGGTTTCGCGTACCCGGCGGTTTGTTATTGCGGATGATGATATTTAAGATCCGCTGGTTgggttcttttcgccctcc                                                                                                               >  1:18590/1‑235 (MQ=255)
                                                                                                                                   gTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGGGAAGTGAGCGCAGAGAGGGATTATCTTTCCTGACTGTACTGGCGATTTAGTTTCCCCGCACCCGGCGGTTTGTTATTTGGGGTGATGGTATTTAAGCTCGGCTGATTGGTTTTTTTTCGCCATCAccc                                                                                                            >  2:245604/1‑250 (MQ=255)
                                                                                                                                       aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGTTCATTTTTTATTTTTCTTCCATTAGAAAGGGGGGCAGAGAGGGATTTTCATGCCGGACTTAAATAGCGATTAAACTTCCGCGCAACCGGCCTTTAGTTATTGCGGATTTTGATATTTGAGCTCCACTGATTGCGTTCTTTTTCCCATCACCGCtcg                                                                                                        >  1:268170/1‑249 (MQ=255)
                                                                                                                                              aaGTTTGCGGTTAACCCTTATAAGCTCCTTTCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCTGTCATTTTTTCGTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCTCGTACCCGGCGTTTAGTTATTGCGGATGATGATAGTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTTGACGgtgt                                                                                                 <  2:196118/251‑1 (MQ=255)
                                                                                                                                                        gTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGAGAGCGCAGAGAGGGATTATCTTGCCTGACTGGAATGGGGGTTAAGGTTCCGCGGACCCGGCGGTTTGGTATTTCGGGTGGTGATATTGAAGCCCCGCTGATTGGTTTATTTTCGCCATCACCGCTGGaagggggatttagtaa                                                                                        >  1:302151/1‑235 (MQ=255)
                                                                                                                                                                     gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGTTTTTTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGGTTAGTTATTGCGGGTGGTGAAATTGAAAGTCGGCTGATTGGTTCATTGTTGCCCTCACCGCTGGGCGGGGGTTTAAGTAAAGCGGCAACGCTgg                                                                          >  2:67354/1‑251 (MQ=255)
                                                                                                                                                                                 taCCAGTCAGGAGTCATTTCAGGATTCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCCCCGCGGGACGGTGTATTAAGTaa                                                                                        <  2:117948/225‑1 (MQ=255)
                                                                                                                                                                                 taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTTCCTGACTGCACCGGGGAGTAAGTTTTCGCGTAACCGGCGTTTTGTTATTGCGGGTGATGGTATTGAAGCTCGGCTGTTTGGTTCATTTTCGCCCTCACCCCTGGACGGGGTATTTAGTaa                                                                                        >  1:117948/1‑225 (MQ=255)
                                                                                                                                                                                       tCAGGAGACATTTCAGGATGCCATGTCGCGGAGCTCGTTGGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGGACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCTATCACCGCTGGACGGGGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAggta                                                         <  2:386151/250‑3 (MQ=255)
                                                                                                                                                                                        cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAAc                                                        <  2:18590/250‑1 (MQ=255)
                                                                                                                                                                                        cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGAAATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGGTATCTTGGCTGTCTTAACTCGCGCTTTAGTTTTCCAGCTGCCCGGGGTTTAGTTGTGCGGGTGAGGATATTGAAACTCGGCTTGGTTGGTTTTTTTAGCCCTCCTCGCTGGGCGGGGTATTATGTAAAGCGTCTACGGTTGCGAATTTTCATCAGGGAct                                                       >  2:255343/1‑251 (MQ=255)
                                                                                                                                                                                                aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGt                                               <  1:257502/251‑1 (MQ=255)
                                                                                                                                                                                                    cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCAtt                                           <  1:52180/251‑1 (MQ=255)
                                                                                                                                                                                                               gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTa                                 <  1:245604/250‑1 (MQ=255)
                                                                                                                                                                                                                  gcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGGGGGCAGAGAGGGATTATCTTGCCTGACGGAACTGGCGATTAAGTTTCCGCGTACCCGGCGGTTAGTTATTGCGGATGATGATATTGAAGCTCGTCTGATTGGTTCATTGTCGCCATCACCTCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAGCGCTTTTTGTCATTAAATGGTGTACGTc                             <  2:204371/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAAc         >  1:301003/1‑250 (MQ=255)
                                                                                                                                                                                                                                             ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTaa                                                                                        <  1:258310/165‑1 (MQ=255)
                                                                                                                                                                                                                                             ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTaa                                                                                        >  2:258310/1‑165 (MQ=255)
                                                                                                                                                                                                                                             ntctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATAAGt  >  2:187781/2‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                            
GCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGT  >  NC_000913/4604105‑4604592

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