Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,414,756 (T)6→5 100% coding (38/810 nt) recT ← Rac prophage; recombination and repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,414,7510T.58.1% 44.4 / 249.6 124coding (43/810 nt)recTRac prophage; recombination and repair protein
Reads supporting (aligned to +/- strand):  ref base T (25/27);  new base . (29/43);  total (54/70)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 133 total aligned reads displayed.

GGAAGAGCTGCTGCCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAGCTGGCCATTCA  >  NC_000913/1414654‑1414849
                                                                                                 |                                                                                                  
ggAAGAGCTGCTGCCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGttt                                                                                                  >  2:899735/1‑100 (MQ=255)
    gagCTGCTGCCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGtttttgc                                                                                              <  1:1211551/100‑3 (MQ=255)
    gagCTGCTGCCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGtttttgc                                                                                              <  1:3939266/100‑3 (MQ=255)
    gagCTGCTGCCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTg                                                                                              <  1:3725401/100‑1 (MQ=255)
     agCTGCTGCCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGtttttgc                                                                                              >  2:1058925/1‑97 (MQ=255)
     agCTGCTGCCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGc                                                                                             <  2:2260576/100‑1 (MQ=255)
      gctgctGCCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGtttttgcag                                                                                            >  2:1186491/1‑96 (MQ=255)
       ctgctgCCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAg                                                                                           >  2:2957518/1‑100 (MQ=255)
         gctgcCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATc                                                                                        <  1:1467641/100‑1 (MQ=255)
            gcCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGAt                                                                                         >  2:2662220/1‑96 (MQ=255)
            gcCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGc                                                                                     >  1:55067/1‑100 (MQ=255)
              cAGTTGCTCTTTCATTGATGGCTGGTTACTAAAACTAATCACTTCGCTATTTTTAACTGCTGCTAGTGAACGGTTTCCCTGAGTTTTTGCAGATCGGCt                                                                                    >  2:755510/1‑99 (MQ=255)
                gTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCttt                                                                                  <  2:3519259/100‑1 (MQ=255)
                 ttGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCtttt                                                                                 <  1:2768697/100‑1 (MQ=255)
                  tGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGc                                                                               <  1:406842/100‑1 (MQ=255)
                  tGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGc                                                                               <  1:2446360/100‑1 (MQ=255)
                  tGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGc                                                                               >  2:2176361/1‑100 (MQ=255)
                   gCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGTACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGc                                                                               <  1:2336532/100‑1 (MQ=255)
                    ctctTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCg                                                                              >  1:3588927/1‑100 (MQ=255)
                    ctctTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGa                                                                             >  1:2848450/1‑100 (MQ=255)
                     tctTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGAt                                                                            <  1:3739784/100‑1 (MQ=255)
                     tctTTCATGGATGGCTGGTTAATAAAACTAATTGGGTCGCGATTTTTAACTGCTGTTGGTTCATGGTTTCCCTGAGTTTTTGCAGATTGGCTTTTGCGa                                                                             <  2:2172059/99‑1 (MQ=255)
                       tttCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTATCTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGAtt                                                                           <  2:4111277/99‑1 (MQ=255)
                       tttCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATtg                                                                          >  2:2915932/1‑100 (MQ=255)
                        ttCATTGATGGCTGGTTAATAAAACTAATCGCGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATtgg                                                                         <  2:902064/100‑1 (MQ=255)
                        ttCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATtg                                                                          >  2:3715524/1‑100 (MQ=255)
                          cATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATtggt                                                                        <  1:1746869/100‑1 (MQ=255)
                          cATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATtggtg                                                                       >  2:2184849/1‑100 (MQ=255)
                          cATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATtggtg                                                                       <  2:279015/100‑1 (MQ=255)
                           aTTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATtggtg                                                                       <  1:1085920/100‑1 (MQ=255)
                           aTTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATtggtg                                                                       <  1:3404438/100‑1 (MQ=255)
                           aTTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATtggtgg                                                                      <  2:551371/100‑1 (MQ=255)
                           aTTGATGGCTGGTTAATAAAACTAATCACGTCGCTAGTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATtggtgg                                                                      <  2:1639475/100‑1 (MQ=255)
                            ttGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATtggtgg                                                                      <  1:2637367/100‑1 (MQ=255)
                            ttGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATtggtgg                                                                      >  2:3864472/1‑100 (MQ=255)
                            ttGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTggtggt                                                                     <  1:1367559/100‑1 (MQ=255)
                             tGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGtt                                                                    <  1:2557537/100‑1 (MQ=255)
                             tGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGtt                                                                    <  2:1543800/100‑1 (MQ=255)
                              gATGTCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGtt                                                                    <  2:3300615/100‑1 (MQ=255)
                              gATGGCTGGTTATTAATACTACTCTCGTCGCTGTTTTCACTGCTGTTGTTGCACGTTTTCCCCGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGc                                                                  <  2:1357801/100‑1 (MQ=255)
                              gATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGtt                                                                    >  1:8835/1‑100 (MQ=255)
                                tGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGc                                                                  >  2:657145/1‑100 (MQ=255)
                                tGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGc                                                                  <  1:3869172/100‑1 (MQ=255)
                                 ggCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCt                                                                 >  1:2399064/1‑100 (MQ=255)
                                 ggCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCtt                                                                <  1:3008664/100‑1 (MQ=255)
                                  gCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCtt                                                                <  1:1743608/100‑1 (MQ=255)
                                   cTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTa                                                               >  1:2001618/1‑100 (MQ=255)
                                       ttAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCAt                                                          >  1:2972807/1‑100 (MQ=255)
                                        taataaAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCAtt                                                         >  2:597924/1‑100 (MQ=255)
                                         aataaAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCAtt                                                         >  1:950498/1‑100 (MQ=255)
                                          ataaAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTg                                                       >  2:4054673/1‑100 (MQ=255)
                                             aaaCTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCat                                                    >  2:2698717/1‑100 (MQ=255)
                                                   aTCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGCTGCTTAGTCATTTGCATATTcc                                               >  1:1781864/1‑99 (MQ=255)
                                                    tCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCtt                                             <  1:178675/100‑1 (MQ=255)
                                                     cACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTa                                            >  1:2366132/1‑100 (MQ=255)
                                                      aCGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAg                                           <  1:2839414/100‑1 (MQ=255)
                                                         tCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGcc                                         <  1:395620/100‑1 (MQ=255)
                                                           gCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAg                                      >  2:669092/1‑100 (MQ=255)
                                                            cTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGc                                     >  1:3414521/1‑100 (MQ=255)
                                                             tATTTTTAACTGCTGCTGGTGCACGGTTTGCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGCTTGCTTAGTCATTTGCATATTCCTTAGCCCAGc                                     <  2:2484735/100‑1 (MQ=255)
                                                              aTTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGc                                     >  2:2499176/1‑99 (MQ=255)
                                                                ttttAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCggg                                  >  1:655046/1‑100 (MQ=255)
                                                                ttttAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCgggg                                 >  1:3335894/1‑100 (MQ=255)
                                                                 tttAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGc                                >  1:1738421/1‑100 (MQ=255)
                                                                   tAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAg                              <  1:2625374/100‑1 (MQ=255)
                                                                    aaCTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGCTTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCTAGCgg                                   >  2:946748/1‑95 (MQ=255)
                                                                    aaCTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGt                             >  1:551227/1‑100 (MQ=255)
                                                                     aCTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGt                             <  2:2945582/100‑1 (MQ=255)
                                                                      ctgctgctgGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTgt                           >  2:3486033/1‑99 (MQ=255)
                                                                       tgctgctgGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGa                           <  1:4124913/100‑1 (MQ=255)
                                                                        gctgctgGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATa                         >  2:2266513/1‑100 (MQ=255)
                                                                        gctgctgGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTCAGCCCAGCGGGGCAGTGATa                         <  2:3812861/100‑1 (MQ=255)
                                                                        gctgctgGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTCAGCCCAGCGGGGCAGTGATa                         <  2:2520842/100‑1 (MQ=255)
                                                                         ctgctgGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATa                         >  1:1922795/1‑100 (MQ=255)
                                                                         ctgctgGTGCACGGTTTCCCTGAGTTTTTGCAGATCTGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATaa                        <  1:2902426/100‑1 (MQ=255)
                                                                         ctgctgGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATaa                        <  1:1269302/100‑1 (MQ=255)
                                                                          tgctgGTGCACGGTTTCCTTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTACTCCTCTGCATGTTCATTTGCCCAGCGGGGCAgcgcat                         >  2:3928063/1‑93 (MQ=255)
                                                                          tgctgGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAAt                       >  1:577520/1‑100 (MQ=255)
                                                                          tgctgGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAAt                       <  1:2134725/100‑1 (MQ=255)
                                                                           gctgGTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGTTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTTATAATg                      <  2:2404314/100‑1 (MQ=255)
                                                                              ggTGCACGGTTTCCCTGAGTTTTTTTTTTTTTTTTTTTGCGAATGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCTGGGCTGTGATAATGTCt                   >  2:3477746/1‑100 (MQ=255)
                                                                              ggTGCACGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCTGTGATAATGTCt                   <  1:2715012/100‑1 (MQ=255)
                                                                               gTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCt                   >  1:1591847/1‑100 (MQ=255)
                                                                                tGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCtt                  <  1:803103/100‑1 (MQ=255)
                                                                                tGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCtt                  >  1:170836/1‑100 (MQ=255)
                                                                                 gCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTa                 >  1:102184/1‑100 (MQ=255)
                                                                                   aCGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAAt               >  2:1276981/1‑100 (MQ=255)
                                                                                   aCGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATa              <  1:2662236/100‑1 (MQ=255)
                                                                                   aCGGTTTCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATa              >  1:2069018/1‑100 (MQ=255)
                                                                                    cGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATa              >  1:1475452/1‑100 (MQ=255)
                                                                                    cGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATa              <  2:345241/100‑1 (MQ=255)
                                                                                    cGGTTTCCCTGATTTTTTGTAGATCGGCTTTTGCGATTGGTTGTTGGTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATa              <  2:3778168/99‑1 (MQ=255)
                                                                                      gTTGCCCTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAGCt           <  2:3195671/100‑1 (MQ=255)
                                                                                          cccTGAGTTTTTTTCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAGCTGGc        <  2:3342219/100‑1 (MQ=255)
                                                                                          cccTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAGCTGGc        <  1:164779/100‑1 (MQ=255)
                                                                                           ccTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAGCTGGcc       >  2:3720596/1‑100 (MQ=255)
                                                                                           ccTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAGCTGGCCa      >  1:2362436/1‑100 (MQ=255)
                                                                                            cTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAGCTGGCCa      <  1:2499190/100‑1 (MQ=255)
                                                                                            cTGAGTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATCATGTCTTAATAGCTGGCCAt     >  2:813516/1‑100 (MQ=255)
                                                                                                gTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAGCTGGCCATTCa  >  1:932113/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
GGAAGAGCTGCTGCCAGTTGCTCTTTCATTGATGGCTGGTTAATAAAACTAATCACGTCGCTATTTTTAACTGCTGCTGGTGCACGGTTTCCCTGAGTTTTTTGCAGATCGGCTTTTGCGATTGGTGGTTGCTTAGTCATTTGCATATTCCTTAGCCCAGCGGGGCAGTGATAATGTCTTAATAGCTGGCCATTCA  >  NC_000913/1414654‑1414849

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: