Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,669,340 A→G 22.6% 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,3400AG22.6% 335.6 / 95.3 189intergenic (‑152/‑252)yhjA/treFputative cytochrome C peroxidase/cytoplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base A (77/68);  new base G (23/20);  total (101/88)
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 189 total aligned reads displayed.

CTGCAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCGGTA  >  NC_000913/3669242‑3669438
                                                                                                  |                                                                                                  
ctgcAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGACATACGCTTTTTGACAGGAATCAATCAATTACATAAATATAAACTATCAATAACAt                                                                                                   >  2:3271735/1‑100 (MQ=255)
  gcAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGTTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAg                                                                                                 >  2:2884629/1‑100 (MQ=255)
    aaaCGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGcc                                                                                               >  2:575911/1‑100 (MQ=255)
     aaCGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCt                                                                                              >  1:2708946/1‑100 (MQ=255)
     aaCGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCt                                                                                              >  1:1834241/1‑100 (MQ=255)
       cGTCTTTGCGAGACAAAGCCATCATAAAAAACAGGGATATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCtg                                                                                             <  2:3079461/99‑1 (MQ=255)
       cGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCtgt                                                                                            <  2:3790109/100‑1 (MQ=255)
        gTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCtgtg                                                                                           <  1:3959988/100‑1 (MQ=255)
          ctctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGcc                                                                                         <  1:3410410/100‑1 (MQ=255)
          ctctGCGAGACAAAGCAATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGcc                                                                                         <  2:1268004/100‑1 (MQ=255)
            ctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTa                                                                                       <  1:15801/100‑1 (MQ=255)
            ctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTa                                                                                       >  1:151659/1‑100 (MQ=255)
            ctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTa                                                                                       >  1:2610306/1‑100 (MQ=255)
               cGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTa                                                                                    >  1:1752470/1‑100 (MQ=255)
                 agaCAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGTTATTAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTAcc                                                                                  >  2:1099226/1‑100 (MQ=255)
                 aaaCAAAACCATCATAAAAAACAGGAAGATAAACTTTTCGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGTCTATTAc                                                                                   <  2:1021828/97‑1 (MQ=255)
                  gatctaGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTAccc                                                                                 <  1:3574128/95‑1 (MQ=255)
                  gaCAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTAccc                                                                                 <  2:3345381/100‑1 (MQ=255)
                  gaCAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTAccc                                                                                 >  1:2931138/1‑100 (MQ=255)
                   aCAAAGCCATGATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCa                                                                                <  2:2708930/100‑1 (MQ=255)
                    ctaaGCCATCATAGAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAg                                                                               >  1:2797570/3‑100 (MQ=255)
                    cAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCAAg                                                                               >  2:3330740/1‑100 (MQ=255)
                     aaaGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGCAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAg                                                                               <  2:1424778/99‑1 (MQ=255)
                     aaaGCCATCATAAAAAACAGGGAGATACGCGTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGa                                                                              <  2:2192360/100‑1 (MQ=255)
                      naGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATAAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCAAGAt                                                                             >  2:882206/2‑100 (MQ=255)
                      aaGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGAt                                                                             >  1:3600893/1‑100 (MQ=255)
                       aGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATc                                                                            <  1:324736/100‑1 (MQ=255)
                       aGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTGTTACCCAGATc                                                                            >  2:3574105/1‑100 (MQ=255)
                         ccATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAg                                                                          <  1:3412754/100‑1 (MQ=255)
                         ccATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAg                                                                          >  1:1264675/1‑100 (MQ=255)
                         caataatAAAAAACAGGGAGAAACGCTTTTTAATACGAATCAATCAATTACATGAATACAAACTATCAATAAGATAGACTGTGCCTATTACCCAGATCAg                                                                          <  2:3599104/95‑1 (MQ=255)
                          catcatAAAAAACAGGGAGATACGCTTTTTGATTTGTTTCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAgg                                                                         >  1:887583/1‑100 (MQ=255)
                           atcatAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGa                                                                        >  2:1676145/1‑100 (MQ=255)
                             catAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAAt                                                                      >  1:1166650/1‑100 (MQ=255)
                             catAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAAt                                                                      >  2:3300689/1‑100 (MQ=255)
                              atAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCACTAAGAAAGCCTGTGCCTATTACCCAGATCAGGaa                                                                       >  2:3181041/1‑98 (MQ=255)
                              atAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATg                                                                     <  1:3707941/100‑1 (MQ=255)
                              atAAAAAACAGGAAAATAAGCTTTTTGATAGTAATCAATCAATTACATGAATATTAACTCTCAATAAGATAGCCTGTGCCTATTACCCAGATc                                                                            <  2:2554877/93‑1 (MQ=255)
                               tAAAAAACAGGGAGATACGCTTTTTTTTTTTTTTTAATCAATTACATGAATATAAACTATCAATATTATAGCCTGTGCCTATTACTCAGATCAGGAATg                                                                     >  2:3550172/1‑99 (MQ=255)
                               tAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAAt                                                                      >  2:2076707/1‑98 (MQ=255)
                               tAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATgg                                                                    <  2:1752458/100‑1 (MQ=255)
                               tAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATgg                                                                    >  2:3168790/1‑100 (MQ=255)
                                 aaaaaTAGGGACATACGCATTATGATAGGAATCAATCAACTACATAAATATAAACTACCAAAAAGATACCCTGTGCCTACTACCCACATCAGGAATTGc                                                                   <  2:1280640/99‑1 (MQ=255)
                                        gggAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATg                                                           >  1:917196/1‑100 (MQ=255)
                                        gggAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATg                                                           >  2:2393300/1‑100 (MQ=255)
                                        gggAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCATGTGCCTATTACCCAGATCAGGAATGGCATATCATg                                                           >  2:3036879/1‑100 (MQ=255)
                                          gagaTACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGAt                                                         >  1:2680968/1‑100 (MQ=255)
                                            gaTACGCTTTTTGATAGGAATTAATCAATTGCATGAGTATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGAtt                                                        <  2:2782735/99‑1 (MQ=255)
                                            gaTACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTg                                                       >  2:1762061/1‑100 (MQ=255)
                                            gaTACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAACATAGCCTGTCCCTATTACCCAGATCAGGAATGGCATATCATGAtt                                                        >  2:1732525/1‑99 (MQ=255)
                                              tACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcg                                                     >  1:110986/1‑100 (MQ=255)
                                               aCGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcgg                                                    >  1:3486378/1‑100 (MQ=255)
                                                cGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcggc                                                   >  1:2973055/1‑100 (MQ=255)
                                                cGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGGTCAGGAATGGCATATCATGATTgcggc                                                   <  2:358547/100‑1 (MQ=255)
                                                cGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcggc                                                   >  1:934090/1‑100 (MQ=255)
                                                cGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcggc                                                   >  2:3930051/1‑100 (MQ=255)
                                                cGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcggc                                                   >  1:4045145/1‑100 (MQ=255)
                                                 gCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcggcg                                                  <  1:135993/100‑1 (MQ=255)
                                                  cTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCga                                                 >  1:3797593/1‑100 (MQ=255)
                                                   tttttGATAGGAATCAATCAAATACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCtag                                                <  2:1655412/100‑1 (MQ=255)
                                                    ttttGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCgaga                                               >  1:1693313/1‑100 (MQ=255)
                                                    ttttGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCgaga                                               <  1:211945/100‑1 (MQ=255)
                                                     tttGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATGACCCAGATCAGGAATGGCATATCATGATTGCGGCGagag                                              <  2:2937221/100‑1 (MQ=255)
                                                      ttGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGc                                             <  1:3689339/100‑1 (MQ=255)
                                                       tGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGCTAGCCTGTGCCGAATACTCAGATCCGGAATGGCATATCATGATTGCGGCGAGAGCa                                            >  2:1205272/1‑100 (MQ=255)
                                                        gATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGc                                             >  2:1866222/1‑98 (MQ=255)
                                                        gATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAg                                           <  1:2684542/100‑1 (MQ=255)
                                                        gATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAg                                           >  2:1163719/1‑100 (MQ=255)
                                                            ggAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGAtt                                       >  1:620391/1‑100 (MQ=255)
                                                            ggAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCATGTCCCTATTACCCAGATCAGGCATGGCACATCATGATTGCGGCGCCATCAGGa                                         >  2:2217336/1‑98 (MQ=255)
                                                             gAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTg                                      <  1:3672075/100‑1 (MQ=255)
                                                                tcaatcaaTTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAAt                                   <  2:584995/100‑1 (MQ=255)
                                                                 caatcaatTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATg                                  <  1:229958/100‑1 (MQ=255)
                                                                  aatcaatTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGt                                 >  2:4039036/1‑100 (MQ=255)
                                                                   atcaatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGtt                                >  1:804505/1‑100 (MQ=255)
                                                                    tcaatTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTCGTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTc                               <  2:2637284/100‑1 (MQ=255)
                                                                    tcaatTACATGAATATAAACTATCAATAAGATAGCATGTGCCTATTAACAAGAACAGCAATGGCATATCATGCTTGCGGCGACCGCAGGAATgcaggtaa                               >  2:3714290/1‑93 (MQ=255)
                                                                     caatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCt                              <  1:2454130/100‑1 (MQ=255)
                                                                     caatTACATGAATATAAACTATCAATAAGATAGCCNGTGCCTANTACCCAGATNAGGANTGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCt                              <  1:2754169/100‑1 (MQ=255)
                                                                       atTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtga                            <  2:507680/100‑1 (MQ=255)
                                                                           cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                        >  2:1789963/1‑100 (MQ=255)
                                                                           cATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                        >  1:3845893/1‑100 (MQ=255)
                                                                             tGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatga                      >  1:2020027/1‑100 (MQ=255)
                                                                                atatAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatgaaaa                   <  1:682541/100‑1 (MQ=255)
                                                                                atatAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatgaaaa                   <  1:858374/100‑1 (MQ=255)
                                                                                 tataAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGaaaaa                  >  1:3714279/1‑100 (MQ=255)
                                                                                  ataAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaac                 >  2:1927164/1‑100 (MQ=255)
                                                                                  ataAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaac                 <  2:1588767/100‑1 (MQ=255)
                                                                                   taAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaaca                >  1:481166/1‑100 (MQ=255)
                                                                                       cTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAacac            <  1:401138/100‑1 (MQ=255)
                                                                                       cTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAacac            >  1:3175329/1‑100 (MQ=255)
                                                                                       cTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAacac            >  2:3516314/1‑100 (MQ=255)
                                                                                        tATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACAcc           >  2:1310182/1‑100 (MQ=255)
                                                                                           cAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCaaa        >  2:1383285/1‑100 (MQ=255)
                                                                                           cAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCaaa        >  1:1400856/1‑100 (MQ=255)
                                                                                            aaTAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACACAACCaaa        >  2:2494092/1‑99 (MQ=255)
                                                                                              tAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCg     <  1:3640823/100‑1 (MQ=255)
                                                                                               aaGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCgg    <  1:3300709/100‑1 (MQ=255)
                                                                                                aGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCGGt   >  1:1102307/1‑100 (MQ=255)
                                                                                                 gATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCGGTa  >  1:251391/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
CTGCAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCGGTA  >  NC_000913/3669242‑3669438

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: