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.C87.1% 275.7 / 36.5 93intergenic (‑11/‑131)yfjM/rnlACP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; RNase LS
Reads supporting (aligned to +/- strand):  ref base . (8/4);  new base C (36/43);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

TTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATT  >  NC_000913/2765688‑2765880
                                                                                                     |                                                                                                 
ttGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTg                                                                                                    <  2:1028206/100‑1 (MQ=255)
ttGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTg                                                                                                    <  1:332756/100‑1 (MQ=255)
 tGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTg                                                                                                    <  1:2846418/99‑1 (MQ=255)
 tGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑c                                                                                                 <  2:1950256/100‑1 (MQ=255)
 tGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑c                                                                                                 >  2:2683485/1‑100 (MQ=255)
 tGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑c                                                                                                 <  2:1493482/100‑1 (MQ=255)
 tGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑c                                                                                                 >  2:664965/1‑100 (MQ=255)
 tGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCACCCATTTTGTGGACTg                                                                                                    >  2:1291986/1‑99 (MQ=255)
  gTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cc                                                                                               >  2:2293542/1‑100 (MQ=255)
   tCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccc                                                                                              <  2:803182/100‑1 (MQ=255)
    ctgctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACCTCATCCCCCATTCTGTGGACTG‑‑ccc                                                                                              >  1:1903767/1‑99 (MQ=255)
     tgctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccc                                                                                          <  2:243200/100‑1 (MQ=255)
      gctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccccc                                                                                         <  1:829608/100‑1 (MQ=255)
      gctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccccc                                                                                         <  1:1335638/100‑1 (MQ=255)
        tgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccccccc                                                                                       <  2:424247/100‑2 (MQ=255)
           cGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGCTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccccggt                                                                                    <  2:1930585/100‑5 (MQ=255)
            gAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccccggta                                                                                   <  2:1217608/100‑6 (MQ=255)
            gAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccccggta                                                                                   <  1:1771385/100‑6 (MQ=255)
            gAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccccggt                                                                                    <  2:1922021/99‑5 (MQ=255)
             aataatTACTAAATATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAg                                                                                   <  1:2474666/100‑1 (MQ=255)
              ataatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGa                                                                                  <  1:2346787/100‑1 (MQ=255)
               taattaatAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGAc                                                                                 <  2:2428226/100‑1 (MQ=255)
               taatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGAc                                                                                 <  1:2560590/100‑1 (MQ=255)
               taatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGAc                                                                                 <  2:71566/100‑1 (MQ=255)
                aatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccccccccggtagac                                                                               <  2:2805868/100‑10 (MQ=255)
                 atTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCTCCATATTATGCACTG‑‑CCCTACCCGGTAg                                                                                   >  2:1622693/1‑96 (MQ=255)
                 atTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACg                                                                                >  1:2003326/1‑99 (MQ=255)
                 atTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTGTGTGGACTG‑‑CGCCCCCCGGTAGACGa                                                                               >  1:871425/1‑100 (MQ=255)
                  ttACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGa                                                                               >  1:1989975/1‑99 (MQ=255)
                   tACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATc                                                                             <  2:1241187/100‑1 (MQ=255)
                    aCTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATcc                                                                            >  1:1509094/1‑100 (MQ=255)
                         gTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGccc                                                                       <  2:1405909/100‑1 (MQ=255)
                           aTCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCtat                                                                    <  1:1166637/100‑1 (MQ=255)
                            tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGGTAGACGATACTGCCCta                                                                     >  2:1462035/1‑100 (MQ=255)
                              gAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAg                                                                  <  1:1211759/100‑1 (MQ=255)
                              gAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGACCCCCCCGGTCGACGATCCTGCCCTCTAg                                                                  >  1:658572/1‑100 (MQ=255)
                                gTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGtt                                                                <  1:1877067/100‑1 (MQ=255)
                                  tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTgg                                                              <  2:393190/100‑1 (MQ=255)
                                   gctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGGTAGACGATCCTGCCCTATAGTTg                                                               <  2:2640863/99‑1 (MQ=255)
                                   gctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGa                                                             <  1:2173835/100‑1 (MQ=255)
                                   gctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATACCCCATTTTGTGGACTAC‑‑CCCCCCCGGTAGACGATCCTGCACTATAGTTg                                                               >  2:2318108/1‑98 (MQ=255)
                                    ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCAGTTTGTGGACTGCCCCCCCCGGTAGACCATCCTGCCCTATAGTTGGAg                                                            >  2:794440/1‑100 (MQ=255)
                                      gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATGATACCCAATTTTGTGGACTG‑‑CCCACCCCGGTAGACGATCCTGCCCTATAGTGTGAGTg                                                          >  1:996972/1‑100 (MQ=255)
                                      gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGCCTG‑‑CCCCGCCCTGAAGACGATCCTGCCCTCTAGTTGGGGTg                                                          >  2:1853978/1‑100 (MQ=255)
                                      gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTg                                                          >  2:1284832/1‑100 (MQ=255)
                                       ctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                         >  2:100949/1‑100 (MQ=255)
                                         tGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGTACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccc                                                       <  1:1002347/100‑1 (MQ=255)
                                         tGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccc                                                       <  2:2160879/100‑1 (MQ=255)
                                         tGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccc                                                       >  1:2594421/1‑100 (MQ=255)
                                          ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcccc                                                      <  2:2358888/100‑1 (MQ=255)
                                           gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccccc                                                     <  1:299234/100‑1 (MQ=255)
                                           gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccccc                                                     <  1:2789775/100‑1 (MQ=255)
                                           gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccccc                                                     >  2:306334/1‑99 (MQ=255)
                                            tGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCt                                                    >  2:2245057/1‑100 (MQ=255)
                                             gagaTCAGGAAGAGGTTACCGAAACGGTGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccccc                                                     <  1:2080063/99‑1 (MQ=255)
                                               gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTTTGGACTGCCCCCCCCGTTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTa                                                   <  1:266906/98‑1 (MQ=255)
                                               gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTAtt                                                 >  2:1205259/1‑100 (MQ=255)
                                               gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTAt                                                  >  1:1590412/1‑99 (MQ=255)
                                                  cAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTAAAGTTGAAGTGCCCACTATTCaa                                              >  2:839706/1‑99 (MQ=255)
                                                  cAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAgg                                             <  1:1623948/100‑1 (MQ=255)
                                                   aGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAgg                                             <  1:908002/100‑1 (MQ=255)
                                                   aGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAgg                                             >  2:1057869/1‑100 (MQ=255)
                                                    ggAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGt                                            >  1:887107/1‑100 (MQ=255)
                                                    ggAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGt                                            >  2:650851/1‑100 (MQ=255)
                                                    ggAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTg                                           >  1:2811409/1‑100 (MQ=255)
                                                     gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTg                                           <  1:1867795/100‑1 (MQ=255)
                                                     gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTg                                           <  2:164546/100‑1 (MQ=255)
                                                      aaGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTgg                                          <  1:1014039/100‑1 (MQ=255)
                                                       agagGTTACCGAAACGATGGACATCATCACCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGc                                         >  2:942141/1‑100 (MQ=255)
                                                        gagGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGc                                         <  1:587638/100‑1 (MQ=255)
                                                          ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGCTAGACGATCCTGCCCTATATCTGGAGCGCCCACGATTcgcgtgtcca                                       >  1:1104346/1‑90 (MQ=255)
                                                           gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGGTAGACGATCCTTCCCTATAGTTGGAGCGCCCCCTATTCAGGTGGCTaa                                      >  1:1220349/1‑100 (MQ=255)
                                                            ttACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTCTTCCGGTGGCCAAAt                                    >  2:2113725/1‑100 (MQ=255)
                                                              aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAAttt                                  >  1:1102755/1‑100 (MQ=255)
                                                               ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTa                                 <  2:1273025/100‑1 (MQ=255)
                                                               ccGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                <  1:1716421/100‑1 (MQ=255)
                                                                 gAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                               <  2:2594432/100‑1 (MQ=255)
                                                                     cGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAAt                           >  2:1341920/1‑100 (MQ=255)
                                                                       aTGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATc                          <  1:2549408/100‑1 (MQ=255)
                                                                       aTGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCa                         >  2:2244212/1‑100 (MQ=255)
                                                                        tGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCa                         <  2:593685/100‑1 (MQ=255)
                                                                        tGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCaa                        >  1:1070299/1‑100 (MQ=255)
                                                                         ggACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCACCCTATTCAGGTGGCTAAATTTAGCAAATCAAt                       <  1:857897/100‑1 (MQ=255)
                                                                          gACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAAGTGGCTACATTTAGCAAATCAATa                      >  1:1104151/1‑100 (MQ=255)
                                                                          gACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGTAGACGCACCTGCCCTACAGTTGGAGTGCCCCCTATTCAGCTGGCCAAATCCTGCAAATCagt                       >  1:235850/1‑97 (MQ=255)
                                                                           aCATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                     <  1:1888527/100‑1 (MQ=255)
                                                                            catcatCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGTAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                    >  1:1295409/1‑100 (MQ=255)
                                                                             atcatcCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATacac                   <  2:2518852/100‑1 (MQ=255)
                                                                              tcatcCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCGATTCAGGTGGCTAAATTTAGCAAATCAATacac                   >  1:1528240/1‑99 (MQ=255)
                                                                               catcCCCCATTTTGTGGACTGCCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATacac                   <  1:428328/99‑1 (MQ=255)
                                                                                atcCCCCATTTTGTGGACTGCCCCCCCCCGGTAGACGCTCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACtt                 >  1:1063329/1‑100 (MQ=255)
                                                                                  cccccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg              >  2:1030455/1‑100 (MQ=255)
                                                                                  cccccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg              >  2:2846918/1‑100 (MQ=255)
                                                                                   ccccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTGGGAGTGCCCCCGATTCAGGTGGCTAGATTTAGCAAATCAATACACGTCAg              >  2:858937/1‑99 (MQ=255)
                                                                                            gtgGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtatta    <  1:457020/100‑1 (MQ=255)
                                                                                            gtgGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtatta    >  1:662813/1‑100 (MQ=255)
                                                                                             tgGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtattat   <  2:2003333/100‑1 (MQ=255)
                                                                                             tgGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtattat   >  1:927464/1‑100 (MQ=255)
                                                                                             tgGACTGACCCCACCGGTAGACGATCCTGCCCTACAGTTGGGGTGCACCCTATTCAGGTGGCTAAATTTAGCAAATCAATACTCTTCAGGGGGGtat      <  1:935655/97‑1 (MQ=255)
                                                                                              ggACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTattatt  >  2:889616/1‑100 (MQ=255)
                                                                                                     |                                                                                                 
TTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATT  >  NC_000913/2765688‑2765880

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: