Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,224,239 A→T 18.4% intergenic (+332/+40) ycgI → / ← minE pseudogene/cell division topological specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,224,2390AT18.4% 338.4 / 60.6 146intergenic (+332/+40)ycgI/minEpseudogene/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base A (51/67);  new base T (12/15);  total (63/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

TGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGC  >  NC_000913/1224140‑1224337
                                                                                                   |                                                                                                  
tgtAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCGTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAgaga                                                                                                    <  1:4079073/100‑1 (MQ=255)
   aaaGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGAt                                                                                                   <  1:1595430/98‑1 (MQ=255)
    aaGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAAc                                                                                                <  2:3021192/100‑1 (MQ=255)
      gggggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAActc                                                                                              >  2:2984899/1‑100 (MQ=255)
       ggggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAActct                                                                                             <  1:3078310/100‑1 (MQ=255)
       ggggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGAGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAActc                                                                                              >  2:2460074/1‑99 (MQ=255)
          gAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGcc                                                                                          <  2:355229/100‑1 (MQ=255)
           aGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCt                                                                                         <  2:2004448/100‑1 (MQ=255)
           aGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCt                                                                                         <  2:4149965/100‑1 (MQ=255)
           aGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCt                                                                                         <  1:210961/100‑1 (MQ=255)
            ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCtt                                                                                        <  1:567096/100‑1 (MQ=255)
            ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCtt                                                                                        >  1:379253/1‑100 (MQ=255)
            ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCtt                                                                                        <  2:3467489/100‑1 (MQ=255)
            ggTGGAAAGAGCAACGACGAGTATCAGCAATAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCATAGATAACTCTGCCtt                                                                                        <  2:842234/100‑1 (MQ=255)
            ggTAGAAATAGCCATGAGGACTATCAACCAAAATACTCGCCGCTTTACCACAATGTAGATGAAAAGAAAGAAAAACTCAAGGCAGAGTTATCTCTGCCt                                                                                         <  2:3436057/99‑1 (MQ=255)
             tctaaaaTAGCTAACAGAAGTAACAGCAACAAAAATCACCGCTTTACCACAACGACGATGAAAGATATGACAACCTCACTTTAGAGCTATCTCTGCCTTg                                                                                       <  2:1014378/96‑1 (MQ=255)
             gTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTg                                                                                       <  1:1119765/100‑1 (MQ=255)
             gTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTg                                                                                       <  2:3209702/100‑1 (MQ=255)
             gTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTg                                                                                       <  2:3415907/100‑1 (MQ=255)
             gTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTg                                                                                       >  1:3328643/1‑100 (MQ=255)
              tGGAAATAGCAATGAGGAGGATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGa                                                                                      <  2:1377873/100‑1 (MQ=39)
               ggAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGaa                                                                                     >  2:2685653/1‑100 (MQ=255)
               ggAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGaa                                                                                     >  1:1446418/1‑100 (MQ=255)
               ggAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGaa                                                                                     <  1:794051/100‑1 (MQ=255)
                gAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAg                                                                                    <  1:3978244/100‑1 (MQ=255)
                gAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAg                                                                                    >  1:3301102/1‑100 (MQ=255)
                gAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAg                                                                                    <  2:2077225/100‑1 (MQ=255)
                gAAATAGCAATAAGGAGTATCAGAAAGAATACTCGCAGTTTTAACACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCGTGagt                                                                                    <  2:3806564/100‑3 (MQ=255)
                  aaTAGCAATGATGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTg                                                                                       <  2:2033559/95‑1 (MQ=255)
                   aTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGagtttt                                                                                 <  2:1108095/100‑6 (MQ=38)
                    tAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaa                                                                                >  2:3973068/1‑100 (MQ=255)
                      gCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGagtttttca                                                                              <  1:3962213/100‑9 (MQ=37)
                       cAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATg                                                                             <  1:4151388/100‑1 (MQ=255)
                       cAACGAAGAGTAAAAGCAAAAAAAATCGCCGCTAGACCACAACCTGGATAAAAAAAAAAAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                                              <  2:679609/99‑1 (MQ=255)
                         aTGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgcg                                                                           <  1:961806/100‑1 (MQ=255)
                          tGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgcgc                                                                          <  2:1470049/100‑1 (MQ=255)
                           gaggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCt                                                                         >  1:2039104/1‑100 (MQ=255)
                             ggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCttt                                                                       >  2:144378/1‑100 (MQ=255)
                             ggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCttt                                                                       >  1:2072317/1‑100 (MQ=255)
                              gagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCtttt                                                                      <  1:2096726/100‑1 (MQ=255)
                               agTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                                                     <  2:3625936/100‑1 (MQ=255)
                               agTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                                                     >  2:2775744/1‑100 (MQ=255)
                               agTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                                                     >  1:2970648/1‑100 (MQ=255)
                               agTATCAGCAAGAATACTCGCCGCGTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                                                     <  2:254330/100‑1 (MQ=255)
                                gTAGCAGCAAGAATACTCGCCGCTTTACCACAACGTGTATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTAc                                                                    <  2:129011/100‑1 (MQ=255)
                                     aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                                               >  2:1632555/1‑100 (MQ=255)
                                      gCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCggg                                                              >  1:1057432/1‑100 (MQ=255)
                                        aaGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCt                                                            >  2:2090276/1‑100 (MQ=255)
                                         aGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGc                                                             <  2:562757/98‑1 (MQ=255)
                                          gAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTa                                                          <  2:1898595/100‑1 (MQ=255)
                                           aaTACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAt                                                         >  2:2329783/1‑100 (MQ=255)
                                            aTACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAtt                                                        <  1:2902817/100‑1 (MQ=255)
                                             tACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAttt                                                       >  1:1628075/1‑100 (MQ=255)
                                              aCTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTc                                                      <  1:828294/100‑1 (MQ=255)
                                              aCTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTc                                                      >  2:616840/1‑100 (MQ=255)
                                               cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                     >  1:2369593/1‑100 (MQ=255)
                                               cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                     >  1:1837828/1‑100 (MQ=255)
                                                tCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAg                                                    <  2:895574/100‑1 (MQ=255)
                                                 cgccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGc                                                   >  2:2480519/1‑100 (MQ=255)
                                                  gccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGct                                                  >  1:2307240/1‑100 (MQ=255)
                                                   ccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGctc                                                 >  1:342612/1‑100 (MQ=255)
                                                      cTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttc                                              <  2:2202484/100‑1 (MQ=255)
                                                        ttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTg                                            <  1:221841/100‑1 (MQ=255)
                                                           ccACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCtt                                         >  2:806994/1‑100 (MQ=255)
                                                           ccACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCtt                                         >  1:3320994/1‑100 (MQ=255)
                                                            cacaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTc                                        <  1:12381/100‑1 (MQ=255)
                                                             acaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTcc                                       <  2:937202/100‑1 (MQ=255)
                                                              caACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATACGCTGTTACGACGGGCTTATTTCAGCTCTTCTGCTTcc                                       >  2:3918384/1‑99 (MQ=255)
                                                               aaCGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCgg                                     <  1:2246259/100‑1 (MQ=255)
                                                                    ggATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAgg                                <  2:1702540/100‑1 (MQ=255)
                                                                      aTGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAgg                                <  2:1615426/98‑1 (MQ=255)
                                                                       tGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCa                             >  1:621741/1‑100 (MQ=255)
                                                                       tGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCa                             <  1:1237135/100‑1 (MQ=255)
                                                                          gagGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGt                          <  1:560249/100‑1 (MQ=255)
                                                                           agGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGtt                         <  1:3641063/100‑1 (MQ=255)
                                                                            gggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTc                        >  1:3839071/1‑100 (MQ=255)
                                                                            gggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTc                        <  2:3269654/100‑1 (MQ=255)
                                                                               catttatcttCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACATCAGGCTTATTTCAGCTCTTCTGCTTCTGGTAAGGTCACGTTCAGc                     >  2:797230/10‑100 (MQ=25)
                                                                                 gAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTc                   >  1:2489200/1‑100 (MQ=255)
                                                                                 gAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTc                   >  2:2862904/1‑100 (MQ=255)
                                                                                    aagcTCAAGGCAGAGATAACTCTGGCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTACAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAg                <  2:970370/97‑1 (MQ=255)
                                                                                     tcttCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGa               >  1:3426951/4‑100 (MQ=39)
                                                                                     aaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGa               <  1:3918694/100‑1 (MQ=255)
                                                                                      aCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGaa              <  1:173041/100‑1 (MQ=255)
                                                                                       ttCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAAt             >  2:3739495/2‑100 (MQ=255)
                                                                                       ttCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAAt             >  2:3579352/2‑100 (MQ=255)
                                                                                        tCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATa            >  2:4190969/1‑100 (MQ=255)
                                                                                         cagggCAGAGTTATGTCTGTCTTGAAGATAAATGCGTTTTTACAGCGGGCTTATTTCAGCTGTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAg           >  2:3591809/4‑100 (MQ=255)
                                                                                         cAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAg           >  1:1494872/1‑100 (MQ=255)
                                                                                          aaGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa          >  1:1784150/1‑100 (MQ=255)
                                                                                          aaGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa          >  2:1880210/1‑100 (MQ=255)
                                                                                           aGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaa         <  1:1610963/100‑1 (MQ=255)
                                                                                            ggCAGAGATACCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaaa        <  2:2861384/100‑1 (MQ=255)
                                                                                             gCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAat       >  1:2308002/1‑100 (MQ=255)
                                                                                             gCAGAGATACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAata      >  2:1678185/1‑100 (MQ=255)
                                                                                              cAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaaa        >  2:3423779/1‑98 (MQ=255)
                                                                                              cAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAata      <  1:1836791/100‑1 (MQ=255)
                                                                                                gagaTAACTGTGACTTGCAGAGAAGTGCGTTTGTACAGCGGGCTGATTTCAGCTCTTATGCTTCCGGTAAGGTCACGTTCAGCTCAAGACTAGAAATATc    <  2:2524789/100‑1 (MQ=255)
                                                                                                gagaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATc    >  1:1687411/1‑100 (MQ=255)
                                                                                                  gaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGc  >  2:4141621/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
TGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGC  >  NC_000913/1224140‑1224337

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