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.C84.1% 316.5 / 62.1 113intergenic (‑11/‑131)yfjM/rnlACP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; RNase LS
Reads supporting (aligned to +/- strand):  ref base . (11/7);  new base C (43/50);  total (55/58)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

TTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAG  >  NC_000913/2765688‑2765886
                                                                                                    |                                                                                                        
ttGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTg                                                                                                           >  2:2776846/1‑100 (MQ=255)
ttGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTg                                                                                                           >  2:3400580/1‑100 (MQ=255)
        tgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccccccc                                                                                             <  1:3267606/100‑2 (MQ=255)
         ggCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCgg                                                                                            >  1:250455/1‑100 (MQ=255)
          gCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccccccc                                                                                             >  1:3785298/1‑97 (MQ=255)
          gCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGt                                                                                           >  2:2820894/1‑100 (MQ=255)
           cGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccccggt                                                                                          <  1:83857/100‑5 (MQ=255)
           cGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccccgg                                                                                           >  2:2117219/1‑96 (MQ=255)
            gAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccccggta                                                                                         <  2:309831/100‑6 (MQ=255)
              ataatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGAc                                                                                       <  1:1248549/100‑1 (MQ=255)
              ataatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGAGGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGa                                                                                        <  1:848340/100‑1 (MQ=255)
               taatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGAc                                                                                       <  1:1619856/100‑1 (MQ=255)
               taatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGAc                                                                                       <  1:744011/100‑1 (MQ=255)
               taatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGAc                                                                                       <  2:191545/100‑1 (MQ=255)
               taatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACg                                                                                      <  2:670601/100‑1 (MQ=255)
                aatTACTAGTTATCGTGTTGCTTCTTGGTGAGATCAGGAAAAGGCTACCGAAACTATGCACATCCTCCCCCATTTTGTTGACTG‑‑CCCCCCTCGGTAGACg                                                                                      <  1:552224/100‑1 (MQ=255)
                aatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACg                                                                                      >  2:3126252/1‑100 (MQ=255)
                 atTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccccggcagac                                                                                      >  2:700146/1‑90 (MQ=255)
                  ttACTAAGTATGGAGTGGCGGCTTGGTGAGATCGGGAGGGGGTTAGCGAGACGATGGACATCATTCCGCATTTTGTGGACTGTCCCCCCCGGTAGACg                                                                                      <  1:334614/98‑1 (MQ=255)
                  ttACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGAt                                                                                    <  2:865159/100‑1 (MQ=255)
                  ttACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGAt                                                                                    <  2:2927383/100‑1 (MQ=255)
                  ttACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGAt                                                                                    <  2:1699957/100‑1 (MQ=255)
                       aaGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGc                                                                               <  1:1168369/100‑1 (MQ=255)
                         gTATCGAGTTGCTGCTTGGTGAGATCAGGCAGAGGTTACCGAAACGATGGACATCATCCCCCATTTGGTGGACTGCCCCCCCCGGTAGACGATCCTGccc                                                                             <  1:2280749/100‑1 (MQ=255)
                         gTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGccc                                                                             <  2:2009016/100‑1 (MQ=255)
                         gTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGAATG‑‑CCCCTCCCGGTCGACGGTCCTGccc                                                                             >  1:1323216/1‑100 (MQ=255)
                         gTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTATGTGGACTGCCCCCCCCGGTAGACGATCCTGccc                                                                             >  1:1997858/1‑100 (MQ=255)
                          tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGGTAGACGATCCTGccc                                                                             <  2:1395830/100‑1 (MQ=255)
                          tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCt                                                                            >  1:3062425/1‑100 (MQ=255)
                          tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGAGCCTGCCCt                                                                            >  2:2101989/1‑100 (MQ=255)
                           aTCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCta                                                                           <  1:3701372/100‑1 (MQ=255)
                            tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGACCCTGcccg                                                                            >  2:689140/1‑97 (MQ=255)
                                   gctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGa                                                                   >  2:128996/1‑100 (MQ=255)
                                   gctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                                  >  2:87360/1‑100 (MQ=255)
                                    ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                                  <  2:2716613/100‑1 (MQ=255)
                                    ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                                  <  2:2062467/100‑1 (MQ=255)
                                    ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                                  <  2:1652391/100‑1 (MQ=255)
                                     tgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTg                                                                <  1:3687523/100‑1 (MQ=255)
                                       ctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                               >  2:3119826/1‑100 (MQ=255)
                                       ctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCCGCCCTACAGTTGGAGTGc                                                               >  1:2812830/1‑100 (MQ=255)
                                       ctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                              >  2:2724280/1‑100 (MQ=255)
                                        ttGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccc                                                             <  1:3192229/100‑1 (MQ=255)
                                           gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccccc                                                           <  1:384075/100‑1 (MQ=255)
                                           gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccccc                                                           <  1:2328393/100‑1 (MQ=255)
                                           gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccccc                                                           >  1:1370695/1‑100 (MQ=255)
                                           gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCt                                                          <  2:2791476/100‑1 (MQ=255)
                                            tGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTCGTGGACTGCCCCCCCCGGTAGGCGAACCCGCCCTATAGTTGGAGTGcccac                                                           >  2:433875/1‑97 (MQ=255)
                                             gagaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTa                                                         >  1:2442247/1‑100 (MQ=255)
                                                aTCAGGAAGAGGTTACCGAAACTATGGACATCATCCCCCCTTTTGCGTCCTGCCCCCCCCCCTAGCCAACCCCGCCCTACAGTTGGAgcaccc                                                             >  1:115968/1‑88 (MQ=255)
                                                aTCAGGAAGAGGTTACCGAAACGATGGACATCATCACCAATTTTGAGGACTG‑‑CCCCAACCGGGAGACGAACCTGCAAAAAAGTCGAAGTGCCCCCTAGTc                                                      >  1:115950/1‑100 (MQ=255)
                                                 tCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAg                                                    >  2:3770763/1‑100 (MQ=255)
                                                   aGGAAAGGGTTACCGAACCGATGGACATTATCCCCCATTTTTTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAgg                                                   <  1:959967/100‑1 (MQ=255)
                                                    ggAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGt                                                  <  2:2051053/100‑1 (MQ=255)
                                                     gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTg                                                 >  2:3132059/1‑100 (MQ=255)
                                                     gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTg                                                 <  2:3301406/100‑1 (MQ=255)
                                                     gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTg                                                 <  2:1691089/100‑1 (MQ=255)
                                                     gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTg                                                 <  1:1872050/100‑1 (MQ=255)
                                                      aaGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTgg                                                <  2:2609298/100‑1 (MQ=255)
                                                           gTTACCGAAACGATGGACATCATCCCCCATTTTGCGGACAGCCCCCCCCGGTTGACGCTCCTGCCCTATAGTTCGTGTGCCCCCTATTCAGGTGTCTaaa                                           >  1:239266/1‑100 (MQ=255)
                                                            ttACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAAt                                          <  1:3474772/100‑1 (MQ=255)
                                                             tagcGAAACGATGGACATGATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAAtt                                         <  1:1717902/97‑1 (MQ=255)
                                                             tACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAAtt                                         >  2:3703550/1‑100 (MQ=255)
                                                              aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCTCTATAGTTGGAGTGCCCCCTCTTCAGGTGGCTAAAttt                                        >  1:3045995/1‑100 (MQ=255)
                                                              aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAAttt                                        >  2:3344244/1‑100 (MQ=255)
                                                              aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAAttt                                        >  2:3211942/1‑100 (MQ=255)
                                                               ccGAAACGATGGACATCATCCCCCATTTTGTGGCCTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTa                                       >  1:2279034/1‑100 (MQ=255)
                                                               ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTa                                       <  1:1403010/100‑1 (MQ=255)
                                                               ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAATGCCCCCTATTCAGGTGGCTAAAttt                                        >  2:2238698/1‑99 (MQ=255)
                                                               ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCCGCGCTATACATGGAGTGCCCCCTATTCAGGTGGCTAAATTTa                                       >  1:2667051/1‑100 (MQ=255)
                                                                cGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                      >  1:1268679/1‑100 (MQ=255)
                                                                cGAAACGATGGACATCATCCCCCATTTCGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTACAGTTGGAGTGCCCCCTATTCAGGTGGCTAAGTTTa                                       >  2:472099/1‑99 (MQ=255)
                                                                 gAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                     <  2:3744866/100‑1 (MQ=255)
                                                                  aaaCGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                    <  2:485133/100‑1 (MQ=255)
                                                                  aaaCGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                    <  2:882693/100‑1 (MQ=255)
                                                                      gATGGACATCATCCCCCATTTTGTGAACTG‑‑CACACCCGGGTAGACGATCCTCCCCTATAGTCGGAGTGCCCCCTATTCAGGTTGCGAAATTTAGCa                                    >  1:942017/1‑96 (MQ=255)
                                                                       aTGGACATCCTCCCCCATTTTGTGGACTGCCCCCCCCGGAAGACGCCCCTGCCCTGTAGTTGGAGTGCCCCCTCTCCAGGTGGCGAAATTTAGCAAATc                                >  1:471882/1‑99 (MQ=255)
                                                                       aTGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCa                               <  2:359486/100‑1 (MQ=255)
                                                                        tGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCaa                              >  2:525160/1‑100 (MQ=255)
                                                                        tGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAAt                             <  1:1051531/100‑1 (MQ=255)
                                                                         ggACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAAt                             >  1:2848418/1‑100 (MQ=255)
                                                                          gACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                            >  2:660161/1‑100 (MQ=255)
                                                                           aCATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                           >  2:1530848/1‑100 (MQ=255)
                                                                            catcatCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                          >  2:787636/1‑100 (MQ=255)
                                                                            catcatCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATacac                         >  2:2123880/1‑100 (MQ=255)
                                                                                  cccccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCTCTATAGTTGGCGTGCCCCCTATTCACGTGTCTAAATTTAACAAATCAATACACtt                       >  1:3596285/1‑97 (MQ=255)
                                                                                  cccccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCa                     <  1:3195251/99‑1 (MQ=255)
                                                                                   caccATTTTGTGGACTG‑‑CCACCCCCGGTAGACGATCCTGCCCTGTAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTCCGAAAATAAATACACTTCAgg                   >  1:949528/3‑100 (MQ=255)
                                                                                    cccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                  >  1:754678/1‑100 (MQ=255)
                                                                                    cccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATAAATACACTTCAggg                  >  1:2622517/1‑100 (MQ=255)
                                                                                     ccATTTTGTGGACTGCCCCCCCCGCTAGACGATCCTGCCCTATAGTTTGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                  <  1:1327626/99‑1 (MQ=255)
                                                                                       aTTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg               <  2:1013531/100‑1 (MQ=255)
                                                                                       aTTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg               <  2:1711068/100‑1 (MQ=255)
                                                                                         tttGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta             <  2:1483453/100‑1 (MQ=255)
                                                                                           tgtgGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtatt           >  1:517698/1‑100 (MQ=255)
                                                                                           tgtgGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtatt           <  1:2151439/100‑1 (MQ=255)
                                                                                             tgGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTTGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtattat         <  1:1861805/100‑1 (MQ=255)
                                                                                             tgGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtattat         <  1:1816716/100‑1 (MQ=255)
                                                                                              ggACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTattatt        >  1:211456/1‑100 (MQ=255)
                                                                                               gACTGCCCCCCCCGGTAGACGTTCCTGCCCCATAGCTGGAGCGCCCCCTGTTCAGCTGGCTAAATTTCGCAAATCAATCCACTCCAgggcgcgattct         >  1:3523987/1‑89 (MQ=255)
                                                                                               gACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTg       >  2:1938811/1‑100 (MQ=255)
                                                                                                    |                                                                                                        
TTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAG  >  NC_000913/2765688‑2765886

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: