Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,841,156 T→A 5.1% L12Q (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,1560TA 5.1% 357.3 / 6.4 113L12Q (CTG→CAG) ccoO‑IIcbb3‑type cytochrome c oxidase subunit
Reads supporting (aligned to +/- strand):  ref base T (56/51);  new base A (2/4);  total (58/55)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

Only 100 of 115 total aligned reads displayed.

TTCTTCGCCAGCGGCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAGGCCGCGATATCTACATCCGCG  >  NC_002947/4841009‑4841305
                                                                                                                                                    |                                                                                                                                                      
ttcttcGCCAGCGGCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCctgc                                                                                                                                                      >  1:1787334/1‑150 (MQ=255)
        cAGCGGCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCcttc                                                                                                                                              >  1:1878276/1‑150 (MQ=255)
        cAGCGGCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCcttc                                                                                                                                              >  1:1411227/1‑150 (MQ=255)
         aGCGGCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCAGCAGGCctatg                                                                                                                                             <  2:1153048/150‑4 (MQ=255)
            ggCATGCGGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCAGCCAATTCAAGTTCTGGTGAAGATTCCCGCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCa                                                                                                                                          <  1:286238/150‑1 (MQ=255)
             gCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCAt                                                                                                                                         <  1:640915/150‑1 (MQ=255)
               aTGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATgg                                                                                                                                       <  2:288680/150‑1 (MQ=255)
                   tgctgATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCAt                                                                                                                                   >  1:1517103/1‑150 (MQ=255)
                         tGGCCTACAACGTGTTCCGTACAGTCATTTGCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGt                                                                                                                            <  1:611739/150‑1 (MQ=255)
                           gCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGa                                                                                                                          >  2:1215834/1‑150 (MQ=255)
                           gCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGa                                                                                                                          >  1:1792811/1‑150 (MQ=255)
                           gCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGGCAGCGCCGTGa                                                                                                                          <  2:1792811/150‑1 (MQ=255)
                             cTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGc                                                                                                                        >  1:1468770/1‑150 (MQ=255)
                             cTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGc                                                                                                                        >  1:1805997/1‑150 (MQ=255)
                                  aCGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCAGCAGGCctat                                                                                                                                              <  2:792363/124‑3 (MQ=255)
                                  aCGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCAGCAGGCctat                                                                                                                                              >  1:792363/1‑122 (MQ=255)
                                   cGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGt                                                                                                                  <  1:852321/150‑1 (MQ=255)
                                    ttgtTCCGGACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTg                                                                                                                 <  1:407165/149‑1 (MQ=255)
                                    gtgtTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTg                                                                                                                 >  1:1793194/1‑150 (MQ=255)
                                       ttCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGcc                                                                                                              <  1:1403913/150‑1 (MQ=255)
                                       ttCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGcc                                                                                                              >  2:1975364/1‑150 (MQ=255)
                                       ttCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGAATTCGCCATCGTGGTCGCCGTGAGCGTCGGTGGcc                                                                                                              >  2:2108755/1‑150 (MQ=255)
                                        tCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCt                                                                                                             <  2:1411227/150‑1 (MQ=255)
                                           gtaccgtacGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACGGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGAc                                                                                                          >  1:1045078/1‑150 (MQ=255)
                                              ccgtacGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCa                                                                                                       >  1:689800/1‑150 (MQ=255)
                                                gtacGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGa                                                                                                     <  1:907599/150‑1 (MQ=255)
                                                gtacGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGa                                                                                                     >  1:599791/1‑150 (MQ=255)
                                                 tacGTGCCGCCAACCCGGCACAGGCTGGAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTTCTTGCGTTGGGCGTCGGTGGCCTGACCCAGAt                                                                                                    >  1:139/1‑150 (MQ=255)
                                                     tGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATc                                                                                                   >  2:1362883/1‑147 (MQ=255)
                                                      gccgccAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTcc                                                                                               <  2:578942/150‑1 (MQ=255)
                                                          ccAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCt                                                                                           >  2:2207991/1‑150 (MQ=255)
                                                          ccAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGACGTGGAGAAGAAAATAGGCCTGCTAGGGTTATTCATTGTAATCGCGGGTAGAGTTGGTGGCTTGACCTAGATTGtaaagct                                                                                           >  2:814192/1‑144 (MQ=255)
                                                               cccGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGttt                                                                                       >  1:1045639/1‑150 (MQ=255)
                                                                         gCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACg                                                                             >  2:589637/1‑150 (MQ=255)
                                                                           tGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTc                                                                           >  1:165841/1‑150 (MQ=255)
                                                                             aaGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCAc                                                                         >  1:1259525/1‑150 (MQ=255)
                                                                                   ccgccAAGATCGTCGTTGTGGTAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGGAGCAGCTGCTGAGCGTTGTTGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACcaacaa                                                                   <  2:1418880/150‑1 (MQ=37)
                                                                                      ccAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGcc                                                                >  1:1078595/1‑150 (MQ=255)
                                                                                          gATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCATCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTc                                                            >  1:1379081/1‑150 (MQ=255)
                                                                                           aTCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGc                                                                                            >  2:2103075/1‑117 (MQ=255)
                                                                                           aTCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGc                                                                                            <  1:2103075/117‑1 (MQ=255)
                                                                                              gtcgtTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAgg                                                        <  1:817234/150‑1 (MQ=255)
                                                                                               tcgtTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCcttc                                                                                                                                              >  2:782094/1‑64 (MQ=255)
                                                                                               tcgtTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCcttc                                                                                                                                              <  1:782094/64‑1 (MQ=255)
                                                                                                 gtTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCa                                                     <  2:1259525/150‑1 (MQ=255)
                                                                                                     tgGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGaa                                                 >  2:327246/1‑150 (MQ=255)
                                                                                                      gggAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAg                                                >  2:866129/1‑150 (MQ=255)
                                                                                                      gggAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCAGGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAg                                                <  2:984044/150‑1 (MQ=255)
                                                                                                            gcACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGGTTTTCCAGGGCGTCACCAACAAGCCGGTCGAGGGCATGAAGCCgcgc                                          >  2:2115052/1‑150 (MQ=255)
                                                                                                             cACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCGGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCa                                         >  1:259714/1‑150 (MQ=255)
                                                                                                             cACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCgg                                                              >  1:1491390/1‑129 (MQ=255)
                                                                                                             cACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGaa                                                 >  1:878300/1‑142 (MQ=255)
                                                                                                             cACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGaa                                                 <  2:878300/142‑1 (MQ=255)
                                                                                                              aCTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGc                                                                                                               <  1:142335/79‑1 (MQ=255)
                                                                                                              aCTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGc                                                                                                               >  2:142335/1‑79 (MQ=255)
                                                                                                              aCTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCCCAc                                        >  2:642086/1‑150 (MQ=255)
                                                                                                              aCTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTGCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCAc                                        >  1:559663/1‑150 (MQ=255)
                                                                                                                tgatgaAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCt                                                                                                             <  2:2399691/79‑1 (MQ=255)
                                                                                                                tgatgaAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCt                                                                                                             >  1:2399691/1‑79 (MQ=255)
                                                                                                                  atgaagcatgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTg                                                                                                                 >  2:1522277/1‑73 (MQ=255)
                                                                                                                  atgaagcatgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTg                                                                                                                 <  1:1522277/73‑1 (MQ=255)
                                                                                                                  atgaagcatgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACcgcg                                    >  2:2260536/1‑150 (MQ=255)
                                                                                                                  atgaagcatgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACcgcg                                    >  1:1093255/1‑150 (MQ=255)
                                                                                                                    gaagcatgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGAAAGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCt                                  >  1:2098979/1‑150 (MQ=255)
                                                                                                                    gaagcatgaagcCGCCGCGACGGACGTCGTGCACCTGGCCTTGGCGGTTGTCAGCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCt                                  <  2:1468770/150‑1 (MQ=18)
                                                                                                                      agcatgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTgg                                >  1:930383/1‑150 (MQ=255)
                                                                                                                       gcatgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGTCATGAAGCCGCGCACCGCGCTGGa                               >  1:2188135/1‑150 (MQ=255)
                                                                                                                       gcatgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGa                               >  2:2225774/1‑150 (MQ=255)
                                                                                                                        catgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGttt                                                                                       >  1:1794164/1‑93 (MQ=255)
                                                                                                                        catgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGttt                                                                                       <  2:1794164/93‑1 (MQ=255)
                                                                                                                        catgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAg                              <  2:1705334/150‑1 (MQ=255)
                                                                                                                        catgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAg                              <  2:1805997/150‑1 (MQ=255)
                                                                                                                          tgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGt                            <  2:1787334/150‑1 (MQ=255)
                                                                                                                          tgaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCAGACTCAGAACGTTGCGCTCTTTATCCATCACGTCACCATCAAGCCTGTCGAGGTCATGAATCCGCGCACCGCGCTGGAGGt                            <  2:1480050/150‑1 (MQ=255)
                                                                                                                           gaagcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTg                           <  2:259714/150‑1 (MQ=255)
                                                                                                                              gcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGaa                        <  1:1719569/150‑1 (MQ=255)
                                                                                                                              gcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGaa                        <  2:1045639/150‑1 (MQ=255)
                                                                                                                              gcCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGaa                        <  2:1045078/150‑1 (MQ=255)
                                                                                                                               ccGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGCACGTCACCAATAAGCCGGTCGAGGGCATGAAGCCGCGCGTCGAGCTTGAGGTgtatt                       <  2:599791/150‑5 (MQ=25)
                                                                                                                                  tCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGa                                                          <  1:1216332/112‑1 (MQ=255)
                                                                                                                                  tCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGa                                                          >  2:1216332/1‑112 (MQ=255)
                                                                                                                                  tCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAGGcc                    >  2:298523/1‑150 (MQ=255)
                                                                                                                                    gagaAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAGGCcgc                  <  1:2108755/150‑1 (MQ=255)
                                                                                                                                     agaagaACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAGGCcgcg                 <  1:327246/150‑1 (MQ=255)
                                                                                                                                      gaagaaCATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCAt                                                    >  1:560552/1‑114 (MQ=255)
                                                                                                                                      gaagaaCATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCAt                                                    <  2:560552/114‑1 (MQ=255)
                                                                                                                                      gaagaaCATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAGGCCGCGa                <  2:1359515/150‑1 (MQ=255)
                                                                                                                                       aagaaCATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCa                                                     >  1:667344/1‑112 (MQ=255)
                                                                                                                                       aagaaCATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCa                                                     <  2:667344/112‑1 (MQ=255)
                                                                                                                                       aagaaCATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGAAAGGCCGCGat               >  2:210239/1‑150 (MQ=255)
                                                                                                                                         gaaCATAGGCCTGCTGGCCTTCTTCCTGGTCAGCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAGGCCGCGatat             >  1:1271497/1‑150 (MQ=255)
                                                                                                                                         gaaCATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGaa                        >  1:141034/1‑139 (MQ=255)
                                                                                                                                         gaaCATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGaa                        <  2:141034/139‑1 (MQ=255)
                                                                                                                                           aCATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAg                                                         >  1:312492/1‑104 (MQ=255)
                                                                                                                                            agaaGGCCAGCAGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAg                       <  2:1480520/134‑1 (MQ=255)
                                                                                                                                            agaaGGCCAGCAGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAg                       >  1:1480520/4‑137 (MQ=255)
                                                                                                                                                ggCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCgcgc                                   <  2:189894/121‑1 (MQ=255)
                                                                                                                                                 gCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTTACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGaa                        >  2:2145342/1‑131 (MQ=255)
                                                                                                                                                 gCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGaa                        <  1:2145342/131‑1 (MQ=255)
                                                                                                                                                    tgctgGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAGGCCGCGATATCTACATCcgcg  >  1:1554033/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TTCTTCGCCAGCGGCATGCTGCTGATGGCCTACAACGTGTTCCGTACCGTACGTGCCGCCAACCCGGCACAGGCTGAAGAGGCCGCCAAGATCGTCGTTGTGGGAGCGCACTGATGAAGCATGAAGCCGTCGAGAAGAACATAGGCCTGCTGGCCTTCTTCATGGTCATCGCCGTGAGCGTCGGTGGCCTGACCCAGATCGTCCCGCTGTTTTTCCAGGACGTCACCAACAAGCCGGTCGAGGGCATGAAGCCGCGCACCGCGCTGGAGGTGGAAGGCCGCGATATCTACATCCGCG  >  NC_002947/4841009‑4841305

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