Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 416,389 (C)5→4 100% coding (564/1203 nt) sbcD ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913416,3850C.96.0% 437.6 / 8.4 100coding (568/1203 nt)sbcDexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (42/54);  total (44/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

Only 100 of 106 total aligned reads displayed.

TATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGG  >  NC_000913/416286‑416480
                                                                                                   |                                                                                               
tatGCCCGAGCGCGATGTAGTCGGTTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGc                                                                                                 <  1:422845/100‑1 (MQ=255)
  tGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGccc                                                                                               >  2:2389258/1‑100 (MQ=255)
  tGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGccc                                                                                               >  1:1619158/1‑100 (MQ=255)
    ctcGAGCTCTATTTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGcccc                                                                                              <  1:1474648/97‑1 (MQ=255)
    cccGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGcccca                                                                                             <  2:1927344/100‑2 (MQ=255)
     ccGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCa                                                                                            >  2:2189445/1‑100 (MQ=255)
      cGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGccccacg                                                                                           <  2:2006511/100‑4 (MQ=255)
         gcgcGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTc                                                                                       <  1:472525/100‑1 (MQ=255)
          cgcgATTTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACTGCTTCACTTTTACTGGCCCCACGGTCg                                                                                      <  2:1810624/100‑1 (MQ=255)
          cgcgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCg                                                                                      >  2:2703862/1‑100 (MQ=255)
           gcgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGt                                                                                     >  2:1720681/1‑100 (MQ=255)
            cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGtt                                                                                    >  1:180269/1‑100 (MQ=255)
            cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGtt                                                                                    >  2:2423361/1‑100 (MQ=255)
            cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGtt                                                                                    <  2:2602043/100‑1 (MQ=255)
             gATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTa                                                                                   <  2:449692/100‑1 (MQ=255)
             gATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTa                                                                                   <  1:2884128/100‑1 (MQ=255)
                gTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAAt                                                                                <  1:2946213/100‑1 (MQ=255)
                 tAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCATCGTGCCAATATAAATGTCACGCACTGCGTCATTTTTACTGGCCTCACGGTCGTTAAATg                                                                               <  2:2892465/100‑1 (MQ=255)
                 tAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCACTATAAATGTCACGCACGGCGTCATTTTTACTGGCTCCACGGTCGTTAAATg                                                                               <  2:2750734/100‑1 (MQ=255)
                   ggcGGCTGGTGGCAAGTTTTGTGCCGGATACGGGGCCAGCGGGCCAATATAGATGGCACTCACTGCGTTACTTTTACTGGCCCCACGGTCGTTAAATg                                                                               <  2:1041107/96‑1 (MQ=255)
                   gTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTc                                                                             >  2:2935344/1‑100 (MQ=255)
                   gTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGAACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTc                                                                             >  1:2581823/1‑100 (MQ=255)
                    tCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTcc                                                                            >  2:1440400/1‑100 (MQ=255)
                      ggCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCg                                                                          >  1:1144416/1‑100 (MQ=255)
                       gCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGt                                                                         >  1:1776094/1‑100 (MQ=255)
                          ggtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTACCGTGGc                                                                      >  2:1750729/1‑100 (MQ=255)
                           gtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCg                                                                     >  1:568121/1‑100 (MQ=255)
                           gtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCg                                                                     >  1:568120/1‑100 (MQ=255)
                           gtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCg                                                                     >  2:2437061/1‑100 (MQ=255)
                           gtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCg                                                                     >  1:206366/1‑100 (MQ=255)
                           gtggAAAGTTTTGTGCCGGAAACGCGTACAGCTTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGGCCCGTGGCg                                                                     <  2:2995176/100‑1 (MQ=255)
                            tggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCga                                                                    >  1:494729/1‑100 (MQ=255)
                             ggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCgat                                                                   >  2:21664/1‑100 (MQ=255)
                             ggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGACCCACGGTCGTTAAATGTCCCGTGGc                                                                      >  1:2599244/1‑97 (MQ=255)
                              gAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCgatg                                                                  <  1:62379/100‑1 (MQ=255)
                                aaGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCgatgat                                                                <  2:2068871/100‑1 (MQ=255)
                                  gTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGTCCCCACGGTCGTTAAATGTCCCTTGGCGATGATgg                                                              >  2:1218568/1‑100 (MQ=255)
                                   ttttGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATggg                                                             >  2:2571437/1‑100 (MQ=255)
                                   ttttGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATggg                                                             >  1:1576493/1‑100 (MQ=255)
                                     ttGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCa                                                           <  1:2437090/100‑1 (MQ=255)
                                     ttGTCCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACTGCGTCACTTTTACTGGCCCCACGGTCTTTAAATGTCCCGTGGCGATGATGGGCa                                                           <  2:270616/100‑1 (MQ=255)
                                      tgtgCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCag                                                          <  1:2658283/100‑1 (MQ=255)
                                        tgCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCagag                                                        >  1:2473695/1‑100 (MQ=255)
                                         gCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAgg                                                       >  2:551581/1‑100 (MQ=255)
                                         gCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACAGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAgg                                                       >  2:1652768/1‑100 (MQ=255)
                                          ccGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGc                                                      >  1:695022/1‑100 (MQ=255)
                                          ccGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGc                                                      >  2:572193/1‑100 (MQ=255)
                                          ccGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGc                                                      <  1:326392/100‑1 (MQ=255)
                                          ccGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGc                                                      <  2:589924/100‑1 (MQ=255)
                                          ccGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGc                                                      >  2:356679/1‑100 (MQ=255)
                                            ggAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTg                                                    >  1:138447/1‑100 (MQ=255)
                                            ggAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTg                                                    <  2:1576475/100‑1 (MQ=255)
                                             ggaaCGCGTCCAGAGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCACGTGGCGATGATAGGCAGAGGCTGa                                                   >  1:2060352/3‑100 (MQ=255)
                                              aaaCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGAt                                                  <  1:572196/100‑1 (MQ=255)
                                               aaCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGTGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATc                                                 <  1:2551713/100‑1 (MQ=255)
                                               aaCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTTACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGTCAGAGGCTGATc                                                 <  1:1996439/100‑1 (MQ=255)
                                               aaCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATc                                                 <  2:2256673/100‑1 (MQ=255)
                                               aaCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATc                                                 <  1:493020/100‑1 (MQ=255)
                                               aaCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATc                                                 <  2:2073159/100‑1 (MQ=255)
                                               aaCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATc                                                 >  1:1050804/1‑100 (MQ=255)
                                                  gcgTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATcgcc                                              <  2:84556/100‑1 (MQ=255)
                                                   cgTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATcgccg                                             <  1:893830/100‑1 (MQ=255)
                                                    gTCCAGCGTGCCAATATAAATGTCACGCACGGCGTTATTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATcgccgc                                            <  1:79834/100‑1 (MQ=255)
                                                     tCCAGCTTGCTAATATAAATGTCACGCACGGTGTCATTTTTACTGGTCCTACGGTCGTTAAGTGTCCGGTGGCGATGATGGGCAGAGGCTGATCGCcgcg                                           <  1:563029/100‑1 (MQ=255)
                                                     tCCAGCGTGCCAATATAAATGTCACGCACGGTGTTACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCcgcg                                           <  1:1227056/100‑1 (MQ=255)
                                                     tCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCcgcg                                           >  1:2925461/1‑100 (MQ=255)
                                                        aGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCgcgcn                                         <  2:111473/100‑2 (MQ=255)
                                                          cGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCTCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGtt                                      <  1:2636298/100‑1 (MQ=255)
                                                             gCCAATATAAATGTCACGCACGGGGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGc                                   <  2:2060329/100‑1 (MQ=255)
                                                             gCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGc                                   <  2:365943/100‑1 (MQ=255)
                                                              ccAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCa                                  <  2:1564569/100‑1 (MQ=255)
                                                              ccAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGATGATAGCCGCGCAGTTTGCa                                  >  1:2730490/1‑100 (MQ=255)
                                                               cAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAg                                 <  2:1605585/100‑1 (MQ=255)
                                                                aaTATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAgg                                <  1:2110487/100‑1 (MQ=255)
                                                                aaTATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATATCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAgg                                >  1:2556487/1‑100 (MQ=255)
                                                                  tataAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCa                              >  1:267051/1‑100 (MQ=255)
                                                                  tataAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCa                              >  2:1281353/1‑100 (MQ=255)
                                                                        tGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCa                        <  2:2497066/100‑1 (MQ=255)
                                                                          tCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATa                      <  2:180267/100‑1 (MQ=255)
                                                                             cgcacgGCGTCACTTTTACTGGCCTCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAgtg                   <  1:1461886/100‑1 (MQ=255)
                                                                                  ggCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGTGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGtgttgt               <  1:2321269/99‑1 (MQ=255)
                                                                                  ggCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTgttgtt              <  2:2557091/100‑1 (MQ=255)
                                                                                  ggCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTgttgtt              <  2:2557139/100‑1 (MQ=255)
                                                                                    cGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGTGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTgg            <  2:1355066/100‑1 (MQ=255)
                                                                                    cGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCTATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTgg            <  1:393755/100‑1 (MQ=255)
                                                                                     gTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGt           <  2:433077/100‑1 (MQ=255)
                                                                                      tCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGta          <  1:2423390/100‑1 (MQ=255)
                                                                                       cACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGta          <  2:1024459/100‑1 (MQ=255)
                                                                                       cACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaa         >  1:603069/1‑100 (MQ=255)
                                                                                        aCTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaat        <  2:2622825/100‑1 (MQ=255)
                                                                                         cTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaata       >  2:462611/1‑100 (MQ=255)
                                                                                          ttttACTGGCCCCACGGTCGTTAATTGTCCCGTGGCGATGTTGGGCAGAGGCTGTTCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaataa      <  2:791245/100‑1 (MQ=255)
                                                                                          ttttACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaataa      >  1:2850498/1‑100 (MQ=255)
                                                                                          ttttACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaataa      >  1:1469113/1‑100 (MQ=255)
                                                                                            ttACTGGCCCCACGGTCGTTAAATGTCCCGTGGCTATGATGCGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGTATAGTGTTGTTGGTAATAATc    <  2:2353053/100‑1 (MQ=255)
                                                                                            ttACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATc    <  1:382011/100‑1 (MQ=255)
                                                                                            ttACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATc    <  2:1475385/100‑1 (MQ=255)
                                                                                             tACTGGCCCCACGGTCGTTAAATGTCCCGTGGTGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCg   <  1:1517630/100‑1 (MQ=255)
                                                                                             tACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCg   >  1:2894464/1‑100 (MQ=255)
                                                                                              aCTGGCCCCACGGTCGTTAAATGTCCCGTGGTGATGATGGGTAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCgg  <  1:214229/100‑1 (MQ=255)
                                                                                                   |                                                                                               
TATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGG  >  NC_000913/416286‑416480

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: