Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,841,159 T→A 5.5% L13Q (CTG→CAG)  ccoO‑II → cbb3‑type cytochrome c oxidase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,841,1590TA 5.5% 284.2 / 5.6 91L13Q (CTG→CAG) ccoO‑IIcbb3‑type cytochrome c oxidase subunit
Reads supporting (aligned to +/- strand):  ref base T (42/44);  new base A (2/3);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

TTCGCCAGCGGCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAGGCCGCGATATCTACATCCGCG  >  NC_002947/4841012‑4841305
                                                                                                                                                   |                                                                                                                                                  
ttcGCCAGCGGCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTgg                                                                                                                                                  >  1:771620/1‑150 (MQ=255)
  cGCCAGCGGCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGcc                                                                                                                                                <  1:488386/150‑1 (MQ=255)
    ccAGCGGCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCctt                                                                                                                                              >  1:946430/1‑150 (MQ=255)
    ccAGCGGCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAACACGTCGAGAAGAACATAGGCCTGCTGGCctt                                                                                                                                              >  2:979505/1‑150 (MQ=255)
     cAGCGGCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCcttc                                                                                                                                             <  2:771620/150‑1 (MQ=255)
     cAGCGGCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCcttc                                                                                                                                             <  2:2034650/150‑1 (MQ=255)
         ggCATGCTGCTGATGGCCTACAAAGTCTTCCGAACAGTACGTCGCACTTACCGCTCACAGTCTTAAGAGTCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCAGCAGGCctatgttc                                                                                                                                         <  2:1105011/150‑7 (MQ=255)
          gCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCAGCAGGCctatg                                                                                                                                            <  1:232410/146‑4 (MQ=255)
          gCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCAGCAGGCctatg                                                                                                                                            >  2:232410/1‑143 (MQ=255)
          gCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCAGCAGGCctat                                                                                                                                             >  1:2168479/1‑143 (MQ=255)
          gCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCAGCAGGCctat                                                                                                                                             <  2:2168479/145‑3 (MQ=255)
            aTGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATgg                                                                                                                                      <  2:798000/150‑1 (MQ=255)
             tgctgcTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATg                                                                                                                                       >  1:682034/1‑148 (MQ=255)
             tgctgcTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATg                                                                                                                                       <  2:682034/148‑1 (MQ=255)
              gctgctGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTc                                                                                                                                    <  2:1712839/150‑1 (MQ=255)
                tgctgATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCAt                                                                                                                                  <  2:2283867/150‑1 (MQ=255)
                tgctgATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTAGTGGGAGCGCAATGATGAAGCATGAAGCCTTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCAt                                                                                                                                  <  1:1993099/150‑1 (MQ=255)
                     aTGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCg                                                                                                                             <  2:1396791/150‑1 (MQ=255)
                        gCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGa                                                                                                                          <  2:972388/150‑1 (MQ=255)
                        gCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGa                                                                                                                          <  2:1848283/150‑1 (MQ=255)
                          cTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGc                                                                                                                        >  1:1033906/1‑150 (MQ=255)
                          cTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGc                                                                                                                        >  2:1365363/1‑150 (MQ=255)
                               aCGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCgg                                                                                                                   <  1:2021037/150‑1 (MQ=255)
                                     tCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGGCATCGCCGTGAGCGTCGGTGGCCt                                                                                                             <  1:65129/150‑1 (MQ=255)
                                       cgtaccgtacGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCttcttc                                                                                                                                          >  2:226056/1‑119 (MQ=255)
                                       cgtaccgtacGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCttcttc                                                                                                                                          <  1:226056/119‑1 (MQ=255)
                                            cgtacGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCg                                                                                                                       <  2:574422/133‑1 (MQ=255)
                                            cgtacGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCg                                                                                                                       >  1:574422/1‑133 (MQ=255)
                                            cgtacGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCa                                                                                                       >  2:330018/1‑149 (MQ=255)
                                            cgtacGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCa                                                                                                       <  1:330018/149‑1 (MQ=255)
                                             gtacGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGa                                                                                                     <  2:923881/150‑1 (MQ=255)
                                              tacGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCGGTGTGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGAt                                                                                                    >  1:623683/1‑150 (MQ=255)
                                                    ccgccAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTccc                                                                                              <  2:812548/150‑1 (MQ=255)
                                                     cgccAACCCGGCACAGGCTGAAGACGCCGCCAAGATCGTCGTTGTGGGAGCGCACTAATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCg                                                                                             >  1:1066216/1‑150 (MQ=255)
                                                      gccAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGa                                                                                                           <  2:794193/135‑1 (MQ=255)
                                                      gccAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGa                                                                                                           >  1:794193/1‑135 (MQ=255)
                                                       ccAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCt                                                                                           <  1:1522870/150‑1 (MQ=255)
                                                        cAACCCGGCACAGGCTGAAGAGGCCGCCAAGAGCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTg                                                                                          <  1:2088026/150‑1 (MQ=255)
                                                         aaCCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGt                                                                                         <  1:150240/150‑1 (MQ=255)
                                                            ccGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATg                                                                                                                                       >  1:1741221/1‑101 (MQ=255)
                                                            ccGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATg                                                                                                                                       <  2:1741221/101‑1 (MQ=255)
                                                            ccGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGtttt                                                                                      >  2:1479693/1‑150 (MQ=255)
                                                            ccGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGtttt                                                                                      >  1:2202758/1‑150 (MQ=255)
                                                             cGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGt                                                                                                                  >  2:2167584/1‑121 (MQ=255)
                                                             cGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGt                                                                                                                  <  1:2167584/121‑1 (MQ=255)
                                                             cGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGttttt                                                                                     <  2:478284/150‑1 (MQ=255)
                                                                 acaGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAg                                                                                 >  2:338304/1‑150 (MQ=255)
                                                                      cTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGt                                                                            >  1:1589500/1‑150 (MQ=255)
                                                                              gccgccAAGATCGTCGTTGTGGGGGCGCACTGATGCGTCATGCAGCCGTCGAGAAGAACATAGGCCTTCTGGCGTTCTTCATGGCCATTGCCGTGAGCGTCGGTGTCCTGACACAGATCGTCCGGCTGTTTTTCCAGGAGGGCACcaaca                                                                    >  2:201550/1‑150 (MQ=255)
                                                                              gccgccAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACcaaca                                                                    >  1:2141928/1‑150 (MQ=255)
                                                                              gccgccAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACcaaca                                                                    >  2:1493464/1‑150 (MQ=255)
                                                                                 gccAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGc                                                                 <  1:338304/150‑1 (MQ=255)
                                                                                 gccAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTcc                                                                                               <  1:1926384/119‑1 (MQ=255)
                                                                                 gccAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTcc                                                                                               >  2:1926384/1‑119 (MQ=255)
                                                                                   cAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTTGTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCg                                                               <  1:2250773/150‑1 (MQ=255)
                                                                                   cAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCa                                                                                                       >  1:1709235/1‑110 (MQ=255)
                                                                                   cAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCa                                                                                                       <  2:1709235/110‑1 (MQ=255)
                                                                                   cAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCg                                                               <  2:946430/150‑1 (MQ=255)
                                                                                          gtcgtTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCcttctt                                                                                                                                           <  2:701659/67‑1 (MQ=255)
                                                                                          gtcgtTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCcttctt                                                                                                                                           >  1:701659/1‑67 (MQ=255)
                                                                                          gtcgtTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAgg                                                        <  1:2059979/150‑1 (MQ=255)
                                                                                          gtcgtTGTGGGAGCGCACTGATGAAGCATGAAGCCGGCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAgg                                                        <  2:787711/150‑1 (MQ=255)
                                                                                             gtTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGGTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCa                                                     <  2:130733/150‑1 (MQ=255)
                                                                                              ttGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCAt                                                    >  1:957451/1‑150 (MQ=255)
                                                                                                  gggAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAg                                                >  1:1851892/1‑150 (MQ=255)
                                                                                                          aCTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATc                                                                                                                                 >  1:2271487/1‑61 (MQ=255)
                                                                                                          aCTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATc                                                                                                                                 <  2:2271487/61‑1 (MQ=255)
                                                                                                                gaagcatgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCt                                  >  1:140650/1‑150 (MQ=255)
                                                                                                                gaagcatgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGAATTATTATTTCCCGCTGTTTTTCCAGGACGTCACCAACAACCCGGTCGAGGGCATGAAGCCGCCCACCGCGCt                                  >  1:148299/1‑150 (MQ=255)
                                                                                                                  agcatgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTgg                                >  2:1038310/1‑150 (MQ=255)
                                                                                                                    catgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAg                              <  1:537589/150‑1 (MQ=255)
                                                                                                                     atgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGc                                                                                                               <  1:1465954/68‑1 (MQ=255)
                                                                                                                     atgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGc                                                                                                               >  2:1465954/1‑68 (MQ=255)
                                                                                                                     atgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACcaacaa                                                                   >  2:1630897/1‑112 (MQ=255)
                                                                                                                     atgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACcaacaa                                                                   <  1:1630897/112‑1 (MQ=255)
                                                                                                                        aagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAggg                                                       >  2:911905/1‑121 (MQ=255)
                                                                                                                        aagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAggg                                                       <  1:911905/121‑1 (MQ=255)
                                                                                                                         agcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGa                         >  1:28553/1‑150 (MQ=255)
                                                                                                                         agcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGa                         >  1:785024/1‑150 (MQ=255)
                                                                                                                          gcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGaa                        <  1:782058/150‑1 (MQ=255)
                                                                                                                                gagaAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAGGCcgc                  >  1:559115/1‑150 (MQ=255)
                                                                                                                                 agaagaACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAGGCcgcg                 >  1:2002641/1‑150 (MQ=255)
                                                                                                                                   aagaaCATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGcc                                              <  1:1250363/119‑1 (MQ=255)
                                                                                                                                   aagaaCATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAGGCCGCGat               >  1:480569/1‑150 (MQ=255)
                                                                                                                                   aagaaCATAGGCCTGCTGGCCTTCTTCAGGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGcc                                              >  2:1250363/1‑119 (MQ=255)
                                                                                                                                     gaaCATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAGGCCGCGatat             <  2:2141928/150‑1 (MQ=255)
                                                                                                                                     gaaCATAGGCCTGCTGGCCTTCCTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGt                            <  1:2294072/135‑1 (MQ=255)
                                                                                                                                     gaaCATAGGCCTGCTGGCCTTCCTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGt                            >  2:2294072/1‑135 (MQ=255)
                                                                                                                                         aTAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAg                                                                                 >  1:2026283/1‑78 (MQ=255)
                                                                                                                                         aTAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAg                                                                                 <  2:2026283/78‑1 (MQ=255)
                                                                                                                                             gCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCAAGATCGTCCCGCAGTTTTTCGAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAGGCCGCAATATCTACATcc     >  2:1329517/1‑150 (MQ=255)
                                                                                                                                               ctgctgGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAGGCCGCGATATCTACATCcgc   <  1:1493464/150‑1 (MQ=255)
                                                                                                                                                tgctgGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAGGCCGCGATATCTACATCcgcg  >  1:89115/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTCGCCAGCGGCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAGGCCGCGATATCTACATCCGCG  >  NC_002947/4841012‑4841305

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