Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,765,794 (C)7→6 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,7880C.97.4% 494.5 / 4.7 113intergenic (‑12/‑130)yfjM/rnlACP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; RNase LS
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (57/53);  total (58/55)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

TGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTA  >  NC_000913/2765689‑2765883
                                                                                                    |                                                                                               
tGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑c                                                                                                 >  1:405245/1‑100 (MQ=255)
tGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑c                                                                                                 >  1:1447907/1‑100 (MQ=255)
  tCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑cct                                                                                               <  1:64411/100‑2 (MQ=255)
    tgctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGGGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑ccccc                                                                                             <  2:752451/100‑1 (MQ=255)
    tgctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑ccccc                                                                                             <  2:2270721/100‑1 (MQ=255)
    tgctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑ccccc                                                                                             >  1:924853/1‑100 (MQ=255)
     gctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑ccccc                                                                                             <  1:75811/99‑1 (MQ=255)
       nggCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑ccccccgg                                                                                          >  2:2776494/2‑98 (MQ=255)
        ggCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑ccccccggt                                                                                         <  1:252548/100‑4 (MQ=255)
        ggCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑ccccccggt                                                                                         >  1:607241/1‑97 (MQ=255)
        ggCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑ccccccggt                                                                                         >  2:193402/1‑97 (MQ=255)
          cGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAg                                                                                      >  2:106523/1‑100 (MQ=255)
          cGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAg                                                                                      <  2:872631/100‑1 (MQ=255)
           gAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACAGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGa                                                                                     <  1:43900/100‑1 (MQ=255)
              taatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGa                                                                                  <  1:1584967/100‑1 (MQ=255)
                  tACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCt                                                                              >  1:704201/1‑100 (MQ=255)
                   aCTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTg                                                                             >  1:3000901/1‑100 (MQ=255)
                   aCTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTg                                                                             >  1:828033/1‑100 (MQ=255)
                    cTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGc                                                                            >  2:3020939/1‑100 (MQ=255)
                    cTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGc                                                                            >  1:2181716/1‑100 (MQ=255)
                     tAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGcc                                                                           >  2:1816967/1‑100 (MQ=255)
                       aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCt                                                                         <  1:2558824/100‑1 (MQ=255)
                             gAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGtt                                                                   >  1:1254353/1‑100 (MQ=255)
                              aGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTg                                                                  <  1:2902237/100‑1 (MQ=255)
                              aGTTGCTGCTTGGTGAGATAAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTg                                                                  >  2:633158/1‑100 (MQ=255)
                               gTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTgg                                                                 <  1:90546/100‑1 (MQ=255)
                               gTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTgg                                                                 <  2:2291231/100‑1 (MQ=255)
                                 tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                               >  1:188947/1‑100 (MQ=255)
                                 tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                               <  2:2461348/100‑1 (MQ=255)
                                 tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                               <  2:364037/100‑1 (MQ=255)
                                 tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                               >  1:297317/1‑100 (MQ=255)
                                  gctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGt                                                              >  2:2801332/1‑100 (MQ=255)
                                    tgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                            >  1:157569/1‑100 (MQ=255)
                                     gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTg                                                             >  2:1684812/1‑100 (MQ=255)
                                     gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                           >  1:2817756/1‑100 (MQ=255)
                                     gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                           >  1:2118438/1‑100 (MQ=255)
                                     gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                           >  2:2986992/1‑100 (MQ=255)
                                     cttccgatctaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                           <  1:138979/90‑1 (MQ=255)
                                       ttGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcccc                                                         <  1:1453679/100‑1 (MQ=255)
                                       ttGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcccc                                                         >  2:227505/1‑100 (MQ=255)
                                        tGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccccc                                                        <  2:1264836/100‑1 (MQ=255)
                                         ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCt                                                       >  1:2505264/1‑100 (MQ=255)
                                         ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCt                                                       >  1:1732734/1‑100 (MQ=255)
                                             agaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTc                                                   >  1:2729651/1‑100 (MQ=255)
                                              gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCa                                                  >  2:2644554/1‑100 (MQ=255)
                                              gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCa                                                  >  1:2704279/1‑100 (MQ=255)
                                               aTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccagatcggaag                                                 >  2:138978/1‑90 (MQ=255)
                                                tCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAgg                                                >  2:2840668/1‑100 (MQ=255)
                                                tCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAgg                                                <  1:581838/100‑1 (MQ=255)
                                                tCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAgg                                                >  1:1533030/1‑100 (MQ=255)
                                                 cAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTaagatcgg                                               >  2:720017/1‑93 (MQ=255)
                                                 cAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGt                                               >  2:405641/1‑100 (MQ=255)
                                                  aGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTg                                              <  2:2546787/100‑1 (MQ=255)
                                                   ggAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTgg                                             >  1:2904039/1‑100 (MQ=255)
                                                    gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGc                                            <  1:426012/100‑1 (MQ=255)
                                                     aaGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCt                                           <  1:2429570/100‑1 (MQ=255)
                                                      agagGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTa                                          <  2:2209941/100‑1 (MQ=255)
                                                       gagGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTaa                                         <  1:428753/100‑1 (MQ=255)
                                                       gagGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTaa                                         >  2:1261636/1‑100 (MQ=255)
                                                           ttACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAAttt                                     >  1:2136183/1‑100 (MQ=255)
                                                            tACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTa                                    <  2:828026/100‑1 (MQ=255)
                                                             aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                   >  2:2479300/1‑100 (MQ=255)
                                                             aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                   <  1:1101029/100‑1 (MQ=255)
                                                             aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                   <  1:1940851/100‑1 (MQ=255)
                                                              ccGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                  >  1:1133962/1‑100 (MQ=255)
                                                              ccGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                  <  2:471441/100‑1 (MQ=255)
                                                              ccGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                  <  1:117837/100‑1 (MQ=255)
                                                              ccGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                  >  1:1619757/1‑100 (MQ=255)
                                                              ccGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                  <  1:2607384/100‑1 (MQ=255)
                                                              ccGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                  <  1:297883/100‑1 (MQ=255)
                                                              ccGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                  <  1:942483/100‑1 (MQ=255)
                                                               cGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                 >  2:2380375/1‑100 (MQ=255)
                                                                       tGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                         <  1:1839824/100‑1 (MQ=255)
                                                                       tGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                         <  2:2655961/100‑1 (MQ=255)
                                                                        ggACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                        <  1:1874854/100‑1 (MQ=255)
                                                                        ggACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                        <  1:2612277/100‑1 (MQ=255)
                                                                        ggACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                        <  2:2517079/100‑1 (MQ=255)
                                                                         gACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                       >  1:2700790/1‑100 (MQ=255)
                                                                          aCATCATCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                        <  2:82398/98‑1 (MQ=255)
                                                                           catcatCCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                     <  1:1424961/100‑1 (MQ=255)
                                                                            atcatcCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACtt                    <  2:2101519/100‑1 (MQ=255)
                                                                              catcCCCCATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCa                  >  2:1550702/1‑100 (MQ=255)
                                                                                  ccccATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg              >  1:872522/1‑100 (MQ=255)
                                                                                   cccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg               <  1:1684826/100‑1 (MQ=255)
                                                                                   cccATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggggg             >  2:1326505/1‑100 (MQ=255)
                                                                                   cccATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggggg             >  2:1519600/1‑100 (MQ=255)
                                                                                    ccATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg            <  1:2867559/100‑1 (MQ=255)
                                                                                    ccATTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg            >  2:778961/1‑100 (MQ=255)
                                                                                      aTTTTGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta          >  1:2766298/1‑100 (MQ=255)
                                                                                       ttttGTGGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtat         <  2:570523/100‑1 (MQ=255)
                                                                                          tgtgGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtattat      <  2:188944/100‑1 (MQ=255)
                                                                                           gtgGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTattatt     >  2:1560789/1‑100 (MQ=255)
                                                                                            tgGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTg    <  1:1858260/100‑1 (MQ=255)
                                                                                            tgGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTg    >  2:3013783/1‑100 (MQ=255)
                                                                                            tgGACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTg    <  1:2644577/100‑1 (MQ=255)
                                                                                             ggACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGt   >  2:2616269/1‑100 (MQ=255)
                                                                                             ggACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGt   <  1:896705/100‑1 (MQ=255)
                                                                                             ggACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGt   >  2:792566/1‑100 (MQ=255)
                                                                                              gACTG‑‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTa  >  1:2088082/1‑100 (MQ=255)
                                                                                              gACTG‑‑CCCCCCAGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGCGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGAATTATTGTa  >  2:1171014/1‑100 (MQ=255)
                                                                                                    |                                                                                               
TGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTA  >  NC_000913/2765689‑2765883

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