Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 758,324 C→A 8.3% C212* (TGC→TGA sdhB → succinate dehydrogenase, FeS subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913758,3240CA 8.3% 1598.9 / 93.5 454C212* (TGC→TGAsdhBsuccinate dehydrogenase, FeS subunit
Reads supporting (aligned to +/- strand):  ref base C (227/188);  new base A (26/13);  total (253/201)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 462 total aligned reads displayed.

TGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCT  >  NC_000913/758175‑758473
                                                                                                                                                     |                                                                                                                                                     
  gtggAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGt                                                                                                                                                     <  2:2896498/150‑1 (MQ=255)
  gtggAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGt                                                                                                                                                     <  1:2153968/150‑1 (MQ=255)
    ggAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCa                                                                                                                                                   >  2:784825/1‑150 (MQ=255)
      aaTCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAgt                                                                                                                                                 <  1:1385149/150‑1 (MQ=255)
      aaTCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAgt                                                                                                                                                 <  2:495480/150‑1 (MQ=255)
       aTCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGAGTCAgtg                                                                                                                                                <  2:2722604/150‑1 (MQ=255)
        tCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAgtgt                                                                                                                                               <  1:5840789/150‑1 (MQ=255)
          ccGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTAt                                                                                                                                             <  2:3921998/150‑1 (MQ=255)
           cGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATg                                                                                                                                            >  2:2459526/1‑150 (MQ=255)
                 gTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGaa                                                                                                                                      >  1:630638/1‑150 (MQ=255)
                    tATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAggg                                                                                                                                   <  1:2679326/150‑1 (MQ=255)
                     aTCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAgggg                                                                                                                                  <  1:1273820/150‑1 (MQ=255)
                       cGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCt                                                                                                                                <  1:853718/150‑1 (MQ=255)
                         gCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGAGTCAGTGTATGTCCGAAGGGGCTGa                                                                                                                              >  1:3657090/1‑150 (MQ=255)
                          cccGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGaa                                                                                                                             <  1:5299230/150‑1 (MQ=255)
                            cGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAAcc                                                                                                                           <  2:4797440/150‑1 (MQ=255)
                                aGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGAc                                                                                                                       >  2:2827382/1‑150 (MQ=255)
                                aGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGAc                                                                                                                       >  1:4218509/1‑150 (MQ=255)
                                    ttGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGAcgcgc                                                                                                                   <  2:5313438/150‑1 (MQ=255)
                                    ttGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGAcgcgc                                                                                                                   >  2:5890315/1‑150 (MQ=255)
                                      gTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACgcgcgc                                                                                                                 >  1:1531943/1‑150 (MQ=255)
                                          gCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATc                                                                                                             >  1:1605797/1‑150 (MQ=255)
                                           cGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCg                                                                                                            <  2:4454391/150‑1 (MQ=255)
                                           cGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCg                                                                                                            >  2:2068597/1‑150 (MQ=255)
                                            ggCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCgg                                                                                                           >  1:3580298/1‑150 (MQ=255)
                                             gCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGc                                                                                                          >  2:4436197/1‑150 (MQ=255)
                                              cATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGcc                                                                                                         <  2:1194880/150‑1 (MQ=255)
                                               atatCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCa                                                                                                        >  1:3347425/1‑150 (MQ=255)
                                                 atCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCata                                                                                                      >  2:5436706/1‑150 (MQ=255)
                                                        ccTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATc                                                                                                             >  2:1859508/1‑136 (MQ=255)
                                                          tgattgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGa                                                                                             >  2:4008672/1‑150 (MQ=255)
                                                          tgattgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGa                                                                                             >  1:3956463/1‑150 (MQ=255)
                                                          tgattgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGa                                                                                             >  1:2826961/1‑150 (MQ=255)
                                                               gatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGAtgttg                                                                                        <  1:4678363/150‑1 (MQ=255)
                                                                 tAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATgttgtt                                                                                      <  1:1451197/150‑1 (MQ=255)
                                                                     cGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCaa                                                                                  >  2:1688679/1‑150 (MQ=255)
                                                                       tGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACg                                                                                >  2:3568796/1‑150 (MQ=255)
                                                                         aTACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGAGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTa                                                                              >  2:5597025/1‑150 (MQ=255)
                                                                         aTACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGAGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGAGGTTGTTGCAACGTa                                                                              >  2:2623244/1‑150 (MQ=255)
                                                                            ccGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGatcaagtcgatgttgttgcaacgtaatgcgtaa                                                                           <  1:4591116‑M1/150‑34 (MQ=255)
                                                                            ccGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATg                                                                           <  1:5659936/150‑1 (MQ=255)
                                                                             cGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGc                                                                          >  2:3310661/1‑150 (MQ=255)
                                                                             cGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGc                                                                          >  2:713359/1‑150 (MQ=255)
                                                                             cGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGc                                                                          >  1:668090/1‑150 (MQ=255)
                                                                              gagaCTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCg                                                                         >  2:5727988/1‑150 (MQ=255)
                                                                                   tGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAAc                                                                    >  2:4288682/1‑150 (MQ=255)
                                                                                   tGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAAc                                                                    >  1:4507288/1‑150 (MQ=255)
                                                                                   tGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAAc                                                                    >  2:1256398/1‑150 (MQ=255)
                                                                                     aCAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCg                                                                  >  1:969654/1‑150 (MQ=255)
                                                                                     aCAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATCAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCg                                                                  <  2:4101090/150‑1 (MQ=255)
                                                                                      cAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGAGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGt                                                                 >  1:3820923/1‑150 (MQ=255)
                                                                                       aGCCGCCTCGACGGTGTGAGTGATGCATTCAGGGTTTTCCGCTGTCACAGCATGATGAACTGCGTCAGTGTATGTCCGTAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTa                                                                <  1:3180200/150‑1 (MQ=255)
                                                                                           gccTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGAGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGcc                                                            >  2:1545187/1‑150 (MQ=255)
                                                                                              tCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGa                                                         <  2:3755898/150‑1 (MQ=255)
                                                                                              tCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGa                                                         >  1:1814413/1‑150 (MQ=255)
                                                                                                   ggTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGa                                                    <  1:732523/150‑1 (MQ=255)
                                                                                                      ttGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGAcgc                                                 <  2:2291642/150‑1 (MQ=255)
                                                                                                          gTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCaa                                             <  2:342967/150‑1 (MQ=255)
                                                                                                          gTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCaa                                             <  1:2026090/150‑1 (MQ=255)
                                                                                                           TGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCgcacg                                            <  2:3232096‑M1/150‑3 (MQ=255)
                                                                                                           TGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAg                                            >  1:1100957‑M1/1‑150 (MQ=255)
                                                                                                           TGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAg                                            <  1:2074853‑M1/150‑1 (MQ=255)
                                                                                                           TGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAg                                            <  1:3501317‑M1/150‑1 (MQ=255)
                                                                                                            GATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGc                                           >  2:1678568‑M1/1‑150 (MQ=255)
                                                                                                            GATGCATTCAGCGTATTCCGCTATCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGc                                           >  1:3919459‑M1/1‑150 (MQ=255)
                                                                                                             ATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCg                                          <  1:4295505‑M1/150‑1 (MQ=255)
                                                                                                              TGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGt                                         >  1:3045871‑M1/1‑150 (MQ=255)
                                                                                                              TGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGt                                         <  1:4081865‑M1/150‑1 (MQ=255)
                                                                                                               gctgtgacagcggaatacgctgaatgcatcaTGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTc                                        >  1:1032354‑M2/32‑150 (MQ=255)
                                                                                                               GCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTc                                        <  1:2433909‑M1/150‑1 (MQ=255)
                                                                                                                CATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCg                                       <  1:5222787‑M1/150‑1 (MQ=255)
                                                                                                                 aTTCAGCGTATTCCGCTGTCACAGCATCATGATCTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGc                                      <  1:5355320/150‑1 (MQ=255)
                                                                                                                 aTTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGc                                      >  1:4787662/1‑150 (MQ=255)
                                                                                                                  ttCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCa                                     >  1:437820/1‑150 (MQ=255)
                                                                                                                  ttCAGCGTATTCCGCTGTCACAGCATCATGAACTGAGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCa                                     <  1:2351299/150‑1 (MQ=255)
                                                                                                                   tCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCAt                                    <  1:1179898/150‑1 (MQ=255)
                                                                                                                    cAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATc                                   >  2:60990/1‑150 (MQ=255)
                                                                                                                    cAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATc                                   >  1:5725705/1‑150 (MQ=255)
                                                                                                                       cGTATTCCGCTGTCACAGCATCATGAACTGAGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCCTATCAAGTCGCTGGTGTTGCAACGTACTGCGTAAACCGGAGGCCGGCTAAGACGCGCACGCGGCGCATCAgg                                >  2:3880329/1‑150 (MQ=255)
                                                                                                                           ttCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAAc                            <  1:3397930/150‑1 (MQ=255)
                                                                                                                             ccGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCa                          >  1:4948611/1‑150 (MQ=255)
                                                                                                                              cGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAg                         >  1:4490846/1‑150 (MQ=255)
                                                                                                                               gCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGt                        >  2:4623583/1‑150 (MQ=255)
                                                                                                                                   tCACAGCATCATGAACTGAGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCg                    >  1:1724588/1‑150 (MQ=255)
                                                                                                                                     acaGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGa                  >  1:4119402/1‑150 (MQ=255)
                                                                                                                                       aGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATg                <  2:3061486/150‑1 (MQ=255)
                                                                                                                                        gcatcaTGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATgc               >  1:3277034‑M2/7‑150 (MQ=255)
                                                                                                                                        gcatcaTGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATgc               >  2:5972761‑M2/7‑150 (MQ=255)
                                                                                                                                           tcaTGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCCACCAGTGCCGGATgcggc            >  2:1069527‑M2/4‑150 (MQ=255)
                                                                                                                                           tcaTGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTGGTTGCAACGTGATGCGTAAACCGGAGGCCTGATAAGACGCGCAAGCGTCGCCTCCGGCAACCAGTGCCGGATgcggc            >  2:847811‑M2/4‑150 (MQ=255)
                                                                                                                                            caTGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATgcggcg           >  2:2397348‑M2/3‑150 (MQ=255)
                                                                                                                                            caTGAACTGCGTCAGTGTATGACCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATgcggcg           <  1:600688‑M2/148‑1 (MQ=255)
                                                                                                                                             aTGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGt          >  1:1662004‑M2/2‑150 (MQ=255)
                                                                                                                                             aTGAACTGAGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGt          >  2:4990462‑M2/2‑150 (MQ=255)
                                                                                                                                              tGAACTGCGTCAGTGTATGTCCGACGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTg         <  1:1898475/150‑1 (MQ=255)
                                                                                                                                              tGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTg         >  2:4742023/1‑150 (MQ=255)
                                                                                                                                                   tGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGc    <  2:5663758/150‑1 (MQ=255)
                                                                                                                                                   tGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGc    <  2:4155525/150‑1 (MQ=255)
                                                                                                                                                    gCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGcc   <  1:294388/150‑1 (MQ=255)
                                                                                                                                                     cGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCt  >  2:5181743/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCT  >  NC_000913/758175‑758473

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: