Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 448,368 A→T 5.3% V98D (GTT→GAT)  cyoC ← cytochrome o ubiquinol oxidase subunit III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913448,3680AT 5.3% 392.9 / 9.1 147V98D (GTT→GAT) cyoCcytochrome o ubiquinol oxidase subunit III
Reads supporting (aligned to +/- strand):  ref base A (65/72);  new base T (6/2);  total (72/75)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.84e-01

Only 100 of 151 total aligned reads displayed.

GACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAAC  >  NC_000913/448219‑448516
                                                                                                                                                     |                                                                                                                                                    
  ccGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAccc                                                                                                                                                    >  2:699319/1‑150 (MQ=255)
   cgtgcgtgCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAAccg                                                                                                                                                   >  2:1664648/1‑149 (MQ=255)
    gtgcgtgcCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTg                                                                                                                                                  >  2:1369890/1‑150 (MQ=255)
         tgcCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCttt                                                                                                                                             >  1:1721756/1‑150 (MQ=255)
            cGACCAGCGCAAAGACCGCTGACAGGAAGCCGCTGCGATCCGGACCCCTTCATTTAACAATCAGTTGATGGAATGCATAGAATTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTttgt                                                                                                                                          <  1:314557/150‑1 (MQ=255)
             gACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGTTATCTCCTGGCtggcgtt                                                                                                                                         <  2:1787979/150‑7 (MQ=255)
             gACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTttgtt                                                                                                                                         <  2:465585/150‑1 (MQ=255)
                     cAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGt                                                                                                                                 <  2:927148/150‑1 (MQ=255)
                     cAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGt                                                                                                                                 <  1:1625256/150‑1 (MQ=255)
                     cAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGt                                                                                                                                 >  1:161409/1‑150 (MQ=255)
                      aaaGAACGCTGACAGGAAGCCGCTGCGATACGGAACCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCAAGATAAATCCGGCACCAAACAACAAGGTCAACTCCAGCCAGGTGCTAACCTGGATTTTGTTGTTTTTGTa                                                                                                                                >  1:1085872/1‑150 (MQ=255)
                       aaGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTAc                                                                                                                               >  1:761760/1‑150 (MQ=255)
                         gAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACAt                                                                                                                             >  1:557454/1‑150 (MQ=255)
                            cGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGc                                                                                                                          <  1:2001821/150‑1 (MQ=255)
                                  cAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAgc                                                                                                                    >  2:498608/1‑150 (MQ=255)
                                   aGGACGCCGCTGCGATCCGGACCCATGCCGTTAATAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGGCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAgcc                                                                                                                   <  1:1948720/150‑1 (MQ=255)
                                    ggAAGCCGCGGCGATCCGGACCCATGCCGTTAACACTCAGGTGAGGGAATTCATAGATTTCCCTCCCGATAAATCCGGCACCAAACATCCAGGTCAAAGCCAGCCAGGAGAAAACCTGGCTTTTGTTGTTTTTGTACATGTCGATAgccg                                                                                                                  <  2:275754/150‑1 (MQ=255)
                                       aGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATACCCTGGCTGTTGTTGTTTTTGTACATGGCGATAGCCGCCa                                                                                                               >  1:1646511/1‑150 (MQ=255)
                                       aGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCa                                                                                                               >  1:142063/1‑150 (MQ=255)
                                        gCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGTTAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTATACATGGCGATAGCCGCCAt                                                                                                              <  1:1520342/150‑1 (MQ=255)
                                         ccGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAACTCCGCCACCAAACAACCAGGTCAACGCCCGCCAGGAGATAACCTGGCTTTTGATGTTTTTGTACATGGCGATAGCCGCCATg                                                                                                             >  1:347095/1‑150 (MQ=255)
                                         ccGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATg                                                                                                             <  1:256909/150‑1 (MQ=255)
                                          cGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGc                                                                                                            <  1:796085/150‑1 (MQ=255)
                                              gCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTa                                                                                                        >  2:32637/1‑150 (MQ=255)
                                              gCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTa                                                                                                        >  2:1832001/1‑150 (MQ=255)
                                                  tCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTg                                                                                                    <  2:1222049/150‑1 (MQ=255)
                                                  tCCGGACCCATGCCGTTAACAATCACGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTg                                                                                                    <  2:103227/150‑1 (MQ=255)
                                                    cGGACCCATGCCGTAAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGAt                                                                                                  >  2:1119545/1‑150 (MQ=255)
                                                        cccATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAg                                                                                              >  1:75040/1‑150 (MQ=255)
                                                        cccATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTATTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAg                                                                                              >  1:59891/1‑150 (MQ=255)
                                                           aTGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTg                                                                                           <  1:1950509/150‑1 (MQ=255)
                                                            tGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGa                                                                                          <  2:1984532/150‑1 (MQ=255)
                                                             gCCGTTAACACGCAGGCGATGGAATTCATCCATTTCCATCCCGATAAATCAGGCACCAAACATCCACGTCCACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATCGCGCTAGCCGCCATGCCGAAGGTGATGGAGCTGaa                                                                                         <  2:349533/150‑1 (MQ=255)
                                                              ccGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGTGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGaac                                                                                        <  2:24495/150‑1 (MQ=255)
                                                              ccGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGaac                                                                                        <  1:1442977/150‑1 (MQ=255)
                                                              ccGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGaac                                                                                        <  2:1059934/150‑1 (MQ=255)
                                                                 ttAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGCTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGaacaac                                                                                     <  1:531824/150‑1 (MQ=255)
                                                                     cAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCa                                                                                 >  2:245766/1‑150 (MQ=255)
                                                                      aaTCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCaa                                                                                >  1:123403/1‑150 (MQ=255)
                                                                           ggTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAg                                                                           <  2:1409728/150‑1 (MQ=255)
                                                                            gTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGt                                                                          >  1:40054/1‑150 (MQ=255)
                                                                              gATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGttt                                                                        <  2:1866557/150‑1 (MQ=255)
                                                                               aTGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTc                                                                       >  1:203030/1‑150 (MQ=255)
                                                                                tGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGATCAACAGCAAGAAAGTTTCa                                                                      >  2:728135/1‑150 (MQ=255)
                                                                                tGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCa                                                                      <  1:584172/150‑1 (MQ=255)
                                                                                tGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCa                                                                      >  2:1156586/1‑150 (MQ=255)
                                                                                tGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCa                                                                      <  1:1879081/150‑1 (MQ=255)
                                                                                tGGAATTCATAGATTTACATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGACAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATCGCGATAGCCGCAAAGCCGAAGGTGATGGAGATGAACAACAGCAAGAAAGTTGCa                                                                      >  2:408948/1‑150 (MQ=255)
                                                                                 ggAATTCATAGATTTCCATCCCGATAAATCCGGCACGACACACCCACGTCAACCCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTAGTCCATGGCGAAAGCCGCCATGCCGTAGGGGATGGAGCTGAACAACAGCAAGAAAGTTTCaa                                                                     <  1:91531/150‑1 (MQ=255)
                                                                                 ggAATTCATAGATTTCAATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCaa                                                                     <  2:557454/150‑1 (MQ=255)
                                                                                   aaTTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAAcc                                                                   <  2:1170552/150‑1 (MQ=255)
                                                                                   aaTTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAAcc                                                                   <  1:1601672/150‑1 (MQ=255)
                                                                                    aTTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCa                                                                  <  1:975621/150‑1 (MQ=255)
                                                                                     ttCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAg                                                                 <  2:1991759/150‑1 (MQ=255)
                                                                                     ttCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAAGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAg                                                                 >  1:1320796/1‑150 (MQ=255)
                                                                                        aTAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAgaac                                                              <  1:463076/150‑1 (MQ=255)
                                                                                         tAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAgaacg                                                             >  1:681153/1‑150 (MQ=255)
                                                                                             tttCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGaacgaacg                                                         >  2:286675/1‑150 (MQ=255)
                                                                                              ttCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACgg                                                        >  2:1772958/1‑150 (MQ=255)
                                                                                               tCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGc                                                       <  2:1070064/150‑1 (MQ=255)
                                                                                                 cATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCATCCAGAACGAACGGCAg                                                     >  1:1326851/1‑150 (MQ=255)
                                                                                                 cATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAg                                                     >  1:1543124/1‑150 (MQ=255)
                                                                                                     ccGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCg                                                 >  1:1729459/1‑150 (MQ=255)
                                                                                                        aTAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGaaa                                              <  2:1126361/150‑1 (MQ=255)
                                                                                                        aTAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGaaa                                              <  1:1891374/150‑1 (MQ=255)
                                                                                                         tAAATCCGTCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGaaaa                                             <  2:622791/150‑1 (MQ=255)
                                                                                                         tAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGaaaa                                             <  1:1450414/150‑1 (MQ=255)
                                                                                                           aaTCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATg                                           <  1:1259943/150‑1 (MQ=255)
                                                                                                            aTCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGt                                          >  1:899920/1‑150 (MQ=255)
                                                                                                            aTCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGGTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGGTCGAAAATgg                                          >  1:325518/1‑149 (MQ=255)
                                                                                                             tCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTc                                         <  1:869951/150‑1 (MQ=255)
                                                                                                                ggCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCtt                                      <  1:620934/150‑1 (MQ=255)
                                                                                                                     cAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTg                                 >  2:859397/1‑150 (MQ=255)
                                                                                                                          aacCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTCTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCggg                            <  2:52015/150‑1 (MQ=255)
                                                                                                                           acCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGATCGGCAGTTCGAAAATGTCCTTACCTGTCGGgc                           >  1:821452/1‑150 (MQ=255)
                                                                                                                             cAGGTCAACGACAGCCAGGAGATATCCTGGCTTTTGTTGTTTTTGTACATGGCGATTGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTGCTTACCTGTCGGgccg                         <  2:1019427/150‑1 (MQ=255)
                                                                                                                               ggTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGgccgcc                       >  2:1284877/1‑150 (MQ=255)
                                                                                                                                  cAACGCCAGCCACGAGATAACCTGGCTTTTGTTGGTGTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTGCGAAAATGTCCTTACCTGTCGGGCCGCCTGc                    <  1:1069655/150‑1 (MQ=255)
                                                                                                                                   aaCGTCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCg                   >  2:841074/1‑150 (MQ=255)
                                                                                                                                   aaCGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCAATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCg                   <  1:1016231/150‑1 (MQ=255)
                                                                                                                                    acaaaagccagGTTATCTCCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCgg                  >  1:1388549/6‑150 (MQ=255)
                                                                                                                                    aCGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCgg                  <  2:1657554/150‑1 (MQ=255)
                                                                                                                                    aCGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCgg                  <  2:40054/150‑1 (MQ=255)
                                                                                                                                     cGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGt                 <  2:1497064/150‑1 (MQ=255)
                                                                                                                                      aaaagccagGTTATCTCCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGCTGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTg                >  1:1869493/4‑150 (MQ=255)
                                                                                                                                      aaaagccagGTTATCTCCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTg                >  2:782116/4‑150 (MQ=255)
                                                                                                                                       aaagccagGTTATCTCCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGc               >  1:879403/3‑150 (MQ=255)
                                                                                                                                       aaagccagGTTATCTCCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGc               >  2:1343083/3‑150 (MQ=255)
                                                                                                                                        cagccagGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTTcc              <  1:1192601/150‑1 (MQ=255)
                                                                                                                                         agccagGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCg             <  2:1321398/150‑1 (MQ=255)
                                                                                                                                          gccagGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGt            >  2:872007/1‑150 (MQ=255)
                                                                                                                                              ggAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCAc        <  1:700158/150‑1 (MQ=255)
                                                                                                                                               gagaTAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCAcc       <  1:761826/150‑1 (MQ=255)
                                                                                                                                               gagaTAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCAcc       <  2:1504806/150‑1 (MQ=255)
                                                                                                                                               cagaTAACCTGGCTTTTGTTGTTTTTGTACATGGCCAAACCCGCCAAGCCGTAGGTGAGGGAGCTGAACAACAGCAAGAAAGTTACAACCAGAACGAACGGCAGTTCGATAATGACCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCAcc       <  1:1272016/149‑1 (MQ=255)
                                                                                                                                                agaTACCCTGGCTTTTGTAGTTTTTGTACTTGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACATCATCAAGAAAGTTTCAACCATAACGAACGGCAGATCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCa      <  1:1791691/150‑1 (MQ=255)
                                                                                                                                                 gaTAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAg     <  2:75040/150‑1 (MQ=255)
                                                                                                                                                  aTAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGa    >  2:1845546/1‑150 (MQ=255)
                                                                                                                                                    aaCCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAAc  >  2:1928831/1‑150 (MQ=255)
                                                                                                                                                    aaCCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGAGGGAGCTGAACAACAGCAAGAAAGTTGCAACCAGAACGAACTGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAAc  <  1:101324/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAAC  >  NC_000913/448219‑448516

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