Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,315,705 T→C 21.6% K22K (AAA→AAG yciF ← putative rubrerythrin/ferritin‑like metal‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,315,7050TC21.6% 308.0 / 85.1 166K22K (AAA→AAGyciFputative rubrerythrin/ferritin‑like metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (61/69);  new base C (20/16);  total (81/85)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

Only 100 of 173 total aligned reads displayed.

ATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATAC  >  NC_000913/1315607‑1315803
                                                                                                  |                                                                                                  
aTCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGtt                                                                                                   >  2:3196219/1‑100 (MQ=255)
 tCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGttt                                                                                                  <  2:1063229/100‑1 (MQ=255)
  cTGTCCATGAGTTTCCTCGAGGTGCGCATGAACAGCCTGACGTAATTTTTCATTTGATGTTGCTCTTGCGCGTTTTGCCAGTGCCCGGGTTAATTGtttt                                                                                                 >  2:2278531/1‑100 (MQ=255)
  cTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGtttt                                                                                                 <  1:3581411/100‑1 (MQ=255)
  cTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGtttt                                                                                                 <  1:1033191/100‑1 (MQ=255)
   tGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGttttt                                                                                                <  2:4002914/100‑1 (MQ=255)
      ccATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTg                                                                                             >  1:2716358/1‑100 (MQ=255)
      ccATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTg                                                                                             >  2:4164129/1‑100 (MQ=255)
      ccATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTg                                                                                             >  2:3702920/1‑100 (MQ=255)
        aTGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGTTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTgcg                                                                                           <  2:2342665/100‑1 (MQ=255)
        aTGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTgcg                                                                                           <  2:4072799/100‑1 (MQ=255)
           aGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTg                                                                                        <  2:3345080/100‑1 (MQ=255)
            gTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGt                                                                                       >  1:1073221/1‑100 (MQ=255)
            gTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGt                                                                                       >  1:1278768/1‑100 (MQ=255)
              ttCCTCGGGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAg                                                                                     <  1:1223219/100‑1 (MQ=255)
              ttCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAg                                                                                     <  1:3790621/100‑1 (MQ=255)
              ttCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAg                                                                                     <  2:1160279/100‑1 (MQ=255)
                ccTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGt                                                                                   >  1:1583138/1‑100 (MQ=255)
                ccTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGt                                                                                   >  2:578641/1‑100 (MQ=255)
                ccTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGt                                                                                   >  1:3896754/1‑100 (MQ=255)
                 cTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCATTGCCCGGGTTAATTGTTTTTCTGCGTTGTAGGTa                                                                                  >  2:1965197/1‑100 (MQ=255)
                 cTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTa                                                                                  >  2:1201240/1‑100 (MQ=255)
                 cTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTa                                                                                  >  1:747630/1‑100 (MQ=255)
                   cGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATc                                                                                <  1:441150/100‑1 (MQ=255)
                    gAGGTGCGCCTGAAACGCCTGCCTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATc                                                                                <  2:3476709/99‑1 (MQ=255)
                     aGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTg                                                                              <  1:3124134/100‑1 (MQ=255)
                      ggTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGa                                                                             <  1:962992/100‑1 (MQ=255)
                       gTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGaa                                                                            <  1:3980407/100‑1 (MQ=255)
                        tGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTg                                                                              >  2:3070621/1‑97 (MQ=255)
                         gcgcATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAg                                                                          <  1:3742845/100‑1 (MQ=255)
                         gcgcATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAg                                                                          <  1:23720/100‑1 (MQ=255)
                           gcATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGCTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCa                                                                        <  1:3869113/100‑1 (MQ=255)
                            cATGAAAAGCCTGACTTAATTTTTTTTTTTTTTTTTTTCTTGCGAGTTTTACCAGTACTTGTTTTTATTGTTTTTTTGCGCTGTAGGTATCTGAAACCa                                                                        >  2:912274/1‑99 (MQ=255)
                            cATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAg                                                                       >  2:2088642/1‑100 (MQ=255)
                            cATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAg                                                                       <  1:3166181/100‑1 (MQ=255)
                             aTGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAgg                                                                      <  1:3206212/100‑1 (MQ=255)
                              tGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGTGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGt                                                                     >  2:1325327/1‑100 (MQ=255)
                               gAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTg                                                                    >  1:3344808/1‑100 (MQ=255)
                                   aGCCTGACTTAATTTTTCATTTGCTGTTGCTCTTGCGAGATTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAACCAGGTGaat                                                                 >  2:1548981/1‑99 (MQ=255)
                                    gCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataa                                                               >  1:3111758/1‑100 (MQ=255)
                                    gCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaata                                                                >  2:3512190/1‑99 (MQ=255)
                                    gCCTGACTTAATTTGTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataa                                                               <  2:33901/100‑1 (MQ=255)
                                     ccTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaa                                                              <  1:3238472/100‑1 (MQ=255)
                                     ccTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaa                                                              >  1:3530261/1‑100 (MQ=255)
                                      cTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaat                                                             <  2:3117770/100‑1 (MQ=255)
                                      cTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaat                                                             >  1:1658018/1‑100 (MQ=255)
                                       tGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaata                                                            >  1:454002/1‑100 (MQ=255)
                                        gACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaata                                                            <  2:3116700/99‑1 (MQ=255)
                                        gACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATAc                                                           >  2:1724978/1‑100 (MQ=255)
                                         aCTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACa                                                          >  1:2910208/1‑100 (MQ=255)
                                         aCTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGATGTAGGTATCTGAAAGCAGGTGAATAAATAc                                                           >  2:1037723/1‑99 (MQ=255)
                                         aCTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACa                                                          >  2:2776432/1‑100 (MQ=255)
                                          cTTAATTTTTTTTTTTTTTTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATATATACAt                                                         >  2:712199/1‑100 (MQ=255)
                                            tAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCt                                                       <  1:3632728/100‑1 (MQ=255)
                                              attataaaattGATGTTGCTCTTGCGAGTTATGCCAGTGACAGGATTAATTGTTTTTCTTCGCTGTAGGTATCTGAAAGCAGGTTAATAAATACATCTTc                                                     <  2:768402/91‑1 (MQ=255)
                                              aTTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTc                                                     >  1:661771/1‑100 (MQ=255)
                                                ttttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCaa                                                   >  2:591917/1‑100 (MQ=255)
                                                  ttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTAGCTGAAAGCAGGTGAATAAATACATCTTCAATg                                                 <  2:2005341/100‑1 (MQ=255)
                                                   tCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATATATCTTCAATgg                                                <  2:2175715/100‑1 (MQ=255)
                                                   tCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATgg                                                <  1:5116/100‑1 (MQ=255)
                                                     aTTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTc                                              >  1:998699/1‑100 (MQ=255)
                                                     aTTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTc                                              >  1:1803281/1‑100 (MQ=255)
                                                     aTTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCACGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTc                                              >  2:2402781/1‑100 (MQ=255)
                                                      tttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCt                                             <  2:437782/100‑1 (MQ=255)
                                                        tGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttc                                           >  1:2975809/1‑100 (MQ=255)
                                                           tatttctctTGCGAGTTCTGCCAGTGCCCGGGTTCATTGCTTTTCTGCGCTGTAGGTTTCTGAATGCAGGTGAATAAATACATCTTCAATGGTCttca                                          <  2:2032843/93‑1 (MQ=255)
                                                           tGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcata                                        <  2:2952922/100‑1 (MQ=255)
                                                              tGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattc                                     <  2:1800408/100‑1 (MQ=255)
                                                                ctctTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattcat                                   >  1:2670710/1‑100 (MQ=255)
                                                                ctctTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattcat                                   >  2:430540/1‑100 (MQ=255)
                                                                   ttGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTtcatattcatat                                 <  2:379661/99‑1 (MQ=255)
                                                                   ttGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTcatattcatatt                                <  1:81835/100‑1 (MQ=255)
                                                                    tGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTATTCatatttatattt                               >  1:4139943/1‑100 (MQ=255)
                                                                     gCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATAtttt                              >  2:4158648/1‑100 (MQ=255)
                                                                         gTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTtctc                          >  1:4088941/1‑100 (MQ=255)
                                                                          ttttGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTcc                         <  2:1979757/100‑1 (MQ=255)
                                                                          ttttGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTcc                         <  1:1462114/100‑1 (MQ=255)
                                                                           tttGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTcc                         >  2:782932/1‑99 (MQ=255)
                                                                            ttGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAg                       <  2:3940669/100‑1 (MQ=255)
                                                                             tGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGt                      <  1:2821977/100‑1 (MQ=255)
                                                                             tGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGt                      <  1:248432/100‑1 (MQ=255)
                                                                              gCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCATGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTg                     <  2:745025/100‑1 (MQ=255)
                                                                              gCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTg                     <  1:1596748/100‑1 (MQ=255)
                                                                              gCCAGTGCCCGGGTTAATTGTTTTTATGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAAAGGTCTTCATATTCATAtttttaaacagt                      >  2:264531/1‑92 (MQ=255)
                                                                                  gTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAAt                 <  1:3667252/100‑1 (MQ=255)
                                                                                   tGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATc                <  1:1986815/100‑1 (MQ=255)
                                                                                     cccGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCAc              >  2:3627887/1‑100 (MQ=255)
                                                                                         ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCg          <  2:3763269/100‑1 (MQ=255)
                                                                                         ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCg          >  1:2541726/1‑100 (MQ=255)
                                                                                         ggTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCg          <  1:3277569/100‑1 (MQ=255)
                                                                                          gTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGa         <  1:1724988/100‑1 (MQ=255)
                                                                                            tAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGc       <  1:1129716/100‑1 (MQ=255)
                                                                                            tAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGc       >  2:1472519/1‑100 (MQ=255)
                                                                                               ttGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCa     <  2:1482939/99‑1 (MQ=255)
                                                                                                tGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATa   <  1:1470700/100‑1 (MQ=255)
                                                                                                tGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATa   >  2:908366/1‑100 (MQ=255)
                                                                                                 gTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATAc  >  2:2824548/1‑100 (MQ=255)
                                                                                                 gTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATAc  <  1:2777534/100‑1 (MQ=255)
                                                                                                 gTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATAc  >  2:1082458/1‑100 (MQ=255)
                                                                                                  tttttCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATAc  >  2:2324114/1‑99 (MQ=255)
                                                                                                  |                                                                                                  
ATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATAC  >  NC_000913/1315607‑1315803

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: