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% 506.3 / 82.2 189P11R (CCG→CGG) dppFdipeptide/heme ABC transporter ATPas
Reads supporting (aligned to +/- strand):  ref base G (70/88);  new base C (17/13);  total (87/102)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-01

Only 100 of 199 total aligned reads displayed.

TTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATT  >  NC_000913/3702738‑3702936
                                                                                                   |                                                                                                   
ttAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCg                                                                                                     >  2:1172002/1‑100 (MQ=255)
 tAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCgg                                                                                                    >  1:2474311/1‑100 (MQ=255)
  aaacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgt                                                                                                   >  1:4003056/1‑100 (MQ=255)
   aacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgtt                                                                                                  >  2:3199148/1‑100 (MQ=255)
     cgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgttgt                                                                                                >  2:2371968/1‑100 (MQ=255)
     cgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgttgt                                                                                                >  1:193079/1‑100 (MQ=255)
     caaaacgCCATCCAGCGCTTTAACCAGACGTTTCTGCTCGCACATGCCTTACTTCACCGCATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgtt                                                                                                  <  1:930875/96‑1 (MQ=255)
      gaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgttgtt                                                                                               >  1:2935167/1‑100 (MQ=255)
       aaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCCAACATGCCTTTCTTCACCGGATAATGTTTTTTCCGGTCGATAGCCTGCAACAGCGGttgttg                                                                                              >  1:1559234/1‑100 (MQ=255)
        aacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGc                                                                                             <  2:184905/100‑1 (MQ=255)
        aacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGc                                                                                             >  1:2522036/1‑100 (MQ=255)
         acgCCATCCAGCGCTTTAACCAGCCGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGttgttg                                                                                              <  1:4187960/98‑1 (MQ=255)
         acgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTACAACAACCGCTGTTGCa                                                                                            <  1:1180714/100‑1 (MQ=255)
          cgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAg                                                                                           <  2:1109961/100‑1 (MQ=255)
           gCCATCCAGCGCTTTAACCAGACGTTCCGGCTCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAgg                                                                                          <  1:2418471/100‑1 (MQ=255)
           gCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATGATGGTTTTGCAGGTCGATAGCCTGCAACAACCGCTGTTGCAgg                                                                                          <  1:4215118/100‑1 (MQ=255)
           gCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAgg                                                                                          <  1:2065405/100‑1 (MQ=255)
           gCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAgg                                                                                          <  2:2998935/100‑1 (MQ=255)
           gCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAgg                                                                                          <  1:554209/100‑1 (MQ=255)
            ccATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAggg                                                                                         <  1:1374204/100‑1 (MQ=255)
             cATCCAGCGCTTTAACCAGATGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGt                                                                                        <  2:3929605/100‑1 (MQ=255)
              aTCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCATTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTg                                                                                       >  2:4291670/1‑100 (MQ=255)
               tCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTgg                                                                                      >  1:277353/1‑100 (MQ=255)
               tCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTgg                                                                                      >  1:3389829/1‑100 (MQ=255)
               tCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAggctat                                                                                      <  1:26864/100‑5 (MQ=38)
                 cAGCGCGTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGcc                                                                                    <  2:1779523/100‑1 (MQ=255)
                  aGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGcc                                                                                    <  1:3676950/99‑1 (MQ=255)
                      cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                               <  2:4003083/100‑1 (MQ=255)
                      cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                               >  2:2879139/1‑100 (MQ=255)
                      cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                               <  2:1871713/100‑1 (MQ=255)
                      cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                               >  2:1946037/1‑100 (MQ=255)
                      cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                               <  1:1951959/100‑1 (MQ=255)
                        ttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                               >  1:778354/1‑98 (MQ=255)
                        ttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCg                                                                             <  2:1385896/100‑1 (MQ=255)
                        ttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCg                                                                             >  1:3321519/1‑100 (MQ=255)
                        ttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCg                                                                             <  2:1191647/100‑1 (MQ=255)
                          aaCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTa                                                                           >  1:3239125/1‑100 (MQ=255)
                          aaCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTa                                                                           >  1:971603/1‑100 (MQ=255)
                           aCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTAc                                                                          <  2:2179372/100‑1 (MQ=255)
                           aCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTAc                                                                          >  2:1608231/1‑100 (MQ=255)
                            ccAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACt                                                                         <  2:1661797/100‑1 (MQ=255)
                              aGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCa                                                                       <  2:2601936/100‑1 (MQ=255)
                                aCGTTCCTTCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATa                                                                     <  1:3948024/100‑1 (MQ=255)
                                 cGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTTTTTTTTTTTACCTGCAACAGAGGTTGTTGCAGGGTGGCCTCTTGCTTACTCATa                                                                     >  1:3069178/1‑99 (MQ=255)
                                 cGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAg                                                                    <  2:2116746/100‑1 (MQ=255)
                                 cGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAg                                                                    <  2:2655428/100‑1 (MQ=255)
                                 cGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTATTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAg                                                                    <  1:3088181/100‑1 (MQ=255)
                                  gTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAgt                                                                   <  2:659850/100‑1 (MQ=255)
                                   ttCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAgtg                                                                  >  2:2019146/1‑100 (MQ=255)
                                      cGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCg                                                               >  2:4239760/1‑100 (MQ=255)
                                       ggCGCGACCATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCgg                                                              <  1:1817436/100‑1 (MQ=255)
                                        gcgcGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGc                                                             >  2:3553457/1‑100 (MQ=255)
                                         cgcgAACATGTCTTTCGTCACCGGATAATTTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGcc                                                            <  1:1748404/100‑1 (MQ=255)
                                         cgcgAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGcc                                                            <  2:4195416/100‑1 (MQ=255)
                                         cgcgAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGcc                                                            >  1:1680182/1‑100 (MQ=255)
                                           cgAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTc                                                          <  1:2382898/100‑1 (MQ=255)
                                             aaCATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTccc                                                        <  1:3376794/100‑1 (MQ=255)
                                               cATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCgg                                                      <  1:1776144/100‑1 (MQ=255)
                                                aTGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGc                                                     <  1:1288248/100‑1 (MQ=255)
                                                 tGCCTTTCTTCACCGGATAATGTTTTTTTTTTTTTTTTTCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGc                                                     >  1:1123608/1‑99 (MQ=255)
                                                 tGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGca                                                    <  1:1532980/100‑1 (MQ=255)
                                                 tGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGca                                                    >  1:2104497/1‑100 (MQ=255)
                                                  gCCTTTCTTCACCGGATAATGTTTTTTTTGTTTGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcat                                                   >  1:3136035/1‑100 (MQ=255)
                                                   ccTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatc                                                  >  2:2343023/1‑100 (MQ=255)
                                                   ccTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatc                                                  <  2:1340023/100‑1 (MQ=255)
                                                   ccTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatc                                                  >  1:1987632/1‑100 (MQ=255)
                                                    cTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatca                                                 <  2:1839274/100‑1 (MQ=255)
                                                      ttcttcACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCatcatc                                               <  1:407745/100‑1 (MQ=255)
                                                        cttcACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCaa                                             >  2:3859647/1‑100 (MQ=255)
                                                             ccGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTgtg                                        <  2:4081463/100‑3 (MQ=255)
                                                              cGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGt                                       >  1:3497156/1‑100 (MQ=255)
                                                                    aTGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTGGTTGCAGGGTGGCCTCTTGAGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGc                                 <  1:2188111/100‑1 (MQ=255)
                                                                     tGTTTTTTTAGGTCGATAGCTTGTAACAGCGTTTGTTGCACAGTCACCTCTTCCGTACTCATAGTGTCAGCCTCCTGGCATCATCAAGTGGGTAATAGc                                 >  1:491966/1‑99 (MQ=255)
                                                                     tGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCa                                <  2:224441/100‑1 (MQ=255)
                                                                      gTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCAt                               <  2:904694/100‑1 (MQ=255)
                                                                         ttttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTg                            <  2:1987640/100‑1 (MQ=255)
                                                                         ttttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTg                            >  2:1196689/1‑100 (MQ=255)
                                                                          tttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTgg                           >  2:323155/1‑100 (MQ=255)
                                                                          tttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTg                            <  2:1288336/99‑1 (MQ=255)
                                                                          tttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTg                            >  2:2056200/1‑99 (MQ=255)
                                                                             cAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCTAGTGGGTAATGGCATTTGGACt                        <  1:2556030/100‑1 (MQ=255)
                                                                             cAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACt                        <  2:2042232/100‑1 (MQ=255)
                                                                                gTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGAc                     <  2:3330701/100‑1 (MQ=255)
                                                                                 tCGAGAGCCTGCAACAGCAGGTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACAGACg                    <  1:3573335/100‑1 (MQ=255)
                                                                                 aggccacccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACg                    >  1:3794028/8‑100 (MQ=38)
                                                                                  ggccacccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGc                   >  2:2515322/7‑100 (MQ=38)
                                                                                   gATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGcc                  <  1:2979627/100‑1 (MQ=255)
                                                                                     cacccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCACGTGGGTAATGGCATTTGGACTGACGccc                 >  1:2088775/4‑99 (MQ=38)
                                                                                     cacccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGAc                     >  1:1652601/4‑95 (MQ=38)
                                                                                     cacccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCg                >  2:889498/4‑100 (MQ=38)
                                                                                       cccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTc              >  1:676700/2‑100 (MQ=255)
                                                                                       cccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTc              >  1:4330869/2‑100 (MQ=255)
                                                                                       cccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTc              >  2:3099024/2‑100 (MQ=255)
                                                                                       cccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTc              >  2:1597732/2‑100 (MQ=255)
                                                                                        ccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCa             >  1:81989/1‑100 (MQ=255)
                                                                                         cTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCACCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAg            >  2:734311/1‑100 (MQ=255)
                                                                                           gCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCg          <  2:3686564/100‑1 (MQ=255)
                                                                                              aCAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGc       <  1:2029022/100‑1 (MQ=255)
                                                                                                  cGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCatat   <  2:3136052/100‑1 (MQ=255)
                                                                                                   ggTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATAtt  >  1:1581994/1‑100 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: