Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 448,367 A→C 5.8% V98V (GTT→GTG 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,3670AC 5.8% 368.4 / 7.6 125V98V (GTT→GTGcyoCcytochrome o ubiquinol oxidase subunit III
Reads supporting (aligned to +/- strand):  ref base A (58/58);  new base C (6/2);  total (65/60)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 126 total aligned reads displayed.

AGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACC  >  NC_000913/448218‑448511
                                                                                                                                                     |                                                                                                                                                
aGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATa                                                                                                                                                  <  2:654822/150‑1 (MQ=255)
   ccGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAAcc                                                                                                                                               >  2:661251/1‑150 (MQ=255)
    cgtgcgtgCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGAGGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCt                                                                                                                                              <  2:1525453/150‑1 (MQ=255)
       gcgtgcCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGc                                                                                                                                           <  1:1028604/150‑1 (MQ=255)
        cgtgcCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCt                                                                                                                                          >  2:1148880/1‑150 (MQ=255)
        cgtgcCGACCAGCGCAAAGAACGCTGACAGGAACCCGCTGCGATCCGGCCCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATACCGATAAATCCGGCACCAAACAACCAGGTCAACGACAGCCAGGAGATAACATAGCt                                                                                                                                          >  2:1294635/1‑150 (MQ=255)
         gtgcCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCtt                                                                                                                                         >  2:816033/1‑150 (MQ=255)
          tgcCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCttt                                                                                                                                        <  2:944191/150‑1 (MQ=255)
          ggcCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATATCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGTTATCTCCTGGCtgg                                                                                                                                        <  1:823303/149‑3 (MQ=255)
           gcCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAACCAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCtttt                                                                                                                                       >  1:126307/1‑150 (MQ=255)
           gcCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGGTAACAATCAGGTGATGGAATTCATAGATTTCCAGCCCGATAACTCCGGCACCAACCACCCAGGTCAACGCCAGCCAGGAGATCACCTGGCgttg                                                                                                                                       >  1:1416838/1‑149 (MQ=255)
            ccGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTttg                                                                                                                                      <  2:230125/150‑1 (MQ=255)
             cGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGTTATCTCCTGGCtggcgt                                                                                                                                     <  1:1318304/150‑6 (MQ=255)
             cGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTttgt                                                                                                                                     <  1:875154/150‑1 (MQ=255)
              gACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTttgtt                                                                                                                                    <  1:538290/150‑1 (MQ=255)
                 cAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTgttgtt                                                                                                                                 >  2:971162/1‑150 (MQ=255)
                  aGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGttt                                                                                                                                >  2:250692/1‑150 (MQ=255)
                     gcAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTg                                                                                                                             <  1:528591/150‑1 (MQ=255)
                     gcAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTg                                                                                                                             <  2:594827/150‑1 (MQ=255)
                      cAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGttttt                                                                                                                              >  2:1568661/1‑148 (MQ=255)
                       aaaGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTa                                                                                                                           <  1:794086/150‑1 (MQ=255)
                       aaaGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTa                                                                                                                           >  2:1418031/1‑150 (MQ=255)
                          gAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACAt                                                                                                                        >  2:96311/1‑150 (MQ=255)
                            aCGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATgg                                                                                                                      >  2:1413616/1‑150 (MQ=255)
                             cGCTGACAGGAAGCCTCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATATCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGGCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTTCATGGc                                                                                                                     <  1:592599/150‑1 (MQ=255)
                               cTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGa                                                                                                                   <  1:1725692/150‑1 (MQ=255)
                                 gACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATa                                                                                                                 <  1:1061535/150‑1 (MQ=255)
                                  aCAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAg                                                                                                                >  1:1375309/1‑150 (MQ=255)
                                   cAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAgc                                                                                                               >  2:1538679/1‑150 (MQ=255)
                                      gAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATCACCGGGCTTTTGTTGTTTTTGTACATGGCGATAgccgc                                                                                                            >  1:120739/1‑150 (MQ=255)
                                           cGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGc                                                                                                       <  1:657181/150‑1 (MQ=255)
                                              tGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGt                                                                                                    <  2:1506229/150‑1 (MQ=255)
                                               gCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTa                                                                                                   <  2:1449514/150‑1 (MQ=255)
                                                 gATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAgg                                                                                                 >  2:517643/1‑150 (MQ=255)
                                                     cGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGAt                                                                                             >  2:989882/1‑150 (MQ=255)
                                                      ggACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATg                                                                                            >  2:395477/1‑150 (MQ=255)
                                                         cccATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAg                                                                                         >  2:1633943/1‑150 (MQ=255)
                                                          ccATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACACCCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGc                                                                                        <  1:431480/150‑1 (MQ=255)
                                                          ccATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGc                                                                                        <  1:737541/150‑1 (MQ=255)
                                                            aTGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTg                                                                                      >  2:754622/1‑150 (MQ=255)
                                                              gCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCGTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGaa                                                                                    >  1:832754/1‑150 (MQ=255)
                                                                cGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGaaca                                                                                  <  1:244813/150‑1 (MQ=255)
                                                                 gTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGaacaa                                                                                 >  1:984463/1‑150 (MQ=255)
                                                                  ttAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGaacaac                                                                                <  1:816033/150‑1 (MQ=255)
                                                                   tAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAacaaca                                                                               <  2:1270757/150‑1 (MQ=255)
                                                                    aaCAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGTATATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTATGTGATGGAGCTGAACAACAg                                                                              <  1:1105628/150‑1 (MQ=255)
                                                                    aaCAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAg                                                                              <  1:923388/150‑1 (MQ=255)
                                                                            ggTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGGCGTAGGTGATGGAGCTGAACAACAGCAAGAAAg                                                                      <  1:1023699/150‑1 (MQ=255)
                                                                            ggTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAg                                                                      <  2:1697364/150‑1 (MQ=255)
                                                                            ggTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAg                                                                      <  2:1690820/150‑1 (MQ=255)
                                                                            ggTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAg                                                                      >  2:697242/1‑150 (MQ=255)
                                                                              tGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGtt                                                                    <  2:944231/150‑1 (MQ=255)
                                                                              ggATGGAAATCTATGAATTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATACCCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGtt                                                                    >  2:1166291/2‑150 (MQ=255)
                                                                               gATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATTAACTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCAAGACGTAGGTGATGGAGCTGAAAAAGAGCAAGAAAGttt                                                                   >  1:296327/1‑150 (MQ=255)
                                                                                aTGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTc                                                                  >  2:1712797/1‑150 (MQ=255)
                                                                                  ggAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCaa                                                                >  1:285274/1‑150 (MQ=255)
                                                                                   gAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAAc                                                               >  1:122197/1‑150 (MQ=255)
                                                                                    aaTTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAAcc                                                              >  2:744002/1‑150 (MQ=255)
                                                                                       tCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAga                                                           <  2:875017/150‑1 (MQ=255)
                                                                                       tCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAga                                                           >  1:1709327/1‑150 (MQ=255)
                                                                                        cATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAgaa                                                          >  1:426587/1‑150 (MQ=255)
                                                                                        cATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAgaa                                                          >  2:1037371/1‑150 (MQ=255)
                                                                                        cATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCGTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGATTCAACCAgaa                                                          >  1:1510682/1‑150 (MQ=255)
                                                                                           aGATTTCCATCCCGATAAATCCGGCCCCAAACAAACAGGTCAACGCCAGCCAGGAGATAACCTGGCTATTGATGTATTTGTACATGGCCATACCCGCCATGCCGTAGGTGATGGAGCAGAACAACAGCAAGAAAGTTTCAACCAgaacga                                                       >  1:922958/1‑150 (MQ=255)
                                                                                            gATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAgaacgaa                                                      >  1:888723/1‑150 (MQ=255)
                                                                                              tttCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATACCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGaacgaacg                                                    >  2:681700/1‑150 (MQ=255)
                                                                                                tccatccaGATAAATCCGGCACCAAACAACCAGGCCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGc                                                  <  2:101602/150‑1 (MQ=255)
                                                                                                tCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGc                                                  <  2:849228/150‑1 (MQ=255)
                                                                                                tCCATCCCGATAAATCCGGCACCAAACAACCAGCTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGc                                                  <  1:1678672/150‑1 (MQ=255)
                                                                                                     cccGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTc                                             <  1:1148880/150‑1 (MQ=255)
                                                                                                       cGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGa                                           >  2:180147/1‑150 (MQ=255)
                                                                                                          tAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGaaaa                                        <  1:1362060/150‑1 (MQ=255)
                                                                                                            aaTCCGGCACCAAACAACCAGGTCAACGCCAGCCTGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATATCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCATAACGAACGGCAGTTCGAAAATg                                      <  1:957211/150‑1 (MQ=255)
                                                                                                            aaTCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAAtt                                      >  1:1000141/1‑149 (MQ=255)
                                                                                                             aTCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGt                                     <  2:1510682/150‑1 (MQ=255)
                                                                                                              tCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTc                                    >  1:370643/1‑150 (MQ=255)
                                                                                                               ccGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTcc                                   <  2:1375309/150‑1 (MQ=255)
                                                                                                                 ggCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCtt                                 >  2:22918/1‑150 (MQ=255)
                                                                                                                   cACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTATGTAGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGCTGGAGCTGCACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTAc                               >  1:287737/1‑150 (MQ=255)
                                                                                                                    aCCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTAcc                              <  2:1117077/150‑1 (MQ=255)
                                                                                                                     ccAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCt                             >  2:1595557/1‑150 (MQ=255)
                                                                                                                     ccAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCt                             <  2:74110/150‑1 (MQ=255)
                                                                                                                      cAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTg                            >  2:392204/1‑150 (MQ=255)
                                                                                                                      cAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTg                            >  1:1440597/1‑150 (MQ=255)
                                                                                                                          caacCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCgg                        >  1:510339/1‑150 (MQ=255)
                                                                                                                                ggTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGgccgcc                  >  2:1345477/1‑150 (MQ=255)
                                                                                                                                 gTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCt                 >  1:1702264/1‑150 (MQ=255)
                                                                                                                                   caacaaaagccagGTTATCTCCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGc               >  2:1321000/8‑150 (MQ=255)
                                                                                                                                   caacaaaagccagGTTATCTCCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCGGc               >  1:78555/8‑150 (MQ=255)
                                                                                                                                   cAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGc               >  2:901265/1‑150 (MQ=255)
                                                                                                                                    aaCGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAATATGTCCTTACCTGTCGGGCCGCCTGCg              >  2:1611055/1‑150 (MQ=255)
                                                                                                                                     aCGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCgg             <  1:1418031/150‑1 (MQ=255)
                                                                                                                                      cGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGt            >  2:1218412/1‑150 (MQ=255)
                                                                                                                                       gcccgccaGGAGATAACCTTGCTTTTGTTGTTTTTGTACAGGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACATCAAGACAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTg           <  2:319854/150‑1 (MQ=255)
                                                                                                                                        ccagccagGAGATAACCTTGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGc          <  2:122197/150‑1 (MQ=255)
                                                                                                                                        ccagccagGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGc          <  1:901836/150‑1 (MQ=255)
                                                                                                                                          agccagGTTATCTCCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCg        >  1:799409/1‑150 (MQ=255)
                                                                                                                                             cagGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTc     >  2:795159/1‑150 (MQ=255)
                                                                                                                                                gagaTAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCAcc  <  2:789264/150‑1 (MQ=255)
                                                                                                                                                gagaTAACCTGGCTTTTGTTGTTTTAGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGACCGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCAcc  <  1:883450/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
AGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACC  >  NC_000913/448218‑448511

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