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.

GCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTG‑GTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGT  >  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)
                                                             cacaGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGTCGATTAAGTTT‑Cc‑gcg                                                                                                                                                                                    <  2:123169/251‑1 (MQ=255)
                                                                                       cGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTTTAGTTATTGCGgatg                                                                                                                                                          <  2:193082/251‑1 (MQ=255)
                                                                                         cACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcgc                                                                                                                                                                                                                                     <  1:192489/176‑1 (MQ=255)
                                                                                                aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTg                                                                                                                                                 <  1:62006/251‑1 (MQ=255)
                                                                                                   aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGCGAGGGGTTATTTTTCCTGACTGAACTGGCGGTTAAGTTT‑TC‑GCGTACCCCGCGTTTATTTTTTGCGGGTTATGATATTGa                                                                                                                                                >  2:10654/1‑249 (MQ=255)
                                                                                                     aaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGGGAGCGCAGAGAGGGATTATTTTGCCTGACTGAACTGGCGATTAAGTTT‑TC‑GCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCt                                                                                                                                            >  1:201172/1‑251 (MQ=255)
                                                                                                          gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGTGAAGTGAGCGCAGAGGGGGGTTATCTTTCCTGACTGTACTGGCGATTTAGTTT‑TC‑GCGGTCCCGGCGGTTTGTTATTTCTGGTGATGGTATTTTAGCTCGGcc                                                                                                                                       >  2:1629/1‑250 (MQ=255)
                                                                                                             aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTTCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTTTTGTTATTTGGGGTGATGATATTGAAGCTCGGCT‑GAt                                                                                                                                   >  1:4599/1‑251 (MQ=255)
                                                                                                                agtagtATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTg‑g                                                                                                                               <  1:315650/251‑1 (MQ=255)
                                                                                                                         ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTG‑GTTCATTGTc                                                                                                                      <  1:166750/251‑1 (MQ=255)
                                                                                                                              tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTG‑GTTCATTGTCGCCa                                                                                                                  <  1:309944/250‑1 (MQ=255)
                                                                                                                               gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTGTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTG‑GTGCATTGTCGCCAt                                                                                                                 <  2:301335/250‑1 (MQ=255)
                                                                                                                                agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCCGAGAGGGGTTATCTTGCCTGACTGAACTGGCGGTTAAGGTT‑TC‑GCGTACCCGGCGGTTTGTTATTGCGGATGATGATATTTAAGATCCGCT‑GGTTg‑ggttcttttcgccctcc                                                                                                               >  1:18590/1‑235 (MQ=255)
                                                                                                                                   gTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGGGAAGTGAGCGCAGAGAGGGATTATCTTTCCTGACTGTACTGGCGATTTAGTTT‑CC‑CCGCACCCGGCGGTTTGTTATTTGGGGTGATGGTATTTAAGCTCGGCT‑GATTG‑GTTTTTTTTCGCCATCAccc                                                                                                            >  2:245604/1‑250 (MQ=255)
                                                                                                                                       aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGTTCATTTTTTATTTTTCTTCCATTAGAAAGGGGGGCAGAGAGGGATTTTCATGCCGGACTTAAATAGCGATTAAACTT‑CC‑GCGCAACCGGCCTTTAGTTATTGCGGATTTTGATATTTGAGCTCCACT‑GATTGCGTTCTTTTTC‑CCATCACCGCtcg                                                                                                        >  1:268170/1‑249 (MQ=255)
                                                                                                                                              aaGTTTGCGGTTAACCCTTATAAGCTCCTTTCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCTGTCATTTTTTCGTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑TCGTACCCGGCGTTTAGTTATTGCGGATGATGATAGTGAAGCTCGGCT‑GATTG‑GTTCATTGTCGCCATCACCGCTTGACGgtgt                                                                                                 <  2:196118/251‑1 (MQ=255)
                                                                                                                                                        gTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGAGAGCGCAGAGAGGGATTATCTTGCCTGACTGGAATGGGGGTTAAGGTT‑CC‑GCGGACCCGGCGGTTTGGTATTTCGGGTGGTGATATTGAAGCCCCGCT‑GATTG‑GTTTATTTTCGCCATCACCGCTGGaagggggatttagtaa                                                                                        >  1:302151/1‑235 (MQ=255)
                                                                                                                                                                     gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGTTTTTTTGCCTGACTGAACTGGCGATTAAGTTT‑TC‑GCGTACCCGGCGGTTAGTTATTGCGGGTGGTGAAATTGAAAGTCGGCT‑GATTG‑GTTCATTGTTGCCCTCACCGCTGGGCGGGGGTTTAAGTAAAGCGGCAACGCTgg                                                                          >  2:67354/1‑251 (MQ=255)
                                                                                                                                                                                 taCCAGTCAGGAGTCATTTCAGGATTCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTG‑GTTCATTGTCGCCATCCCCGCGGGACGGTGTATTAAGTaa                                                                                        <  2:117948/225‑1 (MQ=255)
                                                                                                                                                                                 taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTTCCTGACTGCACCGGGGAGTAAGTTT‑TC‑GCGTAACCGGCGTTTTGTTATTGCGGGTGATGGTATTGAAGCTCGGCT‑GTTTG‑GTTCATTTTCGCCCTCACCCCTGGACGGGGTATTTAGTaa                                                                                        >  1:117948/1‑225 (MQ=255)
                                                                                                                                                                                       tCAGGAGACATTTCAGGATGCCATGTCGCGGAGCTCGTTGGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGGACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTG‑GTTCATTGTCGCTATCACCGCTGGACGGGGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAggta                                                         <  2:386151/250‑3 (MQ=255)
                                                                                                                                                                                        cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTG‑GTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAAc                                                        <  2:18590/250‑1 (MQ=255)
                                                                                                                                                                                        cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGAAATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGGTATCTTGGCTGTCTTAACTCGCGCTTTAGTTTTCCAGC‑TGCCCGGGGTTTAGTT‑GTGCGGGTGAGGATATTGAAACTCGGCTTGGTTG‑GTT‑TTTTTAGCCCTCCTCGCTGGGCGGGGTATTATGTAAAGCGTCTACGGTTGCGAATTTTCATCAGGGAct                                                       >  2:255343/1‑251 (MQ=255)
                                                                                                                                                                                                aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTG‑GTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGt                                               <  1:257502/251‑1 (MQ=255)
                                                                                                                                                                                                    cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTG‑GTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCAtt                                           <  1:52180/251‑1 (MQ=255)
                                                                                                                                                                                                               gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTG‑GTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTa                                 <  1:245604/250‑1 (MQ=255)
                                                                                                                                                                                                                  gcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGGGGGCAGAGAGGGATTATCTTGCCTGACGGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGGTTAGTTATTGCGGATGATGATATTGAAGCTCGTCT‑GATTG‑GTTCATTGTCGCCATCACCTCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAGCGCTTTTTGTCATTAAATGGTGTACGTc                             <  2:204371/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTG‑GTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAAc         >  1:301003/1‑250 (MQ=255)
                                                                                                                                                                                                                                             ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTG‑GTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTaa                                                                                        <  1:258310/165‑1 (MQ=255)
                                                                                                                                                                                                                                             ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTG‑GTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTaa                                                                                        >  2:258310/1‑165 (MQ=255)
                                                                                                                                                                                                                                             ntctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTG‑GTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATAAGt  >  2:187781/2‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                            
GCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTT‑CC‑GCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTG‑GTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGT  >  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: ‑