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.

TGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTA  >  NC_000913/2765689‑2765883
                                                                                                    |                                                                                               
tGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGc                                                                                                 >  1:405245/1‑100 (MQ=255)
tGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGc                                                                                                 >  1:1447907/1‑100 (MQ=255)
  tCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcct                                                                                               <  1:64411/100‑2 (MQ=255)
    tgctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGGGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccc                                                                                             <  2:752451/100‑1 (MQ=255)
    tgctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccc                                                                                             <  2:2270721/100‑1 (MQ=255)
    tgctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccc                                                                                             >  1:924853/1‑100 (MQ=255)
     gctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccc                                                                                             <  1:75811/99‑1 (MQ=255)
       nggCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccccgg                                                                                          >  2:2776494/2‑98 (MQ=255)
        ggCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccccggt                                                                                         <  1:252548/100‑4 (MQ=255)
        ggCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccccggt                                                                                         >  1:607241/1‑97 (MQ=255)
        ggCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccccggt                                                                                         >  2:193402/1‑97 (MQ=255)
          cGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAg                                                                                      >  2:106523/1‑100 (MQ=255)
          cGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAg                                                                                      <  2:872631/100‑1 (MQ=255)
           gAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACAGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGa                                                                                     <  1:43900/100‑1 (MQ=255)
              taatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGa                                                                                  <  1:1584967/100‑1 (MQ=255)
                  tACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCt                                                                              >  1:704201/1‑100 (MQ=255)
                   aCTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTg                                                                             >  1:3000901/1‑100 (MQ=255)
                   aCTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTg                                                                             >  1:828033/1‑100 (MQ=255)
                    cTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGc                                                                            >  2:3020939/1‑100 (MQ=255)
                    cTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGc                                                                            >  1:2181716/1‑100 (MQ=255)
                     tAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGcc                                                                           >  2:1816967/1‑100 (MQ=255)
                       aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCt                                                                         <  1:2558824/100‑1 (MQ=255)
                             gAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGtt                                                                   >  1:1254353/1‑100 (MQ=255)
                              aGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTg                                                                  <  1:2902237/100‑1 (MQ=255)
                              aGTTGCTGCTTGGTGAGATAAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTg                                                                  >  2:633158/1‑100 (MQ=255)
                               gTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTgg                                                                 <  1:90546/100‑1 (MQ=255)
                               gTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTgg                                                                 <  2:2291231/100‑1 (MQ=255)
                                 tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                               >  1:188947/1‑100 (MQ=255)
                                 tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                               <  2:2461348/100‑1 (MQ=255)
                                 tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                               <  2:364037/100‑1 (MQ=255)
                                 tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                               >  1:297317/1‑100 (MQ=255)
                                  gctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGt                                                              >  2:2801332/1‑100 (MQ=255)
                                    tgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                            >  1:157569/1‑100 (MQ=255)
                                     gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTg                                                             >  2:1684812/1‑100 (MQ=255)
                                     gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                           >  1:2817756/1‑100 (MQ=255)
                                     gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                           >  1:2118438/1‑100 (MQ=255)
                                     gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                           >  2:2986992/1‑100 (MQ=255)
                                     cttccgatctaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                           <  1:138979/90‑1 (MQ=255)
                                       ttGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcccc                                                         <  1:1453679/100‑1 (MQ=255)
                                       ttGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcccc                                                         >  2:227505/1‑100 (MQ=255)
                                        tGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccccc                                                        <  2:1264836/100‑1 (MQ=255)
                                         ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCt                                                       >  1:2505264/1‑100 (MQ=255)
                                         ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCt                                                       >  1:1732734/1‑100 (MQ=255)
                                             agaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTc                                                   >  1:2729651/1‑100 (MQ=255)
                                              gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCa                                                  >  2:2644554/1‑100 (MQ=255)
                                              gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCa                                                  >  1:2704279/1‑100 (MQ=255)
                                               aTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccagatcggaag                                                 >  2:138978/1‑90 (MQ=255)
                                                tCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAgg                                                >  2:2840668/1‑100 (MQ=255)
                                                tCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAgg                                                <  1:581838/100‑1 (MQ=255)
                                                tCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAgg                                                >  1:1533030/1‑100 (MQ=255)
                                                 cAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTaagatcgg                                               >  2:720017/1‑93 (MQ=255)
                                                 cAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGt                                               >  2:405641/1‑100 (MQ=255)
                                                  aGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTg                                              <  2:2546787/100‑1 (MQ=255)
                                                   ggAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTgg                                             >  1:2904039/1‑100 (MQ=255)
                                                    gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGc                                            <  1:426012/100‑1 (MQ=255)
                                                     aaGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCt                                           <  1:2429570/100‑1 (MQ=255)
                                                      agagGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTa                                          <  2:2209941/100‑1 (MQ=255)
                                                       gagGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTaa                                         <  1:428753/100‑1 (MQ=255)
                                                       gagGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTaa                                         >  2:1261636/1‑100 (MQ=255)
                                                           ttACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAAttt                                     >  1:2136183/1‑100 (MQ=255)
                                                            tACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTa                                    <  2:828026/100‑1 (MQ=255)
                                                             aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                   >  2:2479300/1‑100 (MQ=255)
                                                             aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                   <  1:1101029/100‑1 (MQ=255)
                                                             aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                   <  1:1940851/100‑1 (MQ=255)
                                                              ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                  >  1:1133962/1‑100 (MQ=255)
                                                              ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                  <  2:471441/100‑1 (MQ=255)
                                                              ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                  <  1:117837/100‑1 (MQ=255)
                                                              ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                  >  1:1619757/1‑100 (MQ=255)
                                                              ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                  <  1:2607384/100‑1 (MQ=255)
                                                              ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                  <  1:297883/100‑1 (MQ=255)
                                                              ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                  <  1:942483/100‑1 (MQ=255)
                                                               cGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                 >  2:2380375/1‑100 (MQ=255)
                                                                       tGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                         <  1:1839824/100‑1 (MQ=255)
                                                                       tGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                         <  2:2655961/100‑1 (MQ=255)
                                                                        ggACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                        <  1:1874854/100‑1 (MQ=255)
                                                                        ggACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                        <  1:2612277/100‑1 (MQ=255)
                                                                        ggACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                        <  2:2517079/100‑1 (MQ=255)
                                                                         gACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                       >  1:2700790/1‑100 (MQ=255)
                                                                          aCATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                        <  2:82398/98‑1 (MQ=255)
                                                                           catcatCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                     <  1:1424961/100‑1 (MQ=255)
                                                                            atcatcCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACtt                    <  2:2101519/100‑1 (MQ=255)
                                                                              catcCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCa                  >  2:1550702/1‑100 (MQ=255)
                                                                                  ccccATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg              >  1:872522/1‑100 (MQ=255)
                                                                                   cccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg               <  1:1684826/100‑1 (MQ=255)
                                                                                   cccATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggggg             >  2:1326505/1‑100 (MQ=255)
                                                                                   cccATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggggg             >  2:1519600/1‑100 (MQ=255)
                                                                                    ccATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg            <  1:2867559/100‑1 (MQ=255)
                                                                                    ccATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg            >  2:778961/1‑100 (MQ=255)
                                                                                      aTTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta          >  1:2766298/1‑100 (MQ=255)
                                                                                       ttttGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtat         <  2:570523/100‑1 (MQ=255)
                                                                                          tgtgGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtattat      <  2:188944/100‑1 (MQ=255)
                                                                                           gtgGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTattatt     >  2:1560789/1‑100 (MQ=255)
                                                                                            tgGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTg    <  1:1858260/100‑1 (MQ=255)
                                                                                            tgGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTg    >  2:3013783/1‑100 (MQ=255)
                                                                                            tgGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTg    <  1:2644577/100‑1 (MQ=255)
                                                                                             ggACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGt   >  2:2616269/1‑100 (MQ=255)
                                                                                             ggACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGt   <  1:896705/100‑1 (MQ=255)
                                                                                             ggACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGt   >  2:792566/1‑100 (MQ=255)
                                                                                              gACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTa  >  1:2088082/1‑100 (MQ=255)
                                                                                              gACTGCCCCCCAGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGCGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGAATTATTGTa  >  2:1171014/1‑100 (MQ=255)
                                                                                                    |                                                                                               
TGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTA  >  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: