Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,002,325 C→T 78.9% I300I (ATC→ATT corA → magnesium/nickel/cobalt transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,002,3250CT78.9% 193.5 / 49.0 104I300I (ATC→ATTcorAmagnesium/nickel/cobalt transporter
Reads supporting (aligned to +/- strand):  ref base C (10/12);  new base T (33/49);  total (43/61)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 106 total aligned reads displayed.

GCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAG  >  NC_000913/4002226‑4002424
                                                                                                   |                                                                                                   
gccgccAACACTGGTTGCTTCCAGCTACGGGATGAACTTTGAGTTTATGCCGGAACTGAAGTGGAGTTTCGGCTACCCTGCCGCGATTATCTTTATGATc                                                                                                     >  2:1876356/1‑100 (MQ=255)
  cgccgACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCt                                                                                                   <  1:1204526/100‑1 (MQ=255)
     cgACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTcgc                                                                                                >  1:2489566/1‑100 (MQ=255)
      gACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTcgcg                                                                                               >  1:2190284/1‑100 (MQ=255)
         gCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTAGCCTGGCGCGATTATCTTTATGATTCTCGCGGGc                                                                                            <  1:1710220/100‑1 (MQ=255)
         gCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGc                                                                                            <  2:517479/100‑1 (MQ=255)
            cGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTg                                                                                         <  2:2264390/100‑1 (MQ=255)
            cGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTg                                                                                         >  1:2357684/1‑100 (MQ=255)
             gTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTgg                                                                                        <  2:233338/100‑1 (MQ=255)
             gTTGCTTCCAGCGATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTgg                                                                                        <  2:2024711/100‑1 (MQ=255)
                gCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCAc                                                                                     >  2:1906056/1‑100 (MQ=255)
                gCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCAc                                                                                     >  1:7358/1‑100 (MQ=255)
                 cTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCAcc                                                                                    >  2:1002656/1‑100 (MQ=255)
                 cTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCAcc                                                                                    >  1:2098058/1‑100 (MQ=255)
                  ttCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAATTGGAGCTTCGGCTACCCTGGCGCGATTATTTTTATGATTCTCGCGGGCCTGGAACCg                                                                                   <  1:1091392/100‑1 (MQ=255)
                  ttCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGATTTTCGGTTACCCTGGCGCGATTTTCTTTATGATTCTCGCGGGACTGGCACCg                                                                                   <  1:327330/100‑1 (MQ=255)
                  ttCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCg                                                                                   <  2:1981145/100‑1 (MQ=255)
                  ttCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCg                                                                                   <  2:2489585/100‑1 (MQ=255)
                      aGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATc                                                                               >  1:1110340/1‑100 (MQ=255)
                       gttctGGCATAAACTCAAAGTTCATGCCAGAACTGAAGTGGCGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCt                                                                              >  1:419263/5‑100 (MQ=25)
                        cTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTg                                                                             <  1:1875087/100‑1 (MQ=255)
                        cTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTg                                                                             >  1:1910627/1‑100 (MQ=255)
                          aTGGCATGAACTTTGAGTTTATGTCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTa                                                                           >  2:398479/1‑100 (MQ=255)
                           tGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTAc                                                                          >  2:2316838/1‑100 (MQ=255)
                             gCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACtt                                                                        <  2:1258817/100‑1 (MQ=255)
                               aTGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTa                                                                      <  2:633471/100‑1 (MQ=255)
                                tGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTaa                                                                     >  1:59991/1‑100 (MQ=255)
                                tGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTaa                                                                     <  2:1330438/100‑1 (MQ=255)
                                 gAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAg                                                                    >  1:1214314/1‑100 (MQ=255)
                                    cTTTGAGTCTATGCCAGAACTGAAGTGGAGCTTCCGCTACCCTGGCGCGATTATCTTTATGAGTCACGCGGGCCTGGCACCGTATCTGTACTTTAAGCgg                                                                 <  1:1630736/100‑1 (MQ=255)
                                       tGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGgaag                                                              <  2:928956/100‑1 (MQ=255)
                                       tGAGTTTATGCCAGAACTGAAGTGGAGCCTCGGCTACCCTGGCGCGATTATCTTTATGATCCGCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGgaag                                                              <  1:440146/100‑1 (MQ=255)
                                         aGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGgaagaa                                                            >  1:1527121/1‑100 (MQ=255)
                                             tATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTgg                                                        >  1:992059/1‑100 (MQ=255)
                                             tATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTgg                                                        >  1:276914/1‑100 (MQ=255)
                                              aTGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGt                                                       <  2:360518/100‑1 (MQ=255)
                                               tGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGtt                                                      >  2:1964760/1‑100 (MQ=255)
                                               tGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGtt                                                      <  2:2337628/100‑1 (MQ=255)
                                                gCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTCTGAGTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTg                                                     <  1:527392/100‑1 (MQ=255)
                                                gCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTg                                                     <  2:1304801/100‑1 (MQ=255)
                                                gCCAGAACTGAAGTGGAGCTTCGGCTACCCCGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTg                                                     <  2:553066/100‑1 (MQ=255)
                                                  cAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTa                                                   <  2:1969251/100‑1 (MQ=255)
                                                   aGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTaa                                                  <  2:1001542/100‑1 (MQ=255)
                                                    gAACTGCAGTGGAGCTTCGGCTCCCCTGCCTCGATCATCTTTATGATTCCCGCGGGCCTGGCACAATATCTTTACTTTAAGCGGAAGAACTGGTTGTaaa                                                 <  1:2422769/100‑1 (MQ=255)
                                                    gAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTaaa                                                 <  2:1358968/100‑1 (MQ=255)
                                                      aCTGAAGTTGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTaaaaa                                               <  2:853975/100‑1 (MQ=255)
                                                      aCTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTaaaaa                                               >  2:1402948/1‑100 (MQ=255)
                                                      aCTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTaaaaa                                               >  2:582976/1‑100 (MQ=255)
                                                       cttccgatctgCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTaaaaaa                                              <  1:2299096/90‑1 (MQ=255)
                                                        tGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAAc                                             <  2:1538295/100‑1 (MQ=255)
                                                        tGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAAc                                             <  1:590483/100‑1 (MQ=255)
                                                        tGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTAGTTGAAAAAAAc                                             <  1:1962331/100‑1 (MQ=255)
                                                         ggaGTGGTGCTTCGGATATGCTGGCGGGATTAGCTGTTTGATTGTCGAGTGCCTGGGACCGTATCTGTACTTTAGGCGGAAGAACTGGTTGTaaa                                                 <  1:2404605/93‑1 (MQ=255)
                                                            gTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACgaga                                         >  2:192302/1‑100 (MQ=255)
                                                            gTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACgaga                                         >  1:22021/1‑100 (MQ=255)
                                                              ggAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGc                                       <  2:1472844/100‑1 (MQ=255)
                                                              ggAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGc                                       <  2:1472966/100‑1 (MQ=255)
                                                                aGCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGATGAACTGGTTGTAAAAAACGACAGCgg                                     <  1:1899417/100‑1 (MQ=255)
                                                                 gCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTaaaaaaagatcggaag                                    >  2:2299110/1‑90 (MQ=255)
                                                                 gCTTCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGt                                    >  1:1972685/1‑100 (MQ=255)
                                                                    tCGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGc                                 <  1:1885253/100‑1 (MQ=255)
                                                                     cGGCTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCt                                <  1:2165169/100‑1 (MQ=255)
                                                                       gCTACCGTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTAGCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCtt                               <  1:1352075/99‑1 (MQ=255)
                                                                        cTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAg                             >  2:2578333/1‑100 (MQ=255)
                                                                        cTACCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAg                             >  2:1202563/1‑100 (MQ=255)
                                                                          aCCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGtn                           <  1:1095539/100‑2 (MQ=255)
                                                                          aCCCTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTc                           >  2:764041/1‑100 (MQ=255)
                                                                            ccTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTg                         >  1:75084/1‑100 (MQ=255)
                                                                             cTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTgg                        >  1:52694/1‑100 (MQ=255)
                                                                             cTGGCGCGATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTgg                        >  1:54835/1‑100 (MQ=255)
                                                                                  gcgATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAg                   <  2:1910637/100‑1 (MQ=255)
                                                                                  gcgATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAg                   <  2:2480404/100‑1 (MQ=255)
                                                                                  gcgATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAg                   >  1:2104766/1‑100 (MQ=255)
                                                                                   cgATTATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGTTGGCTTAGTCTGGCTAAGc                  >  1:2221438/1‑100 (MQ=255)
                                                                                      ttATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCAc               <  2:1043993/100‑1 (MQ=255)
                                                                                      ttATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAc               <  2:956801/100‑1 (MQ=255)
                                                                                      ttATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAc               <  2:1898292/100‑1 (MQ=255)
                                                                                      ttATCTTTATGATCCTCGAGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAc               <  1:2021135/100‑1 (MQ=255)
                                                                                      ttAACTTTATGATTCTCGAGAGCCTGGCACTGTATCTGGACTTTAAGCGGAAAAAAAAAAAAAAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAc               <  1:90457/100‑1 (MQ=255)
                                                                                       tATCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCAcc              <  2:433738/100‑1 (MQ=255)
                                                                                        aTCTTTATGATTCTCGCGGGCCTGGCAACGTATCTGTACTTTAAGCGGAAGAACTGGTTGTCAAACACGAGAGCGGTCGCTTAGTCTGGCTGAGCCActt             >  1:546783/1‑98 (MQ=255)
                                                                                         tCTTTGTGATTCTCCCGGGCCTGGCACCGTATCTGTACTTTGAGCGGAAGATCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCtn            <  1:2018116/100‑2 (MQ=255)
                                                                                         tCTTTATGCTCCTCGCGCGCCTGGCACCATATCTGTACTTTATGCGGAAGAACTGGTTGTAAAAAACGAGATCGGTGGCTTAGTCTGGCTAAGCCACCtc            <  1:574201/100‑2 (MQ=255)
                                                                                         tCTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTg            >  1:1407832/1‑100 (MQ=255)
                                                                                         tCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTg            <  2:1480065/100‑1 (MQ=255)
                                                                                         tCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAccgc            <  2:1105057/100‑3 (MQ=38)
                                                                                          cTTTATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGt           <  1:2316824/100‑1 (MQ=255)
                                                                                          cTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGt           >  1:310541/1‑100 (MQ=255)
                                                                                           tttATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGtt          >  1:802850/1‑100 (MQ=255)
                                                                                            ttATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTa         >  2:2508159/1‑100 (MQ=255)
                                                                                            ttATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTa         >  2:1958498/1‑100 (MQ=255)
                                                                                            ttATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTa         <  2:1668794/100‑1 (MQ=255)
                                                                                             tATGATTCTCGCGGGCCTGGTACCGTATCTGTACTGTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTAt        <  1:2412282/100‑1 (MQ=255)
                                                                                             tATGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTAt        >  2:2103185/1‑100 (MQ=255)
                                                                                              aTGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTAtt       <  2:408498/100‑1 (MQ=255)
                                                                                               tGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTc      <  1:110539/100‑1 (MQ=255)
                                                                                               tGATTCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTc      >  1:562053/1‑100 (MQ=255)
                                                                                               tGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTc      >  2:2292527/1‑100 (MQ=255)
                                                                                                gATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCa     <  2:1214318/100‑1 (MQ=255)
                                                                                                   ccTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGGAAGACTGGTTGTAAAAAACGAGAGCTGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAg  <  1:779885/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAG  >  NC_000913/4002226‑4002424

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: