Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,039,480 T→C 7.5% V62A (GTC→GCC) ‡ appB → cytochrome bd‑II ubiquinol oxidase subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,039,4800TC 7.5% 267.9 / 9.9 93V62A (GTC→GCC) ‡appBcytochrome bd‑II ubiquinol oxidase subunit II
Reads supporting (aligned to +/- strand):  ref base T (45/41);  new base C (5/2);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

CTGCTGATTGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGAAAAATCG  >  NC_000913/1039332‑1039626
                                                                                                                                                    |                                                                                                                                                  
ctgctgATTGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTc                                                                                                                                                   <  2:1311472/150‑1 (MQ=255)
 tgctgATTGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCt                                                                                                                                                  >  2:69969/1‑150 (MQ=255)
      aTTGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTg                                                                                                                                             >  1:901311/1‑150 (MQ=255)
      aTTGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTg                                                                                                                                             >  1:1435974/1‑150 (MQ=255)
        tGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTgcc                                                                                                                                           <  2:1342504/150‑3 (MQ=255)
          gCGTGATCCTGGTGGTCTTTATGATCTCCGACTGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCTTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATcc                                                                                                                                         <  2:135032/150‑1 (MQ=255)
           cGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCt                                                                                                                                        <  1:978256/150‑1 (MQ=255)
           cGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTgccttc                                                                                                                                        <  1:808120/150‑6 (MQ=255)
            gTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTc                                                                                                                                       >  2:578649/1‑150 (MQ=255)
              gagccTGGGGGTCTGTATGATCTCCGACGGATTGGACATGGGGATCGGCGGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGc                                                                                                                                     <  2:1505977/147‑1 (MQ=255)
              gATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGAGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGc                                                                                                                                     >  1:1494005/1‑150 (MQ=255)
                 ccTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCtgg                                                                                                                                  <  2:226919/150‑1 (MQ=255)
                    ggtggtCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCTGTTAGTGATAAAGAGCGTTGGTGCACACTGGGATGGCATCGAGGTCTGGTTGTTCCTCGCtggtgg                                                                                                                               <  2:1076286/150‑1 (MQ=255)
                       ggtCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGc                                                                                                                            <  2:367158/150‑1 (MQ=255)
                       gggcTTTATGATCTCCGACGGATTTGACATGGGGATCGTCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGTTTGATCCTCGCGGGTGGGGc                                                                                                                            <  2:1351999/147‑1 (MQ=255)
                              aTGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTAttt                                                                                                                     >  2:118856/1‑150 (MQ=255)
                               tGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCtggtgg                                                                                                                               >  2:421435/1‑139 (MQ=255)
                               tGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCtggtgg                                                                                                                               <  1:421435/139‑1 (MQ=255)
                               tGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTg                                                                                                                    >  1:685862/1‑150 (MQ=255)
                                  tctcCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTgccg                                                                                                                 <  2:920683/150‑1 (MQ=255)
                                      cGACTGATTTGACATGTGGATCGGCTGTCTTCTGCCGCTGGTGGCGCGTATTGATGATGAACGCCGGATAGTGATAAAGAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTg                                                                                                             <  1:705610/150‑1 (MQ=255)
                                      cGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTg                                                                                                             <  2:915429/150‑1 (MQ=255)
                                       gACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTgg                                                                                                            >  2:1290438/1‑150 (MQ=255)
                                         cGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGcc                                                                                                          >  1:812942/1‑150 (MQ=255)
                                           gATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGcccc                                                                                                        >  1:920628/1‑149 (MQ=255)
                                           gATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCa                                                                                                        <  2:1391947/150‑1 (MQ=255)
                                           gATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCa                                                                                                        >  1:140632/1‑150 (MQ=255)
                                            aTTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCag                                                                                                       >  1:983924/1‑150 (MQ=255)
                                             tttGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCaga                                                                                                      >  2:1150993/1‑150 (MQ=255)
                                                taCATGGGGATCGGCTGTCTGCTGCCGCTGTTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTTGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAgtg                                                                                                   <  2:1494005/149‑1 (MQ=255)
                                                gACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAgtg                                                                                                   >  2:742745/1‑150 (MQ=255)
                                                  caaggggATCGGCTCTCTGCTGCCGCTGGAGGAGCGTAATGATCATGAACGCCGGATAGAGATAAACAGCCTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTa                                                                                                 >  2:579274/4‑150 (MQ=255)
                                                    tGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAATCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATg                                                                                               >  1:548413/1‑150 (MQ=255)
                                                    tGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATg                                                                                               >  2:349836/1‑150 (MQ=255)
                                                    tGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATg                                                                                               >  2:936149/1‑150 (MQ=255)
                                                     ggggATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGc                                                                                              <  2:655012/150‑1 (MQ=255)
                                                      gggATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCa                                                                                             >  2:833318/1‑150 (MQ=255)
                                                       ggATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAg                                                                                            >  2:1468037/1‑150 (MQ=255)
                                                         aTCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCt                                                                                                                                    >  1:783481/1‑108 (MQ=255)
                                                         aTCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCt                                                                                                                                    <  2:783481/108‑1 (MQ=255)
                                                          tCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAgcgg                                                                                         <  1:349836/150‑1 (MQ=255)
                                                            ggCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAgcggcg                                                                                       >  1:210814/1‑150 (MQ=255)
                                                            ggCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAgcgg                                                                                         <  2:185408/148‑1 (MQ=255)
                                                            ggCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAgcgg                                                                                         >  1:185408/1‑148 (MQ=255)
                                                                  ctgctgCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTcc                                                                                 >  1:57878/1‑150 (MQ=255)
                                                                    gctgcCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCgg                                                                               >  2:355684/1‑150 (MQ=255)
                                                                       gcCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCtt                                                                            >  1:371403/1‑150 (MQ=255)
                                                                        ccGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCCTTGCTGCACACTGGGAAGGCAACCAGGTCTGGTTCATCCTCGCAGGTGGGCCATTATTTTACGCCTGGCCCAGAGTGTATGCAGCCGCGGTTTCCGGCtta                                                                           >  2:248875/1‑149 (MQ=255)
                                                                          gCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTa                                                                         >  1:200667/1‑150 (MQ=255)
                                                                          gCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTa                                                                         >  1:251780/1‑150 (MQ=255)
                                                                              gtggCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTAtgtg                                                                     >  1:369811/1‑150 (MQ=255)
                                                                                  cgcgTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCga                                                                 <  1:1161632/150‑1 (MQ=255)
                                                                                        atgatgatgAACGCCGGATAGTGATAAGCAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGTGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATcc                                                           <  2:461607/150‑1 (MQ=255)
                                                                                        atgatgatgAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATcc                                                           <  2:369811/150‑1 (MQ=255)
                                                                                            tgatgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGgt                                                       >  2:525626/1‑150 (MQ=255)
                                                                                              ttgaACGCCGGATTGTGATAAACAGGGTTGGTGCACACTGGGAAGGCAACGAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGgtgc                                                     <  2:532147/149‑1 (MQ=255)
                                                                                                 aaCGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATACTCGCTGGTGTGGCATTATTTGCCGACTTGCCCAGAGTGTATGCAGAGTCGTTTTCCTTCTTTTATTTGGCGATGATCCTGgtgctgt                                                  <  2:1435974/150‑1 (MQ=255)
                                                                                                  aCGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCACTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGgtgctgt                                                  <  1:322271/149‑1 (MQ=255)
                                                                                                  aCGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCACTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGgtgctgt                                                  >  2:322271/1‑149 (MQ=255)
                                                                                                   cGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGgtgctgtgc                                                <  2:812942/150‑1 (MQ=255)
                                                                                                     ccGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTc                                              >  2:59924/1‑150 (MQ=255)
                                                                                                     ccGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTc                                              <  1:90952/150‑1 (MQ=255)
                                                                                                        gATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACt                                           <  2:685862/150‑1 (MQ=255)
                                                                                                        gATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACt                                           <  1:616721/150‑1 (MQ=255)
                                                                                                        gATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACt                                           >  1:181234/1‑150 (MQ=255)
                                                                                                          tAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGGGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGt                                         <  1:297378/150‑1 (MQ=255)
                                                                                                            gTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGttc                                       >  1:1459212/1‑150 (MQ=255)
                                                                                                              gATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGCGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGttctt                                     >  2:1536654/1‑150 (MQ=255)
                                                                                                                tAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTcc                                   <  1:1322267/150‑1 (MQ=255)
                                                                                                                  aaCAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGCTCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTccgc                                 >  2:1159119/1‑150 (MQ=255)
                                                                                                                  aaCAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTccgc                                 <  1:1290438/150‑1 (MQ=255)
                                                                                                                           ggTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTccgcccg                              >  1:476859/1‑144 (MQ=255)
                                                                                                                           ggTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTccgcccg                              <  2:476859/144‑1 (MQ=255)
                                                                                                                           ggTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGcc                        <  1:525626/150‑1 (MQ=255)
                                                                                                                              gCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCttt                     >  1:1355691/1‑150 (MQ=255)
                                                                                                                              gCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCttt                     >  2:39673/1‑150 (MQ=255)
                                                                                                                                acacTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGa                   <  2:371403/150‑1 (MQ=255)
                                                                                                                                acacTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGa                   <  1:936149/150‑1 (MQ=255)
                                                                                                                                 cacTGGGAAGGCAACGAGGTCTGGGTGATCCTCGGTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGAt                  <  1:582152/150‑1 (MQ=255)
                                                                                                                                 cacTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGAt                  >  2:1067174/1‑150 (MQ=255)
                                                                                                                                       aggatcAACCAGACCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATcgc            >  1:295309/6‑150 (MQ=255)
                                                                                                                                       aggatcAACCAGACCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATcgc            >  2:751/6‑150 (MQ=255)
                                                                                                                                        ggatcAACCAGACCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATcgcg           >  1:1408383/5‑150 (MQ=255)
                                                                                                                                        ggatcAACCAGACCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTACTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATcgcg           >  2:1004217/5‑150 (MQ=255)
                                                                                                                                        aaGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATcgcg           <  1:833318/150‑1 (MQ=255)
                                                                                                                                        aaGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATcgcg           <  1:118856/150‑1 (MQ=255)
                                                                                                                                         gatcAACCAGACCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCgg          >  2:1254714/4‑150 (MQ=255)
                                                                                                                                          ggCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGa         >  1:10530/1‑150 (MQ=255)
                                                                                                                                           gCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGaa        <  2:251780/150‑1 (MQ=255)
                                                                                                                                            cAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGaaa       >  1:1252667/1‑150 (MQ=255)
                                                                                                                                            cAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGaaa       <  1:1150993/150‑1 (MQ=255)
                                                                                                                                              acaaGGTCTCGTTTATCCTACCAGGTGGGGAATTATTAGCAGACTGGCCCAGAGTTAATACAGCAAAATATTCAGGCTTTTATGTGGCCATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGaaaaa     >  2:1266586/4‑150 (MQ=255)
                                                                                                                                                 aGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGAAAAATCg  <  1:578649/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTGCTGATTGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGAAAAATCG  >  NC_000913/1039332‑1039626

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: