Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,765,794:1 (C)7→8 100% intergenic (‑18/‑124) yfjM ← / → rnlA CP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; RNase LS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,765,7871.C78.6% 377.7 / 114.1 145intergenic (‑11/‑131)yfjM/rnlACP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; RNase LS
Reads supporting (aligned to +/- strand):  ref base . (18/13);  new base C (52/61);  total (70/75)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

TCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGA  >  NC_000913/2765691‑2765885
                                                                                                 |                                                                                                        
  tgctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccc                                                                                                   <  2:3532161/99‑1 (MQ=255)
     tgGCGAATAATTACTAAGTATCGAGTTGCTTCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATCGACATCATCCCCCCTCTTGTGGACTG‑‑cccccccc                                                                                            >  2:380316/1‑99 (MQ=255)
     tgGCGAATAATTACTAAGTATCGAGTTGATGCTTGGTGAGATAAGGAAGAGATTACCGAAAAGATGGACACCATCGACCAGTTTGTGGACTg                                                                                                           >  2:2940700/1‑92 (MQ=255)
      ggCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccac                                                                                             >  2:3130936/1‑97 (MQ=255)
       gCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccccg                                                                                           <  2:2800607/99‑3 (MQ=255)
          aataatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccccccccggt                                                                                        <  1:573792/99‑6 (MQ=255)
            taatTACTAAGTATCGAGTTGCTGCTTTGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCTTCCCCCATTTTATGGACTGCCCCCCCCGGTAg                                                                                        <  2:3773784/98‑1 (MQ=255)
            taatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGAc                                                                                      >  2:2549257/1‑100 (MQ=255)
              atTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGa                                                                                    <  1:123977/100‑1 (MQ=255)
               ttACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTCTGGACTGCCCCCCCCGGTAAACGa                                                                                    >  2:1317134/1‑99 (MQ=255)
                tACTAAGTATCGAGTTGCTGCTTGGTTATATTAGTAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATc                                                                                  <  1:349216/100‑1 (MQ=255)
                tACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCCTTTTGTGGACTGCCCCCCCCGGTAGACGAt                                                                                   >  2:26488/1‑99 (MQ=255)
                tACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATcc                                                                                 <  2:2823236/100‑1 (MQ=255)
                 tccgatcttCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATcc                                                                                 <  1:99881/92‑1 (MQ=255)
                 aCTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATcc                                                                                 <  1:2126295/100‑1 (MQ=255)
                  cTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCt                                                                                >  2:717284/1‑100 (MQ=255)
                    aaGTATCGTGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACTGAAACGATGGAAATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGc                                                                              <  1:1783689/100‑1 (MQ=255)
                    aaGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGc                                                                              <  2:3294771/100‑1 (MQ=255)
                    aaGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCAATTTGTGGACTGCCCCCCCCGGTAGACGATCCTGc                                                                              >  2:2058966/1‑100 (MQ=255)
                    aaGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATAAACAATTTTGTGGACTGCCCCCCCCGGTAGAAGATCCTGc                                                                              >  1:3044131/1‑100 (MQ=255)
                      gTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTAATGGACTC‑‑CCCCCACCGATAGACGATCCTGcc                                                                             >  1:3553009/1‑99 (MQ=255)
                        aTCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCta                                                                          <  1:1317148/100‑1 (MQ=255)
                         tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCtat                                                                         <  1:2789113/100‑1 (MQ=255)
                            aGTTTCTGCTTAGTGAGATCAGGAAGTGGTTACCGAAACTATTGACATCATCCCCCATTTTATATACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGt                                                                      <  2:1492248/100‑1 (MQ=255)
                             gTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGtt                                                                     <  1:2296670/100‑1 (MQ=255)
                             gTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCACGGTAGACGATCCTGCCCTATAGt                                                                      >  1:3783381/1‑99 (MQ=255)
                              ttGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTg                                                                    <  2:82175/100‑1 (MQ=255)
                               tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTgg                                                                   <  1:2348898/100‑1 (MQ=255)
                                gctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGa                                                                  >  2:846102/1‑100 (MQ=255)
                                 ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                                 <  2:2453684/100‑1 (MQ=255)
                                  tgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGt                                                                <  2:2085248/100‑1 (MQ=255)
                                  tgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTg                                                               <  2:3127447/100‑1 (MQ=255)
                                  tgctTGGTGAGACCAGGAAGAGGGTACCGAAACGATGGACATCATCTCCCATTTTATGGAATA‑‑CCCCGGCGGAAGATGATCATACCATATAGTTGCAGTg                                                               >  2:1361991/1‑100 (MQ=255)
                                      tGGTGAGATCAGTAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccc                                                            <  2:1722268/100‑1 (MQ=255)
                                      tGGTGAGATCAGGAAGAGGTTACCGAAACTATGGACATCATCCTCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccc                                                            <  1:1764058/100‑1 (MQ=255)
                                       ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCACCCCGGAAGAAGATACTGCCCTATAGTTGGAGTGccc                                                            >  1:3171328/1‑99 (MQ=255)
                                        gTGAGATCAGGAAGAGGTTACCGAAACTATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccccc                                                          <  1:1496617/100‑1 (MQ=255)
                                        gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccccc                                                          >  2:1549752/1‑100 (MQ=255)
                                        gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccccc                                                          <  1:123418/100‑1 (MQ=255)
                                        gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccccc                                                          <  1:1852091/100‑1 (MQ=255)
                                         tGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCGCATTTTGTGGACTG‑‑CCCCCGCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCt                                                         >  2:2122984/1‑100 (MQ=255)
                                         tGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCt                                                         <  1:3373821/100‑1 (MQ=255)
                                           agaTCAGGAAGAGGTTACCGAAACGATGGACATCATACCCCATTTTGTAGACTGCCCCCCCCGGTAAACTATCCTGCCCTATAGGTGGAGTGCCCCCTAt                                                       >  2:1404348/1‑100 (MQ=255)
                                            gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTAtt                                                      >  2:3112178/1‑100 (MQ=255)
                                             aTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACAATCCTGCACTATAGATGGAGTGCCCCCTATTc                                                     >  1:596220/1‑100 (MQ=255)
                                              ttaGGTAGAGGTTACTGAAATGATTGACATCATCTCCTATTTTGTGGAATG‑‑CCCCTCACGGTAGACGATCCTGCCCGATAGTTGGAGTGCCCCCTATTc                                                     <  1:758677/97‑1 (MQ=255)
                                              tCAGGAAGAGGTTACCGAAACGATGGACATCATTCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCa                                                    <  1:45466/100‑1 (MQ=255)
                                              tCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCa                                                    <  2:3103539/100‑1 (MQ=255)
                                              tCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCa                                                    <  1:2629083/100‑1 (MQ=255)
                                               cAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAg                                                   >  1:2692814/1‑100 (MQ=255)
                                                aGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAgg                                                  <  2:733234/100‑1 (MQ=255)
                                                aGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAgg                                                  >  2:1907184/1‑100 (MQ=255)
                                                 ggAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTg                                                <  2:448067/100‑1 (MQ=255)
                                                  gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGCCAGACGATCCTGCCCTATCGCTGGAGCGCCCCCTACTCAGGt                                                 >  2:2030934/1‑99 (MQ=255)
                                                   aaGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTgg                                               <  2:1979561/100‑1 (MQ=255)
                                                    agagGTTACCGTAATGATGAACATAATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTg                                                <  1:3572353/98‑1 (MQ=255)
                                                    agagGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATCGTTGGAGTGCCCCCTATTCCGGTGGc                                              >  2:3492355/1‑100 (MQ=255)
                                                    agagGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGc                                              <  2:2023381/100‑1 (MQ=255)
                                                          tACCGAAACGATGGACATCATCCCCCATTTTGTGTACTT‑‑CCCCCCCGTTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAAttt                                       <  1:1111684/100‑1 (MQ=255)
                                                            ctgAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAATTGGAAAGCACACTATTCAGGTGGCTAAATTTa                                      >  1:746249/3‑100 (MQ=255)
                                                            ccGAAACGGTGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTTGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                     <  2:998640/100‑1 (MQ=255)
                                                            ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTa                                      >  1:3613633/1‑100 (MQ=255)
                                                             cGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTACCACCTATTCAGGTGGCTAAATTTa                                      >  1:1373174/1‑99 (MQ=255)
                                                             cGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                     <  2:3278889/99‑1 (MQ=255)
                                                             cGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                    >  2:251564/1‑100 (MQ=255)
                                                              naaaCGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                    >  2:2703625/2‑100 (MQ=255)
                                                              gAAACGATGGACCTCATCCCCCATTTTGTGGACTGTCCCCCCCAGTATACGATGGCGACCTATAATTGGAGTGCCGCCTATTCAGATGGCTAAATTTAGc                                    <  2:36878/100‑1 (MQ=255)
                                                               aaaCGATGTACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                   <  1:2453080/100‑1 (MQ=255)
                                                               aaaCGATGGCCATCATCCGTCATTTTGTGGACTG‑‑CCACCACCGGTTGACGGTCCTGCGCAATGGTTGGAGTGCCACATATCAAGGTGGTTaatgtttgct                                   <  2:2277401/100‑9 (MQ=255)
                                                                aaCGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                   <  2:44689/99‑1 (MQ=255)
                                                                   gATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATc                               <  1:1066165/100‑1 (MQ=255)
                                                                     tGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCaa                             <  2:2447236/100‑1 (MQ=255)
                                                                      ggACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAAt                            <  1:377979/100‑1 (MQ=255)
                                                                       gACATCATCCCCCATTTTGTGGACTG‑‑CCCCCTCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                           <  2:2634025/100‑1 (MQ=255)
                                                                        aCATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                          <  1:908340/100‑1 (MQ=255)
                                                                        aCATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGATTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                         >  1:3638063/1‑100 (MQ=255)
                                                                            catcCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTAGCTAAATTTAACAAATAAATacac                        >  1:3793409/1‑98 (MQ=255)
                                                                               cccccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                   >  2:3855451/1‑100 (MQ=255)
                                                                               cccccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                   >  2:2162883/1‑100 (MQ=255)
                                                                                ccccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgg                  >  1:2295103/1‑100 (MQ=255)
                                                                                ccccATTTTGTGCACTGCCCCCCCCCGGTATACGATCCTGCCCTATATTTGGATTTCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                   <  1:3330322/100‑1 (MQ=255)
                                                                                ccccATGTTGTGGACTGCCCCCCCCGGCAGACGATCCTGCCCTATAATTGGATCGCCGCCGATTCCGGTGGCTAAATTCCGCAAATCAATACCATTCatg                  >  1:1942477/1‑98 (MQ=255)
                                                                                 cccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                 <  1:535590/100‑1 (MQ=255)
                                                                                 cccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTAAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                 >  1:2033373/1‑100 (MQ=255)
                                                                                 cccATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                >  1:3292940/1‑100 (MQ=255)
                                                                                 cccATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                >  1:2596708/1‑100 (MQ=255)
                                                                                  ccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                <  2:1123924/100‑1 (MQ=255)
                                                                                   cATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggggg               >  1:2522406/1‑100 (MQ=255)
                                                                                   cATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggggg               <  1:883267/100‑1 (MQ=255)
                                                                                   cATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                <  2:1829432/99‑1 (MQ=255)
                                                                                   cATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                <  2:1341216/99‑1 (MQ=255)
                                                                                      tttGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATAAATACACTTCAGGGGGGta            >  2:3463027/1‑100 (MQ=255)
                                                                                         gtgGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtatta         >  1:3195971/1‑100 (MQ=255)
                                                                                         gtgGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCATATTAAGGTGGATAAATTTAGCAAATCAATACACTTCAGGGGGGtatta         >  1:183434/1‑100 (MQ=255)
                                                                                          tgGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtattat        >  1:3509659/1‑100 (MQ=255)
                                                                                           ggACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTattatt       >  2:694707/1‑100 (MQ=255)
                                                                                           ggACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTattatt       >  2:1163039/1‑100 (MQ=255)
                                                                                           ggACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAAGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTattatt       >  1:3277453/1‑100 (MQ=255)
                                                                                            gACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTg      >  2:36985/1‑100 (MQ=255)
                                                                                                g‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGATGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTaga  >  2:1820354/1‑99 (MQ=255)
                                                                                                 |                                                                                                        
TCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGA  >  NC_000913/2765691‑2765885

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