Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,669,340 A→G intergenic (‑152/‑252) yhjA ← / → treF putative cytochrome C peroxidase/cytoplasmic trehalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,669,3400AG98.6% 481.5 / ‑3.0 140intergenic (‑152/‑252)yhjA/treFputative cytochrome C peroxidase/cytoplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (58/80);  total (59/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

Only 100 of 140 total aligned reads displayed.

GCAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCGGT  >  NC_000913/3669244‑3669437
                                                                                                |                                                                                                 
gcAAACGTCTCTGCGAGACAAAGCCATCATAAATAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAg                                                                                                >  2:2642213/1‑100 (MQ=255)
gcAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAg                                                                                                >  1:1924333/1‑100 (MQ=255)
  aaaCGTCTCTGCGAGACAAAGACATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGcc                                                                                              <  1:1693445/100‑1 (MQ=255)
     cGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCtgt                                                                                           >  2:3222617/1‑100 (MQ=255)
      gTCTCTGCGAGACAAAGCCATCATAAAAAACAGGTAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCtgtg                                                                                          <  1:2192116/100‑1 (MQ=255)
      gTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCtgtg                                                                                          <  1:3219839/100‑1 (MQ=255)
       tctcgcagagaCAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGc                                                                                         >  2:3396844/8‑100 (MQ=255)
       tctcTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGc                                                                                         <  1:296374/100‑1 (MQ=255)
        ctctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGcc                                                                                        >  2:927161/1‑100 (MQ=255)
        ctctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGcc                                                                                        <  1:3046336/100‑1 (MQ=255)
        ctctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGcc                                                                                        <  1:1634801/100‑1 (MQ=255)
          ctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGcc                                                                                        <  2:1538209/98‑1 (MQ=255)
          ctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTa                                                                                      >  1:160247/1‑100 (MQ=255)
          ctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTa                                                                                      >  1:1309671/1‑100 (MQ=255)
            gCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTAtt                                                                                    >  2:1962843/1‑100 (MQ=255)
            gCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTAtt                                                                                    <  2:1755556/100‑1 (MQ=255)
             cGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTa                                                                                   >  2:3286456/1‑100 (MQ=255)
             cGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTa                                                                                   <  1:2708108/100‑1 (MQ=255)
                 aCAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCa                                                                               <  2:2240521/100‑1 (MQ=255)
                 aCAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCa                                                                               <  1:1099909/100‑1 (MQ=255)
                  cAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTAccc                                                                                <  1:1635627/98‑1 (MQ=255)
                  cAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAg                                                                              >  2:2404388/1‑100 (MQ=255)
                  cAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAg                                                                              >  2:1057228/1‑100 (MQ=255)
                  cAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAg                                                                              >  2:1232754/1‑100 (MQ=255)
                   aaaGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGa                                                                             <  1:106495/100‑1 (MQ=255)
                     aGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATc                                                                           >  2:243534/1‑100 (MQ=255)
                      gacatcatAAAAAACAGGGAGATACGCTTTTTGATAGTAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATc                                                                           <  1:2536176/97‑1 (MQ=255)
                        cctatcaaaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAgg                                                                        <  2:1937075/94‑1 (MQ=255)
                        catcatAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAgg                                                                        >  2:667264/1‑100 (MQ=255)
                         atcatAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGa                                                                       >  1:3252721/1‑100 (MQ=255)
                         atcatAAAAAACAGAACAATACGCTTTCAGAGATGAAACAATCAATTACCTTAATATAAACTATCAATAACGTAGCCTGTGCCTATAAACCAGAACAgg                                                                        <  1:1546637/99‑1 (MQ=255)
                           catAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGaan                                                                     <  2:2644945/100‑2 (MQ=255)
                            atAAAAAACAGGGAGATACGCTTTTTTTTTTTTTTTTATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCTTATTACGCATATCAGGAAt                                                                     >  1:445407/1‑99 (MQ=255)
                            atAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATg                                                                    <  2:1821249/100‑1 (MQ=255)
                            atAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATg                                                                    >  2:712484/1‑100 (MQ=255)
                             tAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATgg                                                                   <  2:1865859/100‑1 (MQ=255)
                             tAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATgg                                                                   <  1:3454969/100‑1 (MQ=255)
                              aaaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGc                                                                  <  2:2016716/100‑1 (MQ=255)
                              aaaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGc                                                                  <  1:2509102/100‑1 (MQ=255)
                               aaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCa                                                                 <  1:927149/100‑1 (MQ=255)
                                  aaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCatat                                                              <  2:1309690/100‑1 (MQ=255)
                                   aCAGGGAGACACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATc                                                             <  2:3312983/100‑1 (MQ=255)
                                    cAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCa                                                            >  1:1112853/1‑100 (MQ=255)
                                    cAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCa                                                            >  2:848480/1‑100 (MQ=255)
                                        gagaTACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGAt                                                        <  2:3330257/100‑1 (MQ=255)
                                          gaTACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTg                                                      >  1:207230/1‑100 (MQ=255)
                                             aCGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcgg                                                   <  1:221055/100‑1 (MQ=255)
                                              cGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcggc                                                  <  2:328211/100‑1 (MQ=255)
                                                  ttttGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCgaga                                              >  1:2477313/1‑100 (MQ=255)
                                                     tGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCa                                           <  2:3200639/100‑1 (MQ=255)
                                                     tGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCa                                           >  1:3194023/1‑100 (MQ=255)
                                                      gATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAg                                          >  2:2694769/1‑100 (MQ=255)
                                                      gATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAg                                          >  1:2559394/1‑100 (MQ=255)
                                                      gATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAg                                          >  2:3436931/1‑100 (MQ=255)
                                                         aGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGAt                                       >  1:1848933/1‑100 (MQ=255)
                                                           gAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTg                                     <  1:718820/100‑1 (MQ=255)
                                                           gAATCAATCAATTACATGAATATAAACTATCAATAAGGGAGCCTGTGCCTATTACCCACATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTg                                     <  1:1894460/100‑1 (MQ=255)
                                                            aatcaatcAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGa                                    >  1:1236677/1‑100 (MQ=255)
                                                              tcaatcaaTTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAAt                                  <  2:682345/100‑1 (MQ=255)
                                                              tcaatcaaTTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAAt                                  <  2:2253068/100‑1 (MQ=255)
                                                              tcaatcaaTTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAAt                                  <  2:2217316/100‑1 (MQ=255)
                                                                aatcaatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGt                                >  2:26219/1‑100 (MQ=255)
                                                                  tcaatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTc                              <  2:57367/100‑1 (MQ=255)
                                                                  tcaatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTc                              <  1:3273620/100‑1 (MQ=255)
                                                                  tcaatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTc                              <  2:2576528/100‑1 (MQ=255)
                                                                  tcaatTACATAAATATAAACAATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGTGCAGGATTGAATGTTc                              <  1:2275723/100‑1 (MQ=255)
                                                                   caatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCt                             <  1:3396792/100‑1 (MQ=255)
                                                                    aatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtg                            >  2:2463242/1‑100 (MQ=255)
                                                                      ttACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaa                          <  2:2160304/100‑1 (MQ=255)
                                                                        aCATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaa                        <  2:1280294/100‑1 (MQ=255)
                                                                         cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                       >  1:2853415/1‑100 (MQ=255)
                                                                         cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                       <  1:1057214/100‑1 (MQ=255)
                                                                         cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                       <  2:137135/100‑1 (MQ=255)
                                                                         cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                       >  2:2857059/1‑100 (MQ=255)
                                                                         cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                       >  2:2700245/1‑100 (MQ=255)
                                                                         cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                       >  1:2244869/1‑100 (MQ=255)
                                                                         cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATAACATGATAGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                       >  1:2442060/1‑100 (MQ=255)
                                                                          aTGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatg                      <  1:1937046/100‑1 (MQ=255)
                                                                          aTGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatg                      <  1:292568/100‑1 (MQ=255)
                                                                            gAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatgaa                    >  2:3414370/1‑100 (MQ=255)
                                                                            gAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatgaa                    >  1:3058439/1‑100 (MQ=255)
                                                                            gAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGAACAAAAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatga                     >  1:1311031/1‑99 (MQ=255)
                                                                             aaTATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatgaaa                   <  2:3359684/100‑1 (MQ=255)
                                                                                 taAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaaca               <  1:2685037/100‑1 (MQ=255)
                                                                                 taAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaaca               >  2:1339475/1‑100 (MQ=255)
                                                                                 taAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaaca               <  2:1949015/100‑1 (MQ=255)
                                                                                   aaCTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaacaac             <  1:726587/100‑1 (MQ=255)
                                                                                   aaCTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaacaac             >  2:2207737/1‑100 (MQ=255)
                                                                                   aaCTATCAATAAGGTAGCCAGTGCCTATGACCCAGATCAGGAAAGGCATATAATGATTGAGGCGAGAGCAGGATTGAAAGTTCTGAAAATGAAAaacatc             <  1:987420/100‑3 (MQ=255)
                                                                                    aCTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAacaaca            <  1:2207705/100‑1 (MQ=255)
                                                                                    aCTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAacaaca            >  1:3097059/1‑100 (MQ=255)
                                                                                      tATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACAcc          <  2:2984065/100‑1 (MQ=255)
                                                                                      tATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACAcc          <  2:291628/100‑1 (MQ=255)
                                                                                       aTCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCa         <  1:1775058/100‑1 (MQ=255)
                                                                                         cAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCaaa       >  2:1583041/1‑100 (MQ=255)
                                                                                          aaTAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAAt      <  1:1070042/100‑1 (MQ=255)
                                                                                          aaTAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAAt      <  1:2097859/100‑1 (MQ=255)
                                                                                           aTAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATc     <  2:2617906/100‑1 (MQ=255)
                                                                                             aaGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCgg   <  2:529245/100‑1 (MQ=255)
                                                                                              aGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCGGt  <  2:207233/100‑1 (MQ=255)
                                                                                                |                                                                                                 
GCAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCGGT  >  NC_000913/3669244‑3669437

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: