Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,702,837 G→C 15.9% P11R (CCG→CGG)  dppF ← dipeptide/heme ABC transporter ATPas

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,702,8370GC15.9% 610.2 / 100.0 225P11R (CCG→CGG) dppFdipeptide/heme ABC transporter ATPas
Reads supporting (aligned to +/- strand):  ref base G (85/104);  new base C (16/20);  total (101/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

Only 100 of 233 total aligned reads displayed.

TTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATT  >  NC_000913/3702738‑3702936
                                                                                                   |                                                                                                   
ttAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCg                                                                                                     <  2:1718745/100‑1 (MQ=255)
 tAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCgg                                                                                                    >  1:1834626/1‑100 (MQ=255)
  aaacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgt                                                                                                   >  1:3028934/1‑100 (MQ=255)
  aaacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgt                                                                                                   >  2:882089/1‑100 (MQ=255)
   aacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgtt                                                                                                  >  2:1886488/1‑100 (MQ=255)
   aacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgtt                                                                                                  >  1:1826622/1‑100 (MQ=255)
    acgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgttg                                                                                                 >  1:1612202/1‑100 (MQ=255)
    acgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgttg                                                                                                 >  1:2748294/1‑100 (MQ=255)
     cgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgttgt                                                                                                >  1:3117111/1‑100 (MQ=255)
       aaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACGGCGGttgttg                                                                                              >  2:2747387/1‑100 (MQ=255)
          cgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCATGTTTTTTGCCTGCAACAGCGGTTGTTGCa                                                                                            >  2:837550/1‑99 (MQ=255)
          cgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCan                                                                                           <  1:669796/100‑2 (MQ=255)
           gCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAgg                                                                                          <  1:861537/100‑1 (MQ=255)
           gCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAgg                                                                                          <  2:3906055/100‑1 (MQ=255)
            ccATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCTGATAATGTGTTTTCAGGTCGATAGCCTGCAACAACCACTGTTGCAggc                                                                                         <  2:3553193/100‑2 (MQ=255)
            ccATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAggg                                                                                         >  1:3790864/1‑100 (MQ=255)
            ccATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAggg                                                                                         >  1:246185/1‑100 (MQ=255)
            ccATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAggg                                                                                         <  1:4014302/100‑1 (MQ=255)
            ccATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAggc                                                                                         <  1:1803162/100‑2 (MQ=255)
            ccATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAggc                                                                                         <  1:3617165/100‑2 (MQ=255)
            ccATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAgg                                                                                          <  1:2390498/99‑1 (MQ=255)
             cATCCAGCGCTTTAAGCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAggc                                                                                         <  2:2453239/99‑2 (MQ=255)
             cATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAggct                                                                                        <  1:2976314/100‑3 (MQ=255)
               tCCAGCGCTTTAACCAGACGTTCCTGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTGCAGGTCGATAGCCTTCAATAGCGGCTGTTGCAGGGTgg                                                                                      <  2:2584530/100‑1 (MQ=255)
                 cAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGcc                                                                                    >  2:1379938/1‑100 (MQ=255)
                   gcgcTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCctc                                                                                  >  2:2792982/1‑100 (MQ=255)
                    cgcTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCctct                                                                                 <  2:731656/100‑1 (MQ=255)
                    cgcTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAggctatcgacc                                                                                 <  2:1742574/100‑10 (MQ=37)
                      cTTTAACCAGGCGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTGTTg                                                                               >  1:336322/1‑100 (MQ=255)
                      cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAgggagatcggaag                                                                               >  2:3032468/1‑90 (MQ=255)
                      cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                               <  1:1370112/100‑1 (MQ=255)
                       tttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGc                                                                              <  2:1469846/100‑1 (MQ=255)
                       tttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGc                                                                              <  2:3588977/100‑1 (MQ=255)
                        ttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCg                                                                             <  1:2657585/100‑1 (MQ=255)
                          aaCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTa                                                                           >  2:3292047/1‑100 (MQ=255)
                           aCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTAc                                                                          >  1:3524594/1‑100 (MQ=255)
                            ccAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACt                                                                         <  1:2115939/100‑1 (MQ=255)
                            ccAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGACGATAGCCTGCAACAGCGGTTGTTGCAGGGTTGCCTCTTGCGt                                                                            >  2:3177356/1‑97 (MQ=255)
                             cAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTc                                                                        >  1:3823179/1‑100 (MQ=255)
                              aGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCa                                                                       <  2:2633944/100‑1 (MQ=255)
                              aGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCa                                                                       <  1:1043943/100‑1 (MQ=255)
                                  gTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAgt                                                                   <  2:2928889/100‑1 (MQ=255)
                                   ttCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAgtg                                                                  <  2:470020/100‑1 (MQ=255)
                                     ccGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTc                                                                <  1:962368/100‑1 (MQ=255)
                                      cGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCg                                                               <  1:2182983/100‑1 (MQ=255)
                                       ggCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCgg                                                              <  2:92159/100‑1 (MQ=255)
                                        gcgcGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGc                                                             <  1:1080793/100‑1 (MQ=255)
                                        gcgcGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGc                                                             <  1:3954748/100‑1 (MQ=255)
                                         cgcgAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGcc                                                            <  1:2168711/100‑1 (MQ=255)
                                         cgcgAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGcc                                                            >  1:881970/1‑100 (MQ=255)
                                              aCATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCg                                                       <  1:3041514/100‑1 (MQ=255)
                                               cATGCCTTTCTTCACCGGATAATGTTTTTTTTTTTTTTTATCTTGCAAGAGCGGTTGTTGGAGGTTGTCCTCTTGCGTGATCATAGGGTTGGTCTtcc                                                        >  2:4140211/1‑98 (MQ=255)
                                               cATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCgg                                                      <  2:3753854/100‑1 (MQ=255)
                                                 tGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGca                                                    <  1:2792989/100‑1 (MQ=255)
                                                 tGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGca                                                    >  1:4002665/1‑100 (MQ=255)
                                                 tGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGca                                                    >  1:4034633/1‑100 (MQ=255)
                                                  gCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcat                                                   <  2:6261/100‑1 (MQ=255)
                                                   ccTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatc                                                  <  2:161414/100‑1 (MQ=255)
                                                    cTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatca                                                 >  1:621464/1‑100 (MQ=255)
                                                    cTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatca                                                 <  2:4012689/100‑1 (MQ=255)
                                                     tttCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatcat                                                <  1:1966620/100‑1 (MQ=255)
                                                      ttcttcACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCatcatc                                               <  1:3482861/100‑1 (MQ=255)
                                                      ttcttcACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCatcatc                                               <  2:2116867/100‑1 (MQ=255)
                                                       tcttcACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCAtcatca                                              <  1:4138926/100‑1 (MQ=255)
                                                       tcttcACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCAtcatca                                              <  1:2512296/100‑1 (MQ=255)
                                                        cttcACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCaa                                             >  2:334387/1‑100 (MQ=255)
                                                          tcACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGt                                           <  2:2463644/100‑1 (MQ=255)
                                                           cACCGGATAATGTTTTTTTTTTTTTTTTGCCAGCAACAGCGGTCATTGCAGGGTGGCCTATTGCGTAGTCATAGTGTCGGCCTCTCGGCATCATCAAg                                            >  2:1529357/1‑98 (MQ=255)
                                                           cACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTg                                          >  1:3284640/1‑100 (MQ=255)
                                                           cACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTg                                          <  2:638998/100‑1 (MQ=255)
                                                            aCCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTgg                                         >  2:4098879/1‑100 (MQ=255)
                                                                     tGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCa                                >  1:1913663/1‑100 (MQ=255)
                                                                      gttgttTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCAt                               <  2:2966463/100‑1 (MQ=255)
                                                                         ttttCAGGTCGATAGCCTGCAACAGCG‑‑‑GTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGAc                         <  2:204085/100‑1 (MQ=255)
                                                                           ttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGa                          <  1:3926658/100‑1 (MQ=255)
                                                                            taaGGTCGAGAGCCTGCGACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTATTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGAc                         <  2:1063789/98‑1 (MQ=255)
                                                                             cAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACt                        <  1:2880193/100‑1 (MQ=255)
                                                                             cAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACt                        <  1:3671232/100‑1 (MQ=255)
                                                                              aGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTg                       >  2:2155217/1‑100 (MQ=255)
                                                                                 tCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACg                    <  2:3065631/100‑1 (MQ=255)
                                                                                  cGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCATGTGGGTAATGGCATTTGGACTGACGc                   <  1:2552377/100‑1 (MQ=255)
                                                                                   gATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGcc                  >  1:321897/1‑100 (MQ=255)
                                                                                     cacccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCCTTTGGACTGACGccc                 >  2:2537603/4‑99 (MQ=38)
                                                                                      acccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGt               >  1:52692/3‑100 (MQ=255)
                                                                                       gCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTc              <  1:505848/100‑1 (MQ=255)
                                                                                       gCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTc              >  1:69982/1‑100 (MQ=255)
                                                                                       cccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTc              >  1:234061/2‑100 (MQ=255)
                                                                                        ccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCa             >  1:142523/1‑100 (MQ=255)
                                                                                         cTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAg            <  2:283926/100‑1 (MQ=255)
                                                                                         cTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAg            >  1:3406965/1‑100 (MQ=255)
                                                                                         cTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAg            >  2:3844378/1‑100 (MQ=255)
                                                                                         cTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAg            >  1:1185315/1‑100 (MQ=255)
                                                                                          tGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGc           >  2:3469606/1‑100 (MQ=255)
                                                                                          tGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGc           <  1:1119049/100‑1 (MQ=255)
                                                                                           gCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCg          <  2:1169758/100‑1 (MQ=255)
                                                                                              aCAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGc       >  2:1037545/1‑100 (MQ=255)
                                                                                                aGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCat     >  1:2464782/1‑100 (MQ=255)
                                                                                                  cGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGGGCGTCAgtccaaatgc   <  1:1785000/100‑10 (MQ=255)
                                                                                                  cGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCatat   >  1:1688684/1‑100 (MQ=255)
                                                                                                   ggTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGGATTTGGACTGACGCCAGTCAGCGAgga      >  2:3260135/1‑94 (MQ=255)
                                                                                                   |                                                                                                   
TTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATT  >  NC_000913/3702738‑3702936

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