Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,765,794 (C)7→8 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.C97.7% 464.6 / 3.4 133intergenic (‑11/‑131)yfjM/rnlACP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; RNase LS
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (66/62);  total (69/64)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 154 total aligned reads displayed.

TTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTG  >  NC_000913/2765688‑2765881
                                                                                                    |                                                                                                 
ttGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTg                                                                                                    >  2:1141546/1‑100 (MQ=255)
  gTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cc                                                                                                >  2:2517600/1‑100 (MQ=255)
  gTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cc                                                                                                >  1:2758196/1‑100 (MQ=255)
   tCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑c                                                                                                 <  2:2527169/98‑1 (MQ=255)
    ctgctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccc                                                                                             >  1:1869757/1‑100 (MQ=255)
      gctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccccc                                                                                           >  1:1646133/1‑100 (MQ=255)
      gctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccccc                                                                                           <  2:2924400/100‑1 (MQ=255)
       ctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccccc                                                                                           >  2:1373994/1‑99 (MQ=255)
       ctgGCAAATAAGTAATAAGTATCGAGTGGATGCTTGGTGAGATCCGGAAGAGGTTACCGAAACGAGGGACATCAGCCCCCATTTTGTGGACTG‑‑cc                                                                                                <  2:2840108/95‑1 (MQ=255)
            gAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccccggta                                                                                    <  2:853364/100‑6 (MQ=255)
            gAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccccggta                                                                                    <  1:1632213/100‑6 (MQ=255)
            gAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccccggta                                                                                    <  1:1755095/100‑6 (MQ=255)
            gAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccccggta                                                                                    <  2:1723397/100‑6 (MQ=255)
             aataatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAg                                                                                    <  1:1479458/100‑1 (MQ=255)
              ataatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGa                                                                                   <  2:393961/100‑1 (MQ=255)
               taatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGTAAGAGGTTACCGAAACGATTGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGAc                                                                                  <  2:1632275/100‑1 (MQ=255)
               taatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATTCCCCATTTTGTGGACTGCCCCCCCCGGTAGAc                                                                                  <  1:2835013/100‑1 (MQ=255)
               taatTACTAAGTATCGAGTTGCGGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGAc                                                                                  >  2:1733007/1‑100 (MQ=255)
                 atTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGa                                                                                >  1:2795568/1‑100 (MQ=255)
                     cTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCCTTTTGTGGACTGCCCCCCCCGGTAGACGATcc                                                                             >  1:331250/1‑99 (MQ=255)
                      tAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCt                                                                            >  2:397055/1‑99 (MQ=255)
                       aaGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGc                                                                          <  1:2605903/100‑1 (MQ=255)
                       aaGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGCAGACGATACTGc                                                                          >  1:2105043/1‑100 (MQ=255)
                        aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGcc                                                                         <  1:2301746/100‑1 (MQ=255)
                         gTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGccc                                                                        <  2:1996731/100‑1 (MQ=255)
                         gTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGccc                                                                        <  2:433052/100‑1 (MQ=255)
                         gTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGccc                                                                        <  1:2913883/100‑1 (MQ=255)
                         gTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGccc                                                                        >  1:941500/1‑100 (MQ=255)
                         gTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGccc                                                                        <  1:1894266/100‑1 (MQ=255)
                         gTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCACTTTGTGGACTGACCCCCCCGGTAGACgcagctgcc                                                                         >  1:2370991/1‑91 (MQ=255)
                          tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCt                                                                       <  2:564314/100‑1 (MQ=255)
                          tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGCGCGACGCTCCTGCCCt                                                                       >  1:2921101/1‑100 (MQ=255)
                           aTCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCt                                                                       >  2:2893269/1‑99 (MQ=255)
                           aTCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACCGCCCCCCCCGGTAGCCGATCCTGCGCta                                                                      >  2:2823465/1‑100 (MQ=255)
                            tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCtat                                                                     >  1:2611700/1‑100 (MQ=255)
                             cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCACTTTGTGGACTGCCCCCCCCGGTAGACGCCCCTGCCCt                                                                       >  1:1024294/1‑97 (MQ=255)
                              gAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAg                                                                   <  2:2288958/100‑1 (MQ=255)
                              gAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAg                                                                   >  1:1920810/1‑100 (MQ=255)
                                 ttGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTg                                                                <  1:267173/100‑1 (MQ=255)
                                 ttGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTg                                                                <  2:1118049/100‑1 (MQ=255)
                                  tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCGCATTTTGTGGACTG‑‑CCCCCCACGGTAAACGATCATGCCCTATAGTTgg                                                               >  1:2802047/1‑100 (MQ=255)
                                  tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATCGTTgg                                                               >  1:740170/1‑100 (MQ=255)
                                  tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATATTTgg                                                               >  1:215666/1‑100 (MQ=255)
                                  tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTgg                                                               >  2:1542280/1‑100 (MQ=255)
                                  tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTgg                                                               <  1:1651268/100‑1 (MQ=255)
                                  tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGGATACGATCCTGccc                                                                        >  1:2802683/1‑91 (MQ=255)
                                   gctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGa                                                              <  2:317283/100‑1 (MQ=255)
                                      gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTg                                                           >  1:2260294/1‑100 (MQ=255)
                                      gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTg                                                           >  1:1253530/1‑100 (MQ=255)
                                      gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTg                                                           >  1:142040/1‑100 (MQ=255)
                                      gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTg                                                           <  1:2535992/100‑1 (MQ=255)
                                      gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTg                                                           <  2:1257807/100‑1 (MQ=255)
                                      gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTg                                                           >  1:2881321/1‑100 (MQ=255)
                                         tGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTg                                                           <  2:1810520/98‑1 (MQ=255)
                                         tGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccc                                                        >  1:2581944/1‑100 (MQ=255)
                                          ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcccc                                                       >  2:1180585/1‑100 (MQ=255)
                                            tGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCt                                                     <  2:617839/100‑1 (MQ=255)
                                               gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTAtt                                                  <  1:2351307/100‑1 (MQ=255)
                                                aTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTc                                                 <  2:893483/100‑1 (MQ=255)
                                                aTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCa                                                <  1:2986798/100‑1 (MQ=255)
                                                     gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTg                                            <  1:1006432/100‑1 (MQ=255)
                                                     gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTg                                            <  1:328460/100‑1 (MQ=255)
                                                       agagGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGc                                          <  1:2546722/100‑1 (MQ=255)
                                                          ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTCGAGGGCCCCCTATTCAGGTGGCTa                                        >  1:284103/1‑99 (MQ=255)
                                                             tACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAAtt                                    <  2:919281/100‑1 (MQ=255)
                                                             tACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAAtt                                    >  1:2416885/1‑100 (MQ=255)
                                                              aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAAttt                                   <  2:1024289/100‑1 (MQ=255)
                                                               ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTa                                  >  2:2976037/1‑100 (MQ=255)
                                                               ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTa                                  >  1:1518229/1‑100 (MQ=255)
                                                               ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTa                                  >  2:1096910/1‑100 (MQ=255)
                                                               ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTa                                  >  2:1715221/1‑100 (MQ=255)
                                                                 gAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                <  1:2691227/100‑1 (MQ=255)
                                                                 gAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                <  2:1969190/100‑1 (MQ=255)
                                                                 gAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                <  2:331248/100‑1 (MQ=255)
                                                                     cGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAAt                            >  2:1330866/1‑100 (MQ=255)
                                                                      gATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATc                           >  2:90221/1‑100 (MQ=255)
                                                                        tGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCaa                         <  1:473544/100‑1 (MQ=255)
                                                                          gACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                       <  1:1267452/100‑1 (MQ=255)
                                                                          gACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                       <  2:2611685/100‑1 (MQ=255)
                                                                          gACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                       >  1:655495/1‑100 (MQ=255)
                                                                           aCATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                      <  1:1944235/100‑1 (MQ=255)
                                                                           aCATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                      <  1:1542289/100‑1 (MQ=255)
                                                                           aCATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                      >  2:1871236/1‑100 (MQ=255)
                                                                            catcatCCCCCATATTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                     >  1:2026428/1‑100 (MQ=255)
                                                                               catcGCCCGTGTTGTGGACTG‑‑CCGCCCCCGGTAGACGATCCTGCCCTAGAGTTGAAGTGCCCCCGATTCAGGTGGCTAAATTTAGCAAATCAAGacac                    <  2:2861451/98‑1 (MQ=255)
                                                                                atcCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTc                 >  1:2287969/1‑100 (MQ=255)
                                                                                   ccccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgg              >  2:763576/1‑100 (MQ=255)
                                                                                    ctcttttgGTGGATTGCCCCCCCCGGTAGACGAGCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg               <  2:2370977/94‑1 (MQ=255)
                                                                                    cccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCCATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg             >  2:427913/1‑100 (MQ=255)
                                                                                     ccATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg             <  1:604655/99‑1 (MQ=255)
                                                                                       aTTTTGTGGACAGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTAGCTAAATTTAGCAAATCAATACACTTCAgggggg          >  2:1575888/1‑100 (MQ=255)
                                                                                        ttttGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg          <  2:1897050/99‑1 (MQ=255)
                                                                                         tttGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta        <  2:2986615/100‑1 (MQ=255)
                                                                                          ttGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtat       >  1:2959321/1‑100 (MQ=255)
                                                                                          ttGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtat       >  1:432523/1‑100 (MQ=255)
                                                                                          ttGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtat       >  1:2235322/1‑100 (MQ=255)
                                                                                            gtgGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtatta     >  2:2683291/1‑100 (MQ=255)
                                                                                            gtgGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtatta     <  2:1253524/100‑1 (MQ=255)
                                                                                              ggACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTattatt   >  1:1708015/1‑100 (MQ=255)
                                                                                               gACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTACCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTg  >  2:2326328/1‑100 (MQ=255)
                                                                                                    |                                                                                                 
TTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTG  >  NC_000913/2765688‑2765881

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