Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,399,929 A→T L17Q (CTG→CAG)  mreC ← cell shape determining protein MreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,399,9290AT98.1% 161.1 / NA 53L17Q (CTG→CAG) mreCcell shape determining protein MreC
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base T (25/26);  minor base G (0/1);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAA  >  NC_000913/3399687‑3400173
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
tGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAgg                                                                                                                                                                                                                                              >  1:9321/1‑251 (MQ=255)
      gTCCAAGCAGCAGCAGTTCAATGATTTTCAGCACCAGTACCTGACGGAACGCCCGGTTTTCAAGTTCTAATTGATCACGCGAGGCAAGCGTCTGCGGTACGCCATCCAGCAATTAACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATCAGTACGGATTTGACTGAACGACCCCACGCGGCTGTCGGCAATAATAATGCCGAACCCCACCTGCCCCGCCTGAATAAGGcaaagc                                                                                                                                                                                                                                        >  2:85227/1‑246 (MQ=255)
           aGCACCAGCAGTCCCCTGTTTTCCAGCACCAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCTAGCGTCTGCGATACTCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAg                                                                                                                                                                                                                                   <  1:404949/251‑1 (MQ=255)
            gCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGc                                                                                                                                                                                                                                  <  1:284708/251‑1 (MQ=255)
                cagcagTTCACTGTTTTTCAGCAACAGTTCCTGCCGTAACGCCCGGTTTTCCAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAACGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCCCCTGAATAAGGCGACTCTGTAGCGACg                                                                                                                                                                                                                              >  2:227023/1‑251 (MQ=255)
                cagcagTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGAc                                                                                                                                                                                                                               <  1:151637/250‑1 (MQ=255)
                           gtttttcCAGCAACAGTTCCTGCACTAACGCCCGTTTTTCAAGTTCTAATTGGTCACGCAAGGCCAGCGTCGTCGATCCGCCATCCATAAAGTCACTAGGAGCATTGGAAACAAATTAGAAAGGACTGACTGCGGTATCCATATAAGTCCGGATTTGACTGAACGTCCCCAGGCGGCTTTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAACGACGGTCCTCGGCTa                                                                                                                                                                                                                   <  2:442031/249‑1 (MQ=255)
                            gTTTTTCAGCAACCGCTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGAGGCTGTCGGCAATAATAATGACGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTa                                                                                                                                                                                                                   >  1:296537/1‑250 (MQ=255)
                                 ccAGCAACAGTTCCTGCCGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCGTCGATACGCCATCCAGCAATGCACGAGGAGCATTGGAAACAAAGGAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAAtt                                                                                                                                                                                                             <  2:24841/250‑1 (MQ=255)
                                  cAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTg                                                                                                                                                                                                            <  1:255503/251‑1 (MQ=255)
                                          gTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGCGGAGCATTGGAAACAAAGGAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCCCCTGAATAAGGCGAATCTGTAGCGCAGGGCCACGGCACAAAATTGACTTcgt                                                                                                                                                                                                     >  2:328210/1‑248 (MQ=255)
                                              cTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATTAGtt                                                                                                                                                                                                >  2:273492/1‑251 (MQ=255)
                                                cgacctAACGCCCGGTTTTCAAGTGCTAATAGGTCCCGCGAGGCCAGCGTCAGCTATACGCCAGCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCGTCATAAGTTa                                                                                                                                                                                               <  2:86843/246‑1 (MQ=255)
                                                   cGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGc                                                                                                                                                                                            <  1:173358/251‑1 (MQ=255)
                                                   cGTAACGCCCGGTTTGCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGc                                                                                                                                                                                            <  1:277512/251‑1 (MQ=255)
                                                      aaCGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTa                                                                                                                                                                                               >  1:86843/1‑245 (MQ=255)
                                                      aaCGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGACAATAATAATGCCGAGCGCCACCAGCACCGCATGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCCTAAGTTATGCGt                                                                                                                                                                                          >  1:467450/1‑250 (MQ=255)
                                                                    aaGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGAc                                                                                                                                                                           >  1:185812/1‑251 (MQ=255)
                                                                           aatgggTCACGCGTGGCCAGCGTCTGCGATACGCCATCCATCAATTCACGAGTAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGCCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTATCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACcaggc                                                                                                                                                                      <  2:296537/245‑1 (MQ=255)
                                                                                  cACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCtt                                                                                                                                                                                                        >  2:130472/1‑208 (MQ=255)
                                                                                  cACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCtt                                                                                                                                                                                                        <  1:130472/208‑1 (MQ=255)
                                                                                      cgAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGTCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGac                                                                                                                                                         <  1:352552/251‑1 (MQ=255)
                                                                                       gaagCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCATCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTCACAATTGGCTTCTTAAGTTATGCGTATTCTCGTATAAAACCAGGCAGGGGAAACCGac                                                                                                                                                         >  2:305652/4‑250 (MQ=255)
                                                                                       gAGGCAAGCGTCTGCCATCCGCCACCCAGCAATGCCCGAGGAGAATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCCCCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTAGGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGaca                                                                                                                                                        <  2:8475/251‑1 (MQ=255)
                                                                                        ggggcAGCGTCTGCGATACGCCAGCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGATGGACGGCGGTATCCATTTAAGTACGGATTGGACTGAACGTCCCCAGGCGGCTGTCGGTAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGacac                                                                                                                                                       <  2:409868/247‑1 (MQ=255)
                                                                                                   tGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCcctgcct                                                                                                                                            >  1:233690/1‑251 (MQ=255)
                                                                                                          cGCCATCCAGCAATGCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCg                                                                                                                                     <  1:348093/251‑1 (MQ=255)
                                                                                                                    caagtCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGAGTTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCTAGCGCCACCAGCACCGCCTGAATATGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATGCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTTCTCTTc                                                                                                                           <  1:458098/247‑1 (MQ=255)
                                                                                                                     aaTTCACGAGGAGAATTGGAAACAAAGTCGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACAGACCGTCTCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGACGAATCTGTAGCGACGGGCCACGGATAAAAATTGGCTTCATAAGTTATGCGAATTCACGTATCAGACCAGGCAGGGCAAACAGACAATTCCCCTGCCTTCCTCCGATTActatcac                                                                                                                          >  2:32214/1‑246 (MQ=255)
                                                                                                                        tcccggggagCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCt                                                                                                                        <  2:185812/244‑1 (MQ=255)
                                                                                                                        tCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCt                                                                                                                        <  1:468341/250‑1 (MQ=255)
                                                                                                                             aggagCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGCCGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAAc                                                                                                                    >  1:75540/1‑249 (MQ=255)
                                                                                                                             aggagCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGAt                                                                                                                                   >  1:442276/1‑234 (MQ=255)
                                                                                                                             aggagCATTGGAAACAAAGGAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGAt                                                                                                                                   <  2:442276/234‑1 (MQ=255)
                                                                                                                                   aTTGGAAACAAGGGAGAAAGGACTGACGGCGGTATCAATATAAGTACGGACTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAAACTGTAGCGGCGGGCCACGGATAAAAATTGGCATCATAAGTTCTGCTTATTCTCGTAGCAGACCCGGCCGGGCAAACAAACTCTTCCCCTTCCTTCTTCCGGTTATTCTTCCCTGCACATGTcgcc                                                                                                            >  2:172655/1‑251 (MQ=255)
                                                                                                                                        aaccAAAGTAGAAAGGACTGACGGCGGTACCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAAGAATGCCGAGGGCCCCCAGCCCCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCTTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGt                                                                                                       <  1:469676/248‑1 (MQ=255)
                                                                                                                                         aaCAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCCCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTg                                                                                                      >  1:40923/1‑251 (MQ=255)
                                                                                                                                           cAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTCTCAGACCAGGCAGGGGAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTga                                                                                                     >  1:107222/1‑249 (MQ=255)
                                                                                                                                            aaaGGAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGACATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGACAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCCGAACAGGTCACCGCCGTGctt                                                                                                   >  2:217550/1‑249 (MQ=255)
                                                                                                                                               gtagaaagtccggacgggtggagccATTTAAGTGCGCATTTGATGGAACGTCCCAATGCGGTGGTCGGAAATAATAATTCCGATCGCCACCAGCTCCGCCTGAATAAGGCGATTCTGTAGCGACGGTCCACGTCTAAAACTTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGTCACCTCCCCTGTCTGCATCCGATTACTCTGCGCTGAACAGGTCGCCTCCGTGCATGTc                                                                                                <  2:296229/228‑1 (MQ=255)
                                                                                                                                                  gAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCCCCGCCTGAATAAGGCGAATCTGTAGCAACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACCCTTCCCCTGCCTGCCTCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGAAAGTCGAt                                                                                             >  1:363251/1‑251 (MQ=255)
                                                                                                                                                    aaGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGTATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAAAACAGGCAGGGTAAACAGACCCTTCCCCTGCCTGCATCCGCTTCCTCTTCGCTGAACAGGTCGCCCCCGTGCCTGTCGATc                                                                                            >  2:404984/1‑250 (MQ=255)
                                                                                                                                                                       cATATAAGTGCGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGAGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATTCGTATTCTCGTATCATACCAGGCAGGGAAAACAGCCCCTTCCCCTGCCTGCATCCGCTTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCGTTTCCAGCGCTTTtccgcc                                                                        >  2:199630/1‑251 (MQ=255)
                                                                                                                                                                       cATATAAGTACGGATGTGACTGAACGTCCCCAGGCGTCGGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATGGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTGCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgccgcc                                                                        <  2:851/251‑1 (MQ=255)
                                                                                                                                                                          agaagtccgtagttgaccggacgaccccAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACTAGCACCGCCTGAAAAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGAAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCAc                                                                      <  1:35800/226‑1 (MQ=255)
                                                                                                                                                                                                  ccccAGGCGGCTGGCGGCAAAAATAAGGCCGAGCGCCACAAGCAACGCCGGAAGAAGGCGAATCTGTAGCGCCGGGCCACGGCTATAAATGGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCAGCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTc                                             <  2:231165/251‑1 (MQ=255)
                                                                                                                                                                                                  cccagggcgaatGTCGGCAAAAAAAATGCCGAGCGCCCCAAGCACCGACTAAAAAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATCCTCTTATCAGACCAGGCAGGGTAAACAGACACGTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTc                                             <  2:9321/240‑1 (MQ=255)
                                                                                                                                                                                                                 ggCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGAc                              >  1:36319/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gcCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGCACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTttcttac           >  2:243612/1‑249 (MQ=255)
                                                                                                                                                                                                                                    gcCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTAGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACAATTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTcc           >  2:439108/1‑251 (MQ=255)
                                                                                                                                                                                                                                      caccaGCACGGCCTAAATACGGCGAGTCTGTAGCGACCGGCCACGGCTAAACATTGGCTTCATCAATTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGGACCGGTCGCCGCCGTTCATGTCGCTCATTTCCAGCGCTTTGCCGCCCCAGCGCGCCCCACAGGCCAGCGGGTCTTCAGCAACAACACCTGGAATGCCGGTttcttacgt         >  2:236348/1‑247 (MQ=255)
                                                                                                                                                                                                                                          aGCACCGACTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTAGCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTTTTCCATTAAc     <  2:199421/251‑1 (MQ=255)
                                                                                                                                                                                                                                             cccGCCTGAATAAGGCGAACCTGTAGCGACGGGCCACGGCTAAAAATTGCCTTCATAAGTTATGCGTATTCTCGTATCAGAACAGGCAGGGTAAACATACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACaaa  <  1:1449/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
TGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAA  >  NC_000913/3399687‑3400173

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