Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 758,306 T→G 15.8% C206W (TGT→TGG sdhB → succinate dehydrogenase, FeS subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913758,3060TG15.8% 227.1 / 35.5 94C206W (TGT→TGGsdhBsuccinate dehydrogenase, FeS subunit
Reads supporting (aligned to +/- strand):  ref base T (43/36);  new base G (7/8);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 7.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

CCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCC  >  NC_000913/758158‑758454
                                                                                                                                                    |                                                                                                                                                    
ccTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGt                                                                                                                                                      >  2:1457972/1‑149 (MQ=255)
  tctTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTca                                                                                                                                                    >  1:1132715/1‑149 (MQ=255)
  tctTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTca                                                                                                                                                    <  2:416392/149‑1 (MQ=255)
  tctTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTca                                                                                                                                                    <  2:573455/149‑1 (MQ=255)
   ctTGTCCGTCTTTCTGGTGGAATCCCGATAAGGTTATCGGCCCGGCAGGCTTGTTATCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTTATGCATTCAGCGTATTCCGCTGTcac                                                                                                                                                   <  1:163045/149‑1 (MQ=255)
    ttGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTcaca                                                                                                                                                  <  2:124780/149‑1 (MQ=255)
       tCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGca                                                                                                                                               >  2:1502242/1‑149 (MQ=255)
        ccGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGcat                                                                                                                                              <  2:1048869/149‑1 (MQ=255)
        ccGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGcat                                                                                                                                              <  2:1054204/149‑1 (MQ=255)
         cGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCTTATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGcatc                                                                                                                                             <  2:678287/149‑1 (MQ=255)
         agTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGcatc                                                                                                                                             <  2:840893/148‑1 (MQ=255)
          gTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGcatca                                                                                                                                            <  2:545648/149‑1 (MQ=255)
          gTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGcatca                                                                                                                                            >  1:916821/1‑149 (MQ=255)
           tCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGcatcat                                                                                                                                           >  2:524912/1‑149 (MQ=255)
                  ggtggAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCAAATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGc                                                                                                                                    <  2:856793/149‑1 (MQ=255)
                   gtggAATCCCGATAATTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCg                                                                                                                                   >  2:934948/1‑149 (MQ=255)
                   gtggAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTACTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCg                                                                                                                                   >  1:1344489/1‑149 (MQ=255)
                      gAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCa                                                                                                                                >  1:16461/1‑149 (MQ=255)
                        aacccGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAgt                                                                                                                              >  2:1085672/3‑149 (MQ=255)
                        aTCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAgt                                                                                                                              >  2:1250688/1‑149 (MQ=255)
                        aTCCCGATAAGTTTATCGGCCCGGCAGGCGTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAgt                                                                                                                              <  1:1221870/149‑1 (MQ=255)
                           ccGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTa                                                                                                                           >  1:515212/1‑149 (MQ=255)
                            cGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTAt                                                                                                                          <  1:1397930/149‑1 (MQ=255)
                             gATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATg                                                                                                                         >  1:709350/1‑149 (MQ=255)
                               tcaggttATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTGCCTGTTTGATAGCCGTGATACCGAGACTGACCGCCGCCTCGACGGTTTGAGGGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTc                                                                                                                       <  1:304924/144‑1 (MQ=255)
                                  gTTTATCGGCCCGACAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGGCACAGCATCATGAACTGCGTCAGTGTATGTCCGa                                                                                                                    >  2:1458882/1‑149 (MQ=255)
                                   tttATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGaa                                                                                                                   <  1:164850/149‑1 (MQ=255)
                                   tttATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGaa                                                                                                                   >  1:1316763/1‑149 (MQ=255)
                                   tttATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCAGaa                                                                                                                   >  1:116638/1‑149 (MQ=255)
                                   tttATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGGCACAGCATCATGAACTGCGTCAGTGTATGTCCGaa                                                                                                                   >  1:45844/1‑149 (MQ=255)
                                   tttATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGGCACAGCATCATGAACTGCGTCAGTGTATGTCCGaa                                                                                                                   >  1:1218679/1‑149 (MQ=255)
                                       tCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCGCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGGCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAgggg                                                                                                               >  2:191600/1‑149 (MQ=255)
                                         ggCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCt                                                                                                             >  1:19759/1‑149 (MQ=255)
                                             cGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGGCACAGCATCATGAACTGCGTCAGTGTACGTCCGAAGGGGCTGAAc                                                                                                         <  1:130521/149‑1 (MQ=255)
                                                  ggCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGGCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGAc                                                                                                    <  1:1458882/149‑1 (MQ=255)
                                                   gCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGAcg                                                                                                   >  2:1138645/1‑149 (MQ=255)
                                                   gCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGGCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGACCCCGAcg                                                                                                   >  2:1069832/1‑149 (MQ=255)
                                                              gCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCgg                                                                                        <  2:656149/149‑1 (MQ=255)
                                                               cATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGc                                                                                       >  2:1410311/1‑149 (MQ=255)
                                                                     gTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGGCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATc                                                                                 <  2:1218679/149‑1 (MQ=255)
                                                                      tttCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCa                                                                                >  1:1415317/1‑149 (MQ=255)
                                                                      tttCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCa                                                                                >  2:1283730/1‑149 (MQ=255)
                                                                          cTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCGGCGTATTCCGCTGGCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTc                                                                            >  2:206790/1‑149 (MQ=255)
                                                                           tgattgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGATCCCGACGCGCGCCATCGGCCATATCAAGTCg                                                                           >  1:136642/1‑149 (MQ=255)
                                                                           tgattgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCg                                                                           >  2:963762/1‑149 (MQ=255)
                                                                           tgattgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATACCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCCTATCAAGTCg                                                                           >  1:131620/1‑149 (MQ=255)
                                                                            gattgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCAGAGCATCATGAACTGCGTCAGTGTATTTCCGAATGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGa                                                                          <  2:1020569/149‑1 (MQ=255)
                                                                             attgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGAATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCATGTCGaa                                                                         >  1:1504219/1‑148 (MQ=255)
                                                                              ttgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGAtg                                                                        <  1:524912/149‑1 (MQ=255)
                                                                                 atAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGAtgttg                                                                     <  1:1502242/149‑1 (MQ=255)
                                                                                  tAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGAtgttgt                                                                    >  1:940588/1‑149 (MQ=255)
                                                                                  tAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGGCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGt                                                                             >  2:699672/1‑140 (MQ=255)
                                                                                  tAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGGCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGt                                                                             <  1:699672/140‑1 (MQ=255)
                                                                                   aGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATgttgtt                                                                   <  2:1132715/149‑1 (MQ=255)
                                                                                   aGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATgttgtt                                                                   <  2:1344489/149‑1 (MQ=255)
                                                                                       gTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATgttgtttaaa                                                               >  2:319993/1‑145 (MQ=255)
                                                                                       gTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCaa                                                               <  2:1116681/149‑1 (MQ=255)
                                                                                       gTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGGCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCaa                                                               <  2:1266183/149‑1 (MQ=255)
                                                                                        tGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAAc                                                              >  2:1497252/1‑149 (MQ=255)
                                                                                        tGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAAc                                                              <  1:1520471/149‑1 (MQ=255)
                                                                                        tGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAAc                                                              <  1:1526515/149‑1 (MQ=255)
                                                                                          aTACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATgttgtt                                                                   >  1:954345/1‑142 (MQ=255)
                                                                                          aTACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATgttgtt                                                                   <  2:954345/142‑1 (MQ=255)
                                                                                            aCCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTaa                                                          >  2:123194/1‑149 (MQ=255)
                                                                                            aCCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCCTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTaa                                                          >  2:1176967/1‑149 (MQ=255)
                                                                                             ccGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAAt                                                         >  1:659199/1‑149 (MQ=255)
                                                                                              cGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGATGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATg                                                        >  1:250665/1‑149 (MQ=255)
                                                                                               gagaCTGCCAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATAATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAATTCGATGTTGTTGCAACGTAATGc                                                       >  1:274765/1‑149 (MQ=255)
                                                                                                    tGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTaaa                                                  <  1:1250688/149‑1 (MQ=255)
                                                                                                    tGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTaaa                                                  >  2:1445352/1‑149 (MQ=255)
                                                                                                    tGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTaaa                                                  >  1:927959/1‑149 (MQ=255)
                                                                                                    tGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTACGCTGTCACAGCATCATGAACTGCTTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTaaa                                                  >  2:517134/1‑149 (MQ=255)
                                                                                                        aGCCGCCTCGTCGGTTTGTGTGATGCATTCAGCGTATTCCGCTGGCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTTAACCCGACGCGCTCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGt                                              <  1:603204/149‑1 (MQ=255)
                                                                                                          ccgccTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAg                                            >  2:523765/1‑149 (MQ=255)
                                                                                                          ccgccTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAg                                            <  1:883496/149‑1 (MQ=255)
                                                                                                              cTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCt                                        <  2:808457/149‑1 (MQ=255)
                                                                                                              cTCGACGGTTTGAGTGATGCATGCAGCGTATTCCGCTGGCACAGCAGCAGGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCt                                        <  2:784627/149‑1 (MQ=255)
                                                                                                                cGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGa                                      <  2:116638/149‑1 (MQ=255)
                                                                                                                cGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGa                                      <  2:1419362/149‑1 (MQ=255)
                                                                                                                 gACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGAt                                     <  1:764115/149‑1 (MQ=255)
                                                                                                                   cGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATaa                                   <  1:523765/149‑1 (MQ=255)
                                                                                                                          aGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCGACGTAATGCGTAAACCGTAGGCCTGATAAGACgcgc                            >  2:334558/1‑149 (MQ=255)
                                                                                                                                  aTTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCg                    >  1:1521142/1‑149 (MQ=255)
                                                                                                                                   ttCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGc                   >  1:1260874/1‑149 (MQ=255)
                                                                                                                                      aGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGCCGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATc                >  1:751258/1‑149 (MQ=255)
                                                                                                                                      aGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATc                >  2:466012/1‑149 (MQ=255)
                                                                                                                                      aGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATc                >  1:816861/1‑149 (MQ=255)
                                                                                                                                       gCGGATTCCGCTGTCACAGCAGCATGAACTGCGTCAGTGTATGTCCGAAGGGGCGGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTGGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCa               <  2:108142/149‑1 (MQ=255)
                                                                                                                                         gTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAgg             <  2:1521142/149‑1 (MQ=255)
                                                                                                                                          tATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGc            >  1:289436/1‑149 (MQ=255)
                                                                                                                                           aTTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATGCGACGCTTGCGCGTCTTATCAGGcc           <  2:1316763/149‑2 (MQ=255)
                                                                                                                                             tCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAAc         <  2:659199/149‑1 (MQ=255)
                                                                                                                                             tCCGCTGGCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAAc         <  1:191600/149‑1 (MQ=255)
                                                                                                                                               cGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg              <  2:61431/142‑1 (MQ=255)
                                                                                                                                               cGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATAACAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg              >  1:61431/1‑142 (MQ=255)
                                                                                                                                                    tCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGcc  >  1:804102/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCTCTTGTCCGTCTTTCTGGTGGAATCCCGATAAGTTTATCGGCCCGGCAGGCTTGTTAGCGGCATATCGTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCC  >  NC_000913/758158‑758454

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