Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,224,239 A→T 15.1% 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,2390AT15.1% 354.6 / 50.5 140intergenic (+332/+40)ycgI/minEpseudogene/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base A (61/57);  new base T (14/7);  total (76/64)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

Only 100 of 146 total aligned reads displayed.

GTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCC  >  NC_000913/1224141‑1224338
                                                                                                  |                                                                                                   
gtAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGAt                                                                                                    >  1:3501755/1‑100 (MQ=255)
ggaaaGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAgaga                                                                                                     <  2:1463731/97‑1 (MQ=255)
  aaaGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATaa                                                                                                  >  1:2871751/1‑100 (MQ=255)
      ggggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAActct                                                                                              >  1:1198684/1‑100 (MQ=255)
          aGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCt                                                                                          <  1:3185618/100‑1 (MQ=255)
           ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGATAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCtt                                                                                         <  1:3377957/100‑1 (MQ=255)
           ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGACGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCtt                                                                                         >  2:2909268/1‑100 (MQ=255)
            gTGGAAATAGCAATGAGAAGTATCAGCAAGAATACTCGACGCTTTACCACAACGTGGATGAGAGCGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTg                                                                                        <  2:409616/100‑1 (MQ=255)
            gTGGAAATACCAATGAGGAGTATCAGCAAGCATACTCGACGCTTTACCACAACGTGGATGAGAGGGATGAAAATCTCAAGGCAGAGTTATCTCTTCCTTg                                                                                        <  2:2631686/100‑1 (MQ=255)
              ggAAATAGCAATGAGGAGTATCAGCAAGAATACTCGTCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGaa                                                                                      >  2:3450160/1‑100 (MQ=255)
               gAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGagt                                                                                     <  1:1558944/100‑3 (MQ=255)
               gAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAg                                                                                     <  1:803864/100‑1 (MQ=255)
               gAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAg                                                                                     <  1:1704763/100‑1 (MQ=255)
                aaaTAGCAATGAGGAGTATCAGCAAGAAGACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGagt                                                                                     <  2:3180930/99‑3 (MQ=38)
                  aTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATa                                                                                  <  1:1878407/100‑1 (MQ=255)
                  aTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATa                                                                                  >  1:461097/1‑100 (MQ=255)
                  aTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATa                                                                                  <  1:602090/100‑1 (MQ=255)
                     gCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                                               >  1:1141553/1‑100 (MQ=255)
                       aaTGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgc                                                                             <  1:1607247/100‑1 (MQ=255)
                         tGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgcgc                                                                           <  1:1653780/100‑1 (MQ=255)
                           aggataaTCAGCAAGATTACTCGCCGCTGCACCACAGCTCAGGTGGGATGGATGAAATCCTCAAGTTAGAGATAACTCTGCCTAAAAGAGAAATGCGCtt                                                                         <  2:1692385/94‑1 (MQ=255)
                           aggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCACGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCtt                                                                         <  1:2688893/100‑1 (MQ=255)
                           aggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCtt                                                                         >  1:2351133/1‑100 (MQ=255)
                             gagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCtttt                                                                       <  2:1931157/100‑1 (MQ=255)
                             gactaacAGAAAGAATACTCGCCTGTTTAACCCATCGTTGATAAAAAAAATAAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgcg                                                                            <  2:1705626/89‑1 (MQ=255)
                              agTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                                                      <  1:3463109/100‑1 (MQ=255)
                              agTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                                                      <  1:693536/100‑1 (MQ=255)
                                 aTCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                                                   >  2:1351552/1‑100 (MQ=255)
                                 aTCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                                                   <  1:983632/100‑1 (MQ=255)
                                 aTCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                                                   >  2:3174787/1‑100 (MQ=255)
                                  tCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGc                                                                  >  1:1349163/1‑100 (MQ=255)
                                   taGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCg                                                                 <  2:1800337/99‑1 (MQ=255)
                                    aGCAAGAATACTCGCCGCTTTACCCCAACGTGGATTAGAGGGATGAAAAACTCAATGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                                                >  1:3104946/1‑100 (MQ=255)
                                    aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                                                >  1:2808003/1‑100 (MQ=255)
                                    aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCga                                                                >  1:1800347/1‑99 (MQ=255)
                                      cAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGc                                                              >  1:1955493/1‑100 (MQ=255)
                                      cAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGc                                                              <  1:3174811/100‑1 (MQ=255)
                                       aaGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCt                                                             >  1:2057940/1‑100 (MQ=255)
                                       aaGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCt                                                             <  1:150905/100‑1 (MQ=255)
                                        aGAATACTCGCCGCTTTACCACAACGGGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCtt                                                            >  1:2332102/1‑100 (MQ=255)
                                           aTACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAtt                                                         <  1:609300/100‑1 (MQ=255)
                                            tACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAttt                                                        <  1:1236602/100‑1 (MQ=255)
                                              cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                      >  2:3740727/1‑100 (MQ=255)
                                              cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                      <  2:2894346/100‑1 (MQ=255)
                                              cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                      >  1:1142921/1‑100 (MQ=255)
                                              cTCGCCCCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                      <  2:495314/100‑1 (MQ=255)
                                                cgccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGc                                                    <  1:2943853/100‑1 (MQ=255)
                                                cgccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGc                                                    <  1:1811413/100‑1 (MQ=255)
                                                   cgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGctct                                                 >  2:1947823/1‑100 (MQ=255)
                                                   cgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGctct                                                 >  1:2398944/1‑100 (MQ=255)
                                                   cgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGctct                                                 >  1:2396536/1‑100 (MQ=255)
                                                   cgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGctct                                                 <  2:2656637/100‑1 (MQ=255)
                                                     cTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttc                                               <  1:3604527/100‑1 (MQ=255)
                                                       ttACCACAACGTGGATGCGAGGGATGAAAAACTCAGGGCAGAGATCACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCCGCTCTTCTg                                             <  2:834995/100‑1 (MQ=255)
                                                       ttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTg                                             <  2:486277/100‑1 (MQ=255)
                                                           cacaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTc                                         >  2:653609/1‑100 (MQ=255)
                                                               aCGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGt                                     <  1:3430045/100‑1 (MQ=255)
                                                                cGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTa                                    >  1:2533334/1‑100 (MQ=255)
                                                                   ggATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAg                                  <  1:201050/99‑1 (MQ=255)
                                                                    gATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGt                                >  1:1740181/1‑100 (MQ=255)
                                                                     aTGAGGGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTc                               >  1:3339830/1‑100 (MQ=255)
                                                                     aTGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTc                               >  2:3398356/1‑100 (MQ=255)
                                                                     aTGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCACCTCTTCTGCTTCCGGTAAGGTc                               <  2:2217578/100‑1 (MQ=255)
                                                                      tGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCggtcaggtca                              <  2:3580362/100‑1 (MQ=255)
                                                                        agagGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACg                            >  1:3062072/1‑100 (MQ=255)
                                                                         gagGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGt                           >  1:1104159/1‑100 (MQ=255)
                                                                         gagGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGt                           >  2:2363736/1‑100 (MQ=255)
                                                                            ggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAGGCGCTGTGACAGCGTGCTTATTTCGGCTCTTCTGCTTCCGGTAAGGTCACGTTCa                        <  2:1635270/100‑1 (MQ=255)
                                                                              aTGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGc                      >  2:332425/1‑100 (MQ=255)
                                                                              aTGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGc                      >  1:1095236/1‑100 (MQ=255)
                                                                               tGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCt                     >  1:927754/1‑100 (MQ=255)
                                                                               tGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCt                     >  2:2632576/1‑100 (MQ=255)
                                                                                gAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTc                    >  1:236196/1‑100 (MQ=255)
                                                                                gAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTc                    >  1:201383/1‑100 (MQ=255)
                                                                                 aaaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTc                    <  2:1768083/99‑1 (MQ=255)
                                                                                  aaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCa                   <  2:1319755/99‑1 (MQ=255)
                                                                                   aaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAg                 <  2:1369733/100‑1 (MQ=255)
                                                                                   aaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAg                 <  2:1095227/100‑1 (MQ=255)
                                                                                    tcttCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGa                >  1:2859760/4‑100 (MQ=39)
                                                                                    aaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGa                <  1:2732205/100‑1 (MQ=255)
                                                                                    aaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGa                <  2:3651713/100‑1 (MQ=255)
                                                                                     cttCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGaa               >  2:2166950/3‑100 (MQ=255)
                                                                                     aCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGaa               <  1:1458529/100‑1 (MQ=255)
                                                                                     aCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGaa               >  2:2801739/1‑100 (MQ=255)
                                                                                      ttCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAAt              >  2:2186335/2‑100 (MQ=255)
                                                                                       tCAAGGCAGAGATAACTCTGCCTTGAAGCTAACTGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATa             <  2:2392084/100‑1 (MQ=255)
                                                                                         aaGGCAGAGTTATGTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa           >  1:2570292/1‑100 (MQ=255)
                                                                                         aaGGCAGAGTTATCTCTGCCTTGAAGATAAATGGGCTTTTACAGCGGGCTTATGTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa           >  2:1968808/1‑100 (MQ=255)
                                                                                         aaGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa           >  2:2589362/1‑100 (MQ=255)
                                                                                         aaGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa           >  1:2998878/1‑100 (MQ=255)
                                                                                         aaGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa           >  1:515043/1‑100 (MQ=255)
                                                                                         aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa           >  1:1734995/1‑100 (MQ=255)
                                                                                         aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa           <  2:149283/100‑1 (MQ=255)
                                                                                          aGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaa          <  1:2164608/100‑1 (MQ=255)
                                                                                           ggCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaaa         <  1:1015194/100‑1 (MQ=255)
                                                                                           ggCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaaa         >  2:1095838/1‑100 (MQ=255)
                                                                                            gCAGAGATAACTCTGCCTTGAAGATATATGCTCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGTTCATGTTCATCTCAAGAATAGaaa         <  2:3659935/99‑1 (MQ=255)
                                                                                               gagaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGAGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaaa         <  2:1900258/96‑1 (MQ=255)
                                                                                                agaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCg    <  2:2587887/100‑1 (MQ=255)
                                                                                                 gaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGc   <  2:3582913/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCC  >  NC_000913/1224141‑1224338

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