Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,224,0210AG37.1% 43.8 / 11.6 36intergenic (+572/+258)ymgJ/minEuncharacterized protein YmgJ/Z‑ring positioning protein MinE
Reads supporting (aligned to +/- strand):  ref base A (4/18);  new base G (13/0);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACGTGGAATATACCAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTAC  >  NC_000913/1223772‑1224271
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                               
aCGTGGAATATACCAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGgag                                                                                                                                                                                                                                                                <  1:659430/251‑1 (MQ=255)
           aCCAGGCAGTTATTAAAGCCTCATTCAGCAGTTCGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGaat                                                                                                                                                                                                                                                      <  1:732703/250‑1 (MQ=255)
                      aTTAAAGACTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGTCGGTGTGGAATTCCCGTGGAACGCGGTGGCTGTTGTGAAATGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGTAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGAGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTg                                                                                                                                                                                                                                           <  2:1038849/249‑1 (MQ=255)
                                   tCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAACCTCATCCGATTGGTGGACAATGAATAATAAAATTTTGGGGTTattaattt                                                                                                                                                                                                                           >  1:419464/1‑245 (MQ=37)
                                   tCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATCCCACGGGGGGACAAAGAATAATTAAAATGCGGGGGTatcattcc                                                                                                                                                                                                                           >  2:943688/1‑248 (MQ=37)
                                        agTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATCCTATTGGGGAAAAAAGAAAAATGAAAATGCGGGGGTATCATCTcccccg                                                                                                                                                                                                                      >  2:221614/1‑248 (MQ=37)
                                             gTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATc                                                                                                                                                                                                                 <  1:779555/251‑1 (MQ=255)
                                               ttAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAAACCAATCCCATGTGGGAACAATAAAAAATGAAAATGGGGGGGTattattctcccctatatta                                                                                                                                                                                                                >  2:1227009/1‑233 (MQ=25)
                                                 aGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCAt                                                                                                                                                                                                              <  1:943688/250‑1 (MQ=255)
                                                   cGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAAACTCATCCTATTGGGGAAAAATGAATAAAGAAATTGTGGGGGTGTCATCTCCCAGTATAACCATACt                                                                                                                                                                                                           >  1:724112/1‑251 (MQ=37)
                                                             ggggtggggtATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACCCCTCCTATGGGGGAAAAAAGAATAATGAAATTGCGGGGGTATCATCTCCCCGTatatt                                                                                                                                                                                                                 >  1:466008/1‑234 (MQ=37)
                                                             ggggtggggtATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACGGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATc                                                                                                                                                                                                                 <  2:466008/235‑1 (MQ=255)
                                                                  ggggtATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAACCTCATCCTATGTGGGGACAATGACTAATGAAAATTCGGGGGTTTTCTCTCCCACTATATTCACTCTTACAATAAGGTTAtt                                                                                                                                                                                            >  1:48287/1‑251 (MQ=25)
                                                                     gtATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGgaga                                                                                                                                                                                                                                                               <  1:395477/183‑1 (MQ=255)
                                                                             aTGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACGGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAg                                                                                                                                                                                <  2:419464/251‑1 (MQ=255)
                                                                               gggtccGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGTTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCTCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATcc                                                                                                                                                                                                                <  2:790715/214‑1 (MQ=255)
                                                                                  gCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAAACCCCCTCTTTGGGGGGAAAAATAATAATGAAATTTCGGGGGTTTTTTTTCCCCGAATATCCATACTAAAAAAAAGGTTAttttttctaacaagctc                                                                                                                                                                             >  2:213408/1‑236 (MQ=25)
                                                                                   ccgttgtgGAATCCCCGTTGAACGCGGTGGCTGGTGAGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTGCTGCGCATTTTTCGTTTACAGGATGTGGAAAGGGAAAATTCTGGGGCAAAAAAAGCTCGCCACTTACGTCGGGAAACCTCATCCTATGTGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCATTATATCCATACTAACAATAAGGTTATTTACTCATCCAGGCATaa                                                                                                                                                                          <  2:1024493/247‑1 (MQ=255)
                                                                                               ccccGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCAGCCTATGGGAGAACAATGAATAATGaaa                                                                                                                                                                                                                                               <  1:286170/173‑1 (MQ=255)
                                                                                                                         aaCTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGTGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCCTCCTATGGGGGAACAAAGAATAATTAAAATGCGGGGGTATCATCTCCCCCGATACCCCAAATAAACATAAAGTTTTTTACTCAACCCGGCCAAAACATTTTTTTTTGTTCGGGGGAACAGCCTTAAGgtggg                                                                                                                                  >  1:1220676/1‑249 (MQ=255)
                                                                                                                              gTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAAACACATCCTCTGGGGGAACAAAGAATAATGAAATTTCGGGGGTATTATTTCCCCATATATCCATAATAACAATAAAGTTATTTTTTCAACCCAGCACAAACATTTTTTTTTTTGTGGGGGTACAGCCTTTAGGGGGAAAggg                                                                                                                             >  2:661134/1‑251 (MQ=255)
                                                                                                                              gTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCACCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAggg                                                                                                                             <  1:795376/251‑1 (MQ=255)
                                                                                                                              gTCGTTCTGACAATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAACAAAAGCCCGCCAGTTACGGCGGGCAACCTCCCCCTTTGGGGAAAAAAAGAAATATTAAAATGCGGGGGTTTTCTCTCCCCCTATCTCCATACTAACAATAAGGGTTTTTACTCACACAAGCCTAAACCATTTTTTTTTTGCGTGGGGACAATCTTAAGGGGTAAAggg                                                                                                                             >  2:795376/1‑251 (MQ=255)
                                                                                                                                cGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACCCATCCTATGGGGGAACAATGAATAATGAAATTGGGGGGTTATCATTTCCCAGAATATCCATACTAAAAATAAGGTTATTTACTTAACCAAGCATAAACATTTTTTTTTTTGCGGGGGAACAGCCTTAAGGGGGAAAggggg                                                                                                                           >  1:252274/1‑251 (MQ=255)
                                                                                                                                   tCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAAGCCCGCCAGTTACGGCGGGAAAACCCATACTATAGGAGAAAAATGAATAATTAAATTGCGGGGGTATCATTTCCCAGTAAATCCATAATAACAATAAGGGTGTTTAATCCACCAAGGCTAAAAAATTTTTTTTTTGCGGGGGGAAAGCCTTTAGGGGTAAAggggg                                                                                                                           >  1:263501/1‑249 (MQ=255)
                                                                                                                                    cTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGGGAAAAAAGAAAAATGAAATTGCGGGGGTATCATCTCCCAGCATATCCATACTAAAAATAAAGTTATTTTTTCCACCACGCACAAAAATTTTTTTTTTTTGGTGGGGACAGACTTTAGGGGGAAAggggggggg                                                                                                                       >  2:911806/1‑247 (MQ=255)
                                                                                                                                                      aCTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCa                                                                                                                                                                                  <  1:1175895/177‑1 (MQ=255)
                                                                                                                                                                      gTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACt                                                                                     <  1:221614/251‑1 (MQ=255)
                                                                                                                                                                           cGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTcgcc                                                                                 <  1:661134/250‑1 (MQ=255)
                                                                                                                                                                                   ggATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGGGAACAAAGAAAAATGAAATTGCGGGGGTATCATCTCCCAGTATATCCATACTAACAATAAGGGTATTTACTCAACCAAGCATAAACATTTTTTTTTTTGCGTGGGAACAGCCTTAAGGGGTAAAGGGGGAGGGGGGAATAGCAATGAGGAGTATCAGCAAGAATACACGCCGCTTTAcc                                                                         >  1:1106467/1‑250 (MQ=255)
                                                                                                                                                                                     atggtGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACcaca                                                                      <  2:48287/247‑1 (MQ=255)
                                                                                                                                                                                                            gCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTTTCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGaaaaa                                               <  1:213408/251‑1 (MQ=255)
                                                                                                                                                                                                                      gCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCa                                      <  1:1149227/251‑1 (MQ=255)
                                                                                                                                                                                                                         cGCCAGTTACGGCGGGAAACCTCATCCTATTGGAGAACAATGAATAATGAAATTACGGGGTTATCATCTCCCAGTATATTCATAATAACAATAAGGGTATTTACTCAACCAAGCATAAACATTTTTTTTTTTGCGTGGGAACAGCCTTTAAGGGGAAAGGGGGGGGGGGGAATAGCAATGAGGGGTATCAGCAAGAATAATCGCCGCTTTTCCACAACGGGGGTGAGAGGGGGGAAAAACTCaaaggag                                     >  1:14988/1‑244 (MQ=255)
                                                                                                                                                                                                                                       ggAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTg                     >  1:172620/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaa              <  2:252274/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   ccTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgc          >  2:188314/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      aTGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCt       <  2:1106467/251‑1 (MQ=255)
                                                                                                                                                                                                                                                           agaACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTAc  <  1:982058/251‑1 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                               
ACGTGGAATATACCAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTAC  >  NC_000913/1223772‑1224271

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