Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,669,340 A→G 100% 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,3400AG99.4% 377.3 / ‑6.3 111intergenic (‑152/‑252)yhjA/treFputative cytochrome C peroxidase/cytoplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base G (55/53);  minor base T (1/1);  total (56/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

GCAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCGG  >  NC_000913/3669244‑3669436
                                                                                                |                                                                                                
gcAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAg                                                                                               <  1:216047/100‑1 (MQ=255)
gcAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAg                                                                                               >  1:1999737/1‑100 (MQ=255)
 cAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGc                                                                                              <  1:2303594/100‑1 (MQ=255)
   aaCGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCt                                                                                            <  1:1606036/100‑1 (MQ=255)
   aaCGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCt                                                                                            >  1:833602/1‑100 (MQ=255)
   aaCGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCt                                                                                            >  2:550943/1‑100 (MQ=255)
     tgtctcgcagagaCAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCtgt                                                                                          >  1:2091258/10‑100 (MQ=255)
     tgtctcgcagagaCAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCtgt                                                                                          >  1:339987/10‑100 (MQ=255)
     cGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCtgt                                                                                          >  1:1316011/1‑100 (MQ=255)
      gTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCACTTACATGAATATAAACTATCAATAAGGTAGCCtgtg                                                                                         <  1:1811434/100‑1 (MQ=255)
        ctctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGcc                                                                                       >  1:300524/1‑100 (MQ=255)
        ctctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGc                                                                                        <  1:1241729/99‑1 (MQ=255)
         tgttcGAGACAAATCCATAATAAAAAACAGTGATATATGCTTTCTAATCGGAATCAATCAATTACACGAATATAAACTATCAATAAGGTAGCCTGTGcc                                                                                       <  2:1975579/95‑1 (MQ=255)
          ctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTa                                                                                     >  1:450463/1‑100 (MQ=255)
          ctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTa                                                                                     >  2:797806/1‑100 (MQ=255)
           tGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGcc                                                                                       >  2:1426076/1‑97 (MQ=255)
             cGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCACTCAATTACATGCATATAAACTATCAATAAGTTAGCCTGTCTCTATTa                                                                                  >  2:72307/1‑100 (MQ=255)
             cGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTa                                                                                  <  2:2446339/100‑1 (MQ=255)
              gagaCAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTa                                                                                  <  2:2383145/99‑1 (MQ=255)
              gagaCAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTAc                                                                                 <  2:1671745/100‑1 (MQ=255)
              gagaCAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTAc                                                                                 <  1:962399/100‑1 (MQ=255)
                gaCAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTAccc                                                                               >  1:1572651/1‑100 (MQ=255)
                gaCAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTAccc                                                                               >  1:2398587/1‑100 (MQ=255)
                  cAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTAc                                                                                 <  2:1671752/96‑1 (MQ=255)
                  cAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAg                                                                             >  1:562860/1‑100 (MQ=255)
                  cAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAg                                                                             >  1:571686/1‑100 (MQ=255)
                  cAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAg                                                                             <  1:2127400/100‑1 (MQ=255)
                   aaaGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGa                                                                            >  1:1645296/1‑100 (MQ=255)
                    aaGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGCAGCCTGTGCCTATTACCCAGAt                                                                           >  1:1104344/1‑100 (MQ=255)
                     aaccATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACTCAGATc                                                                          >  2:832108/3‑100 (MQ=255)
                      gCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCa                                                                         >  2:823293/1‑100 (MQ=255)
                        catcatAAAAAACAGGGAGTTACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAgg                                                                       <  2:739238/100‑1 (MQ=255)
                        catcatAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAgg                                                                       <  2:27066/100‑1 (MQ=255)
                         atcatAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAg                                                                        <  2:2113472/98‑1 (MQ=255)
                           catAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGaa                                                                     >  2:484438/1‑99 (MQ=255)
                           catAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAAt                                                                    >  1:722328/1‑100 (MQ=255)
                            atAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATg                                                                   <  1:1568162/100‑1 (MQ=255)
                            atAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATg                                                                   >  2:402598/1‑100 (MQ=255)
                             tAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATgg                                                                  <  1:1601405/100‑1 (MQ=255)
                             tAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATgg                                                                  >  2:2136397/1‑100 (MQ=255)
                             tAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATg                                                                   >  1:2202157/1‑99 (MQ=255)
                              aaaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGc                                                                 >  1:1143262/1‑100 (MQ=255)
                               aaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCa                                                                <  2:4132/100‑1 (MQ=255)
                                aaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCat                                                               <  1:293210/100‑1 (MQ=255)
                                aaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCat                                                               <  1:287758/100‑1 (MQ=255)
                                  aaCAGGGAGATACGCTTTTTGATAGGAATCAATGAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGCAATGGCatat                                                             >  2:1685167/1‑100 (MQ=255)
                                   aCAGGGAGATACGCTTTTTTTTTTTTTTTAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATc                                                            >  2:72423/1‑100 (MQ=255)
                                    catggAGATACGCTTTTTGATAGGAATCAATCAATTACATTAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACTCAGATCAGGAATGGCATATc                                                            <  2:609185/96‑1 (MQ=255)
                                          gaTACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGAtt                                                      >  2:2271407/1‑99 (MQ=255)
                                          gaTACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATAAGGAATGGCATATCATGATTg                                                     >  2:2467326/1‑100 (MQ=255)
                                            tACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcg                                                   >  2:1690992/1‑100 (MQ=255)
                                             aCGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcgg                                                  >  1:845121/1‑100 (MQ=255)
                                              cGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcggc                                                 <  1:2136411/100‑1 (MQ=255)
                                               gCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcggcg                                                <  1:123722/100‑1 (MQ=255)
                                               gCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcggcg                                                <  2:735998/100‑1 (MQ=255)
                                                  ttttGATAGGAATCAATCAATTACATGAATCTAAACGATCAATAAGGTAGCCTGTGCCTATTACCCAGATTAGGAATGGCATATCATGATTGCGGCgag                                              <  1:1821052/99‑1 (MQ=255)
                                                  ttttGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCgaga                                             <  2:2279423/100‑1 (MQ=255)
                                                  ttttGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCgaga                                             <  1:141044/100‑1 (MQ=255)
                                                     tGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCa                                          >  2:1052237/1‑100 (MQ=255)
                                                     tGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCa                                          >  1:1790251/1‑100 (MQ=255)
                                                      gATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGTCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAg                                         >  2:1451518/1‑100 (MQ=255)
                                                      gATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAg                                         >  2:2085766/1‑100 (MQ=255)
                                                        tAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGa                                       >  2:2540907/1‑100 (MQ=255)
                                                        tAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGa                                       <  1:2426145/100‑1 (MQ=255)
                                                         agtaatcaatcAATTACATAAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGAt                                      <  2:1857023/97‑1 (MQ=255)
                                                          ggAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGAtt                                     >  1:1088015/1‑100 (MQ=255)
                                                          ggAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGAtt                                     >  1:2263851/1‑100 (MQ=255)
                                                            aatcaatcAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGa                                   >  2:1426125/1‑100 (MQ=255)
                                                            aatcaatcAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGa                                   <  1:1665604/100‑1 (MQ=255)
                                                            aatcaatcAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGa                                   >  1:1716357/1‑100 (MQ=255)
                                                            aatcaatcAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGa                                   >  1:2059436/1‑100 (MQ=255)
                                                            aatcaatcAATTACATGAATATAAACTATCAATAAGGTAGCCTGTACCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGa                                   >  2:405441/1‑100 (MQ=255)
                                                              tctatcaatcaCATGAATATAAACTATCAATAAGTTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCTGGATTGaa                                  <  2:1077744/96‑1 (MQ=255)
                                                              tcaatcaaTTACATGAATTTAAACTATCAATAAGGTAGCCTGTTCCAATTAACCAGATCAGGAATGGAATATCATGATTGCGGCGAGAGCAGGATTGAAt                                 >  2:372385/1‑100 (MQ=255)
                                                              tcaatcaaTTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAAt                                 <  1:134623/100‑1 (MQ=255)
                                                               caatcaatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATAAGGAATGGCATATCATGATTGCGGCGAGAGCAGGACTGAATg                                >  2:746048/1‑100 (MQ=255)
                                                                 atcaatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGtt                              <  1:1852192/100‑1 (MQ=255)
                                                                       tACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaa                        >  2:238907/1‑100 (MQ=255)
                                                                       tACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaa                        <  1:1924390/100‑1 (MQ=255)
                                                                       gaCAGGAAAATAAACGATCAATAAGGGAGACTTTGCCTATTACCCAGATCGGGAATGGCATCTCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaa                        <  2:2392073/99‑1 (MQ=255)
                                                                         cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                      >  2:854289/1‑100 (MQ=255)
                                                                         cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                      >  2:274407/1‑100 (MQ=255)
                                                                         cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                      <  1:1245074/100‑1 (MQ=255)
                                                                         cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                      >  1:991395/1‑100 (MQ=255)
                                                                         cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTAATACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                      >  2:1365798/1‑100 (MQ=255)
                                                                          aTGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatg                     <  1:275475/100‑1 (MQ=255)
                                                                          aTGAATATAAACTATCAATAAGGTAGCCAGTGCCTATTACCCAGATCAGAAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaa                        >  2:2147141/1‑97 (MQ=255)
                                                                            gAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatgaa                   <  1:1117216/100‑1 (MQ=255)
                                                                               tataAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGaaaaa                <  1:990296/100‑1 (MQ=255)
                                                                                  aaaCTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaacaa             <  1:238909/100‑1 (MQ=255)
                                                                                   cgatctCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaacaac            <  1:575441/95‑1 (MQ=255)
                                                                                   aaCTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGCGCAGGATTGAATGTTCTGAAAATGAAAaacaac            <  2:1533355/100‑1 (MQ=255)
                                                                                   aaCTATCAATAAGGTAGCCTGTGCCTATTACACAGATCAGGAATGGCATATCATGACTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaacaac            <  2:296809/100‑1 (MQ=255)
                                                                                    aCTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAGTGAAAAacaaca           <  1:1688217/100‑1 (MQ=255)
                                                                                     cTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAacac          <  2:1844507/100‑1 (MQ=255)
                                                                                       aTCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCa        <  1:1576139/100‑1 (MQ=255)
                                                                                         cAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAacaacagatcg      >  2:575434/1‑95 (MQ=255)
                                                                                          aaTAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCaaa      >  2:450953/1‑99 (MQ=255)
                                                                                           aTAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATc    <  1:1673318/100‑1 (MQ=255)
                                                                                             aaGGTAGCCTGTGCATATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTAAAAAAAAAAAAAAACAACACCAAATCgg  <  2:1464584/100‑1 (MQ=255)
                                                                                                |                                                                                                
GCAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCGG  >  NC_000913/3669244‑3669436

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: