Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,911,824 C→G 5.8% Q134H (CAG→CAC ccmH ← cytochrome c heme lyase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,911,8240CG 5.8% 220.6 / 4.9 69Q134H (CAG→CACccmHcytochrome c heme lyase subunit
Reads supporting (aligned to +/- strand):  ref base C (33/32);  new base G (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

GCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCA  >  NC_002947/4911684‑4911972
                                                                                                                                            |                                                                                                                                                     
gCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGTCGCCCTGCACGGcc                                                                                                                                             >  2:40732/1‑150 (MQ=255)
gCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGcc                                                                                                                                             >  2:696616/1‑150 (MQ=255)
gCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGcc                                                                                                                                             <  2:901849/150‑1 (MQ=255)
 cAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGAGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCg                                                                                                                                            <  1:949315/150‑1 (MQ=255)
      tCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCa                                                                                                                                       >  1:1505332/1‑150 (MQ=255)
       cGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCAc                                                                                                                                      <  2:1281929/150‑1 (MQ=255)
       cGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCAc                                                                                                                                      <  2:1624586/150‑1 (MQ=255)
        gggAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACg                                                                                                                                     <  1:305820/150‑1 (MQ=255)
         ggAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGCCGGCACGCTCGCGTTCTTAGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACgg                                                                                                                                    >  2:387452/1‑150 (MQ=255)
                   aaaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCCTCCGCAGACAACTCATCGGCGCTCCGTTCGGCCGTAGCACGCTGCCTGCGtg                                                                                                                          <  2:602443/150‑2 (MQ=255)
                          aGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCAc                                                                                                                   <  1:1256502/150‑1 (MQ=255)
                           gggCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCAcc                                                                                                                  >  1:221323/1‑150 (MQ=255)
                             gCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGc                                                                                                                >  2:224327/1‑150 (MQ=255)
                             gCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTCTATCCCGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGc                                                                                                                >  2:222563/1‑150 (MQ=255)
                               aGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCa                                                                                                              >  1:113424/1‑150 (MQ=255)
                                  gAGCAGCAGAAGGCCCGCGCTAAGCCAGAAGTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGca                                                                                                           >  1:965084/1‑150 (MQ=255)
                                       gcagAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACg                                                                                                      >  1:996709/1‑150 (MQ=255)
                                        cagAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGTGCAGGGCGcc                                                                                                                                >  1:716211/1‑123 (MQ=38)
                                        cagAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGTGCAGGGCGcc                                                                                                                                <  2:716211/123‑1 (MQ=38)
                                         agAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGaa                                                                                                    <  1:696616/150‑1 (MQ=255)
                                                  gcgcTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCACCATCGTACGATCACCGCCAGCACCACGAAGCCACCGGc                                                                                           >  1:1577486/1‑150 (MQ=255)
                                                          ccAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCAc                                                                                               <  2:293421/138‑1 (MQ=255)
                                                          ccAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCAc                                                                                               >  1:293421/1‑138 (MQ=255)
                                                             gAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATcc                                                                                >  2:1331501/1‑150 (MQ=255)
                                                               aTTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATcccc                                                                              <  2:1624594/150‑1 (MQ=255)
                                                                   aaTCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAACCGCCCGGCCATCAAGATCCCCGGGc                                                                          <  2:965084/150‑1 (MQ=255)
                                                                    aTCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATc                                                                                 <  2:889859/142‑1 (MQ=255)
                                                                    aTCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATc                                                                                 >  1:889859/1‑142 (MQ=255)
                                                                    aTCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGcc                                                                         >  2:240419/1‑150 (MQ=255)
                                                                     tCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCa                                                                        >  1:796516/1‑150 (MQ=255)
                                                                           gTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCa                                                                                                              <  2:1525730/106‑1 (MQ=255)
                                                                           gTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCa                                                                                                              >  1:1525730/1‑106 (MQ=255)
                                                                              tGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAg                                                                                                                                         >  2:564218/1‑76 (MQ=255)
                                                                              tGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAg                                                                                                                                         <  1:564218/76‑1 (MQ=255)
                                                                               gTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagca                                                              >  1:1120893/1‑150 (MQ=255)
                                                                               gTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagca                                                              >  2:534091/1‑150 (MQ=255)
                                                                                ttCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGcc                                                                                                                                             <  1:1258837/70‑1 (MQ=255)
                                                                                ttCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGcc                                                                                                                                             >  2:1258837/1‑70 (MQ=255)
                                                                                ttCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagcag                                                             >  1:91318/1‑150 (MQ=255)
                                                                                  cTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGcc                                                           >  1:759028/1‑150 (MQ=255)
                                                                                   tttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAAGTCAGCGGCGCCCTGCACGGCCGTAGCAAGGCGCCGGCGTACGATCACCGCCAGCACGACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCa                                                          >  1:785466/1‑150 (MQ=255)
                                                                                   tttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAg                                                                                                             >  2:938115/1‑99 (MQ=255)
                                                                                   tttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAg                                                                                                             <  1:938115/99‑1 (MQ=255)
                                                                                   tttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagca                                                              <  1:1039759/146‑1 (MQ=255)
                                                                                   tttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagca                                                              >  2:1039759/1‑146 (MQ=255)
                                                                                      ttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCCTAGCACGGCGCCGGCGTCCGATCACCTCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGt                                                       <  2:108013/150‑1 (MQ=255)
                                                                                         cagcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAg                                                                                        >  2:95381/1‑114 (MQ=255)
                                                                                         cagcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAg                                                                                        <  1:95381/114‑1 (MQ=255)
                                                                                         cagcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCaa                                                                                     >  1:1092514/1‑117 (MQ=255)
                                                                                         cagcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCaa                                                                                     <  2:1092514/117‑1 (MQ=255)
                                                                                                     gagaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGc                                                                                           <  2:917805/99‑1 (MQ=255)
                                                                                                     gagaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGc                                                                                           >  1:917805/1‑99 (MQ=255)
                                                                                                     gagaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCt                                              <  1:604949/144‑1 (MQ=255)
                                                                                                     gagaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCt                                              >  2:604949/1‑144 (MQ=255)
                                                                                                            tCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTg                                 <  2:104397/150‑1 (MQ=255)
                                                                                                              gcgTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTa                               <  2:796516/150‑1 (MQ=255)
                                                                                                                gTTCTTCGGCAGACAACTCATCGGCGACCTGCCCGGCCCTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAgc                             <  2:1365095/150‑1 (MQ=255)
                                                                                                                          aGACAACTCATCGGCGCCCTGCACGGGCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTc                   <  2:785466/150‑1 (MQ=255)
                                                                                                                          aGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAACCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTGGCGCACGAAGTc                   <  2:221323/150‑1 (MQ=255)
                                                                                                                           gACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagca                                                              <  2:959020/106‑1 (MQ=255)
                                                                                                                           gACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagca                                                              >  1:959020/1‑106 (MQ=255)
                                                                                                                              aaCTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCa               >  2:1463699/1‑150 (MQ=255)
                                                                                                                                cTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATa             >  2:439185/1‑150 (MQ=255)
                                                                                                                                 tCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAg            <  2:996709/150‑1 (MQ=255)
                                                                                                                                    tCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCgg         <  1:439185/150‑1 (MQ=255)
                                                                                                                                     cGGCGCCGTGCTACGGCCGTGCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCgg                                                                            >  2:606370/1‑82 (MQ=25)
                                                                                                                                     cGGCGCCGTGCTACGGCCGTGCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCgg                                                                            <  1:606370/82‑1 (MQ=25)
                                                                                                                                     cGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGt        <  1:534091/150‑1 (MQ=255)
                                                                                                                                         gcCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGAc    <  2:1505332/150‑1 (MQ=255)
                                                                                                                                           ccTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCa  >  1:953857/1‑150 (MQ=255)
                                                                                                                                           ccTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCACCAAGGCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTAGACCa  >  1:928059/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
GCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCA  >  NC_002947/4911684‑4911972

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