Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 448,363 A→T 5.7% S100T (TCC→ACC)  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,3630AT 5.7% 417.5 / 7.4 143S100T (TCC→ACC) cyoCcytochrome o ubiquinol oxidase subunit III
Reads supporting (aligned to +/- strand):  ref base A (66/68);  new base T (7/2);  total (73/70)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

GTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCA  >  NC_000913/448214‑448512
                                                                                                                                                     |                                                                                                                                                     
gTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGga                                                                                                                                                       >  2:1693350/1‑150 (MQ=255)
 ttcAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATTGAATTCATTGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGgag                                                                                                                                                      <  2:1316236/148‑1 (MQ=255)
  gCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGCCCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGgaga                                                                                                                                                     >  1:983198/1‑150 (MQ=255)
  gCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGgaga                                                                                                                                                     >  2:899239/1‑150 (MQ=255)
  gCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGgaga                                                                                                                                                     <  1:1365915/150‑1 (MQ=255)
     gACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATaa                                                                                                                                                  >  1:907617/1‑150 (MQ=255)
       ccGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAccc                                                                                                                                                >  2:699319/1‑150 (MQ=255)
        cgtgcgtgCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAAccg                                                                                                                                               >  2:1664648/1‑149 (MQ=255)
        cgtgcgtgCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCt                                                                                                                                               <  2:1941092/150‑1 (MQ=255)
         gtgcgtgcCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTg                                                                                                                                              >  2:1369890/1‑150 (MQ=255)
              tgcCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCttt                                                                                                                                         >  1:1721756/1‑150 (MQ=255)
                  gACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGTTATCTCCTGGCtggcgtt                                                                                                                                     <  2:1787979/150‑7 (MQ=255)
                  gACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTttgtt                                                                                                                                     <  1:1243996/150‑1 (MQ=255)
                  gACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTttgtt                                                                                                                                     <  2:465585/150‑1 (MQ=255)
                         gcAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGGTAACAATCCGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTg                                                                                                                              >  1:111286/1‑150 (MQ=255)
                          cAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGt                                                                                                                             <  2:927148/150‑1 (MQ=255)
                          cAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGt                                                                                                                             <  1:1625256/150‑1 (MQ=255)
                           aaaGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTa                                                                                                                            >  2:761826/1‑150 (MQ=255)
                           aaaGAACGCTGACAGGAAGCCGCTGCGATACGGAACCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCAAGATAAATCCGGCACCAAACAACAAGGTCAACTCCAGCCAGGTGCTAACCTGGATTTTGTTGTTTTTGTa                                                                                                                            >  1:1085872/1‑150 (MQ=255)
                             aGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACa                                                                                                                          <  1:1623120/150‑1 (MQ=255)
                              gAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACAt                                                                                                                         >  1:557454/1‑150 (MQ=255)
                                 cGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGc                                                                                                                      <  1:2001821/150‑1 (MQ=255)
                                 cGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGCTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGc                                                                                                                      >  2:142282/1‑150 (MQ=255)
                                 cGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATACCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGc                                                                                                                      >  1:1896266/1‑150 (MQ=255)
                                    tGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGTTAACCTGGCTTTTGTTGTTTTTGTACATGGCGaa                                                                                                                   >  1:1957742/1‑149 (MQ=255)
                                    tGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGAt                                                                                                                   <  1:572482/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)
                                            aGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCa                                                                                                           >  1:142063/1‑150 (MQ=255)
                                            aGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCa                                                                                                           >  2:208267/1‑150 (MQ=255)
                                             gCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGTTAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTATACATGGCGATAGCCGCCAt                                                                                                          <  1:1520342/150‑1 (MQ=255)
                                             gCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCAt                                                                                                          >  2:1360982/1‑150 (MQ=255)
                                              ccGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATg                                                                                                         <  1:256909/150‑1 (MQ=255)
                                              ccGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATg                                                                                                         >  2:853071/1‑150 (MQ=255)
                                                gCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGcc                                                                                                       >  1:1740273/1‑150 (MQ=255)
                                                   gCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTa                                                                                                    >  2:1832001/1‑150 (MQ=255)
                                                   gCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTa                                                                                                    >  2:32637/1‑150 (MQ=255)
                                                     gATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAACCAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTTTTTTTGTACATGGCGATAGCCGCCATGCCGTAgg                                                                                                  >  2:28984/1‑150 (MQ=255)
                                                       tCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTg                                                                                                <  2:1222049/150‑1 (MQ=255)
                                                       tCCGGACCCATGCCGTTAACAATCACGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTg                                                                                                <  2:103227/150‑1 (MQ=255)
                                                        ccGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGa                                                                                               <  1:608883/150‑1 (MQ=255)
                                                         cGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGAt                                                                                              >  1:359485/1‑150 (MQ=255)
                                                         cGGACCCATGCCGTAAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGAt                                                                                              >  2:1119545/1‑150 (MQ=255)
                                                             cccATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAg                                                                                          >  1:75040/1‑150 (MQ=255)
                                                             cccATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTATTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAg                                                                                          >  1:59891/1‑150 (MQ=255)
                                                               cATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCt                                                                                        >  2:1856747/1‑150 (MQ=255)
                                                               cATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCt                                                                                        <  2:1547068/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)
                                                                           aaTCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCaa                                                                            >  1:123403/1‑150 (MQ=255)
                                                                                cgggATGGAAATCTATGAATTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAg                                                                       >  1:472006/4‑150 (MQ=255)
                                                                                 gTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGt                                                                      >  1:1318342/1‑150 (MQ=255)
                                                                                    aTGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTc                                                                   >  1:1657554/1‑150 (MQ=255)
                                                                                     tGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATTACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGTACAACAGCAAGTAAGTTTCa                                                                  >  1:1030461/1‑150 (MQ=255)
                                                                                     tGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGATCAACAGCAAGAAAGTTTCa                                                                  >  2:728135/1‑150 (MQ=255)
                                                                                     tGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCa                                                                  >  2:1156586/1‑150 (MQ=255)
                                                                                     tGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCa                                                                  <  1:1879081/150‑1 (MQ=255)
                                                                                     tGGAATTCATAGATTTACATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGACAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATCGCGATAGCCGCAAAGCCGAAGGTGATGGAGATGAACAACAGCAAGAAAGTTGCa                                                                  >  2:408948/1‑150 (MQ=255)
                                                                                      ggAATTCATAGATTTCAATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCaa                                                                 <  2:557454/150‑1 (MQ=255)
                                                                                        aaTTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAAcc                                                               <  1:1601672/150‑1 (MQ=255)
                                                                                        aaTTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAAcc                                                               <  2:1170552/150‑1 (MQ=255)
                                                                                         aTTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCa                                                              <  1:975621/150‑1 (MQ=255)
                                                                                            cATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAgaa                                                           >  1:762432/1‑150 (MQ=255)
                                                                                            cATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGATTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAgaa                                                           >  1:317619/1‑150 (MQ=255)
                                                                                             aTAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAgaac                                                          <  1:463076/150‑1 (MQ=255)
                                                                                              tAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAgaacg                                                         >  1:681153/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)
                                                                                                            gATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGaa                                           >  1:859318/1‑150 (MQ=255)
                                                                                                             aTAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGaaa                                          <  2:1126361/150‑1 (MQ=255)
                                                                                                              tAAATCCGTCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGaaaa                                         <  2:622791/150‑1 (MQ=255)
                                                                                                              tAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGaaaa                                         <  1:1450414/150‑1 (MQ=255)
                                                                                                                 aTCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGGTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGGTCGAAAATgg                                      >  1:325518/1‑149 (MQ=255)
                                                                                                                          cAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTg                             <  2:248719/150‑1 (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)
                                                                                                                                     ggcAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCt                  <  1:208267/148‑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:40054/150‑1 (MQ=255)
                                                                                                                                           gccagccaGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTg            <  1:757545/150‑1 (MQ=255)
                                                                                                                                           aaaagccagGTTATCTCCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGAAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTg            >  1:1877014/4‑150 (MQ=255)
                                                                                                                                            aaagccagGTTATCTCCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGc           >  1:879403/3‑150 (MQ=255)
                                                                                                                                              agccagGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCg         <  1:899239/150‑1 (MQ=255)
                                                                                                                                              agccagGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCg         <  2:1321398/150‑1 (MQ=255)
                                                                                                                                               gccagGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGt        >  2:872007/1‑150 (MQ=255)
                                                                                                                                                   ggAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCAc    <  1:699319/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)
                                                                                                                                                     |                                                                                                                                                     
GTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCA  >  NC_000913/448214‑448512

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: