Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,702,839 T→C 16.6% Q10Q (CAA→CAG 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,8390TC16.6% 406.7 / 71.2 182Q10Q (CAA→CAGdppFdipeptide/heme ABC transporter ATPas
Reads supporting (aligned to +/- strand):  ref base T (64/86);  new base C (17/13);  total (82/100)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-01

Only 100 of 197 total aligned reads displayed.

AAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCA  >  NC_000913/3702740‑3702938
                                                                                                   |                                                                                                   
aaacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgt                                                                                                     >  1:4003056/1‑100 (MQ=255)
aaacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgt                                                                                                     <  2:2256876/100‑1 (MQ=255)
 aacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgtt                                                                                                    >  2:112768/1‑100 (MQ=255)
     aaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCCAACATGCCTTTCTTCACCGGATAATGTTTTTTCCGGTCGATAGCCTGCAACAGCGGttgttg                                                                                                >  1:1559234/1‑100 (MQ=255)
      aacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGc                                                                                               <  2:184905/100‑1 (MQ=255)
       acgCCATCCAGCGCTTTAACCAGCCGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGttgttg                                                                                                <  1:4187960/98‑1 (MQ=255)
       acgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTACAACAACCGCTGTTGCa                                                                                              <  1:1180714/100‑1 (MQ=255)
        cgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCTAACATGCCTTTCTCCACCGGATAATGTTTTCTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAg                                                                                             <  1:2519487/100‑1 (MQ=255)
        cgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAg                                                                                             >  2:4150688/1‑100 (MQ=255)
        cgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAg                                                                                             <  2:2531670/100‑1 (MQ=255)
         gCCATCCAGCGCTTTAACCAGACGTTCCGGCTCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAgg                                                                                            <  1:2418471/100‑1 (MQ=255)
         gCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATCATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAgg                                                                                            <  1:2602363/100‑1 (MQ=255)
          ccATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAggg                                                                                           >  1:457359/1‑100 (MQ=255)
          ccATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAggg                                                                                           <  1:1374204/100‑1 (MQ=255)
          ccATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAgg                                                                                            <  2:2346608/99‑1 (MQ=255)
          ccATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAggc                                                                                           <  2:2146136/100‑2 (MQ=255)
          ccATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAggc                                                                                           <  2:892387/100‑2 (MQ=255)
             tCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTATTGCAGGGTgg                                                                                        >  1:3389700/1‑100 (MQ=255)
              ccAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGc                                                                                       <  2:717736/100‑1 (MQ=255)
               cAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGcc                                                                                      >  2:3038430/1‑100 (MQ=255)
               cAGCGCGTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGcc                                                                                      <  2:1779523/100‑1 (MQ=255)
                  cgcTTTAACCAGACGTTCCGGCGCGGACATGCCTTTCTTCACCGGATAATGTTTTTGCAGGTCGATAGCCTGCAAGAGCGGTTGTTGCAGGGTGGCctct                                                                                   <  1:1950806/100‑1 (MQ=255)
                  cgcTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCctct                                                                                   >  2:2634277/1‑100 (MQ=255)
                   gcTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCtt                                                                                  <  2:282557/100‑1 (MQ=255)
                   gcTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCtt                                                                                  >  1:89698/1‑100 (MQ=255)
                    cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                                 <  1:1725245/100‑1 (MQ=255)
                    cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                                 >  2:1946037/1‑100 (MQ=255)
                    cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                                 <  2:1871713/100‑1 (MQ=255)
                    cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                                 <  1:2075232/100‑1 (MQ=255)
                      ttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCg                                                                               <  2:1191647/100‑1 (MQ=255)
                      ttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCg                                                                               <  2:1385896/100‑1 (MQ=255)
                      ttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCg                                                                               >  1:3321519/1‑100 (MQ=255)
                        aaCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTa                                                                             >  1:3239125/1‑100 (MQ=255)
                         aCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTAc                                                                            >  2:1608231/1‑100 (MQ=255)
                         aCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTAc                                                                            >  2:2188119/1‑100 (MQ=255)
                          ccAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACt                                                                           <  2:1661797/100‑1 (MQ=255)
                           cAGACGTTCCGGCGCGAACATGCCTTTCTTCCCCGGATAATGTTTTTTCAGTTCTATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTc                                                                          >  2:3748480/1‑100 (MQ=255)
                            aGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCa                                                                         <  1:803453/100‑1 (MQ=255)
                              aCGTTCCTTCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATa                                                                       <  1:3948024/100‑1 (MQ=255)
                               cGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTTTTTTTTTTTACCTGCAACAGAGGTTGTTGCAGGGTGGCCTCTTGCTTACTCATa                                                                       >  1:3069178/1‑99 (MQ=255)
                               cGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAg                                                                      <  2:3615281/100‑1 (MQ=255)
                               cGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTATTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAg                                                                      <  1:3088181/100‑1 (MQ=255)
                                gTTCCGGCGCGAACATGCCTTGCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAgt                                                                     <  1:2782144/100‑1 (MQ=255)
                                  tCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAgtgt                                                                   >  1:4286375/1‑100 (MQ=255)
                                  tCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAgtgt                                                                   >  2:3957486/1‑100 (MQ=255)
                                     ggCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCgg                                                                >  2:2100008/1‑100 (MQ=255)
                                     ggCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCgg                                                                >  1:1924675/1‑100 (MQ=255)
                                       cgcgAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGcc                                                              <  2:4235024/100‑1 (MQ=255)
                                       cgcgAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGcc                                                              <  2:3871640/100‑1 (MQ=255)
                                        gcgAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCt                                                             <  1:523524/100‑1 (MQ=255)
                                         cgAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTc                                                            <  1:2382898/100‑1 (MQ=255)
                                           aaCATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTccc                                                          <  1:3376794/100‑1 (MQ=255)
                                              aTGCCTTTCTTCACCGGATAATGTTTTTTTTTTTTTTTTTCTTCCAACAGCGGTTGTTTCATTGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGc                                                       >  1:1988503/1‑100 (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:2104618/1‑100 (MQ=255)
                                                gCCTTTCTTCACCGGATAATGTTTTTTTTGTTTGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcat                                                     >  1:3136035/1‑100 (MQ=255)
                                                 ccTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatc                                                    <  2:1340023/100‑1 (MQ=255)
                                                  cTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatca                                                   <  2:1839274/100‑1 (MQ=255)
                                                  cTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatca                                                   <  2:623001/100‑1 (MQ=255)
                                                    ttcttcACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCatcatc                                                 <  1:407745/100‑1 (MQ=255)
                                                     tcttcACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCAtcatca                                                <  2:4041671/100‑1 (MQ=255)
                                                           ccGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTggg                                          <  2:996952/100‑1 (MQ=255)
                                                             ggATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTa                                        <  1:2536975/100‑1 (MQ=255)
                                                              gATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTaa                                       <  2:3219236/100‑1 (MQ=255)
                                                                tAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATg                                     <  1:3751617/100‑1 (MQ=255)
                                                                tAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACGCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATg                                     <  1:2062010/100‑1 (MQ=255)
                                                                  aTGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGc                                   <  2:3131658/100‑1 (MQ=255)
                                                                   tGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCa                                  <  2:3305769/100‑1 (MQ=255)
                                                                     gtttttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCAtt                                >  1:3498529/2‑100 (MQ=255)
                                                                      tttttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCAttt                               >  2:2047648/1‑100 (MQ=255)
                                                                       ttttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTg                              >  2:1196689/1‑100 (MQ=255)
                                                                        tttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTg                              <  2:1288336/99‑1 (MQ=255)
                                                                           cAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACt                          <  2:100564/100‑1 (MQ=255)
                                                                             ggTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGa                        >  1:1753268/1‑100 (MQ=255)
                                                                              gTCGATAGCCTGCAACAGCGGTGGTTGCAGGGTGGCCGCTTGCTTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGAc                       <  1:3449936/100‑1 (MQ=255)
                                                                               tCGATAGTCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACg                      <  2:3151455/100‑1 (MQ=255)
                                                                                 gccacccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGcc                    >  2:3430292/6‑100 (MQ=38)
                                                                                 gATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGcc                    <  1:2979627/100‑1 (MQ=255)
                                                                                   tAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCg                  <  2:2024093/100‑1 (MQ=255)
                                                                                   cacccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGAc                       >  1:1652601/4‑95 (MQ=38)
                                                                                     gCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGt                 >  1:583858/1‑99 (MQ=255)
                                                                                     cccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTc                >  2:1597732/2‑100 (MQ=255)
                                                                                     cccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTc                >  2:4271021/2‑100 (MQ=255)
                                                                                     cccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTc                >  1:1796493/2‑100 (MQ=255)
                                                                                      ccTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCa               >  1:1952413/1‑100 (MQ=255)
                                                                                      ccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCa               >  1:2755516/1‑100 (MQ=255)
                                                                                       cTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAg              <  1:2088799/100‑1 (MQ=255)
                                                                                       cTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAg              >  2:1980076/1‑100 (MQ=255)
                                                                                        tacAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGc             <  2:799754/98‑1 (MQ=255)
                                                                                          cAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGa           <  2:2000455/100‑1 (MQ=255)
                                                                                               gccgttgttGCAGGGTGGCCTCTTGCGTGCTCAGAGTGTCGGCCTCACGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCa        <  1:2270827/95‑1 (MQ=255)
                                                                                               gCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCata      >  2:2389806/1‑100 (MQ=255)
                                                                                               gCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCa        <  1:692139/98‑1 (MQ=255)
                                                                                                 ggTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATAtt    >  1:3659166/1‑100 (MQ=255)
                                                                                                 ggTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATAtt    >  1:1581994/1‑100 (MQ=255)
                                                                                                 ggTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATAtt    >  2:2928117/1‑100 (MQ=255)
                                                                                                 ggTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATAtt    >  1:944010/1‑100 (MQ=255)
                                                                                                   ttgttgCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCa  <  1:3839301/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCA  >  NC_000913/3702740‑3702938

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: