Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 233,545 C→T G548G (GGC→GGT PP_0180 → cytochrome c family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947233,5450CT98.4% 193.0 / NA 63G548G (GGC→GGTPP_0180cytochrome c family protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (34/27);  minor base G (1/0);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCGAATTGTTCGAAGTGATCCTGTTCTACGAGACCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGCTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGCCGGTGTGCTGGGCACACGGCCAGTGGCGTTCTTT  >  NC_002947/233407‑233683
                                                                                                                                           |                                                                                                                                          
ccGCGAATTGTTCGAAGTGATCCTGTTCTACGAGACCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCt                                                                                                                                 >  2:1065623/1‑150 (MQ=255)
ccGCGAATTGTTCGAAGTGATCCTGTTCTACGAGACCCTGTGGCTGCAGGCAGGCCCCACCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCt                                                                                                                                 >  1:1152778/1‑150 (MQ=255)
  gcgAATTGTTCGAAGTGATCCTGTTCTACGAGACCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGc                                                                                                                               >  1:94538/1‑150 (MQ=255)
              aaGTGATCCTGTTCTACGAGACCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGt                                                                                                                   >  2:74839/1‑150 (MQ=255)
               aGTGATCCTGTTCTACGAGACCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCTCGACGGGCACGGTGGTGCTCGCTCCCCTGGCCTTGGTGAGTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGtt                                                                                                                  >  1:1123594/1‑150 (MQ=255)
                   aTCCTGTTCTACGAGACCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGAGGGGTTCGGCCAAGCTGCCGCTATCGCTGttcttc                                                                                                              >  1:26981/1‑150 (MQ=255)
                    tCCTGTTCTACGAGACCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCa                                                                                                             <  2:22888/150‑1 (MQ=255)
                      cTGTTCTACGAGACCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCg                                                                                                                                         <  2:1268007/120‑1 (MQ=255)
                      cTGTTCTACGAGACCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCg                                                                                                                                         >  1:1268007/1‑120 (MQ=255)
                       tGTTCTACGAGACCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCAGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCa                                                                                                          >  2:497510/1‑150 (MQ=255)
                        gTTCTACGAGACCCTGTGGCTGCAGGCAGGCCCCGCGGGGCATCAGGCGCTGCTGTCGGGTGGCGCGACGGCACGGGTGCTGTTGGTGGCCCAGGCCTGGGGGACTCGGAGGGTGTCTGCCAAGCTGCCGCTATCGCTGGTCTTCAGCAt                                                                                                         >  2:602937/1‑150 (MQ=255)
                            tACGAGACCCTGTGGCTGCAGGTAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAAc                                                                                                     >  1:707192/1‑150 (MQ=255)
                             aCGAGACCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTCGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAAcg                                                                                                    >  2:94143/1‑150 (MQ=255)
                               gagaCCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAAc                                                                                                     <  2:129115/147‑1 (MQ=255)
                               gagaCCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAAc                                                                                                     >  1:129115/1‑147 (MQ=255)
                                 gaCCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGATGGCACTGGTGCTGTTGGTGGGCCTGTCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAAcgccgc                                                                                                <  2:250161/150‑1 (MQ=255)
                                   cccTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCgcgc                                                                                              >  2:802612/1‑150 (MQ=255)
                                    ccTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCgct                                                                                             >  2:9210/1‑150 (MQ=255)
                                         ggCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCtgt                                                                                        <  1:74839/150‑1 (MQ=255)
                                           cTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTgc                                                                                      <  1:74926/150‑1 (MQ=255)
                                             gcaggtagGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTgcgc                                                                                    <  2:707192/150‑1 (MQ=255)
                                               aggcaggcCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTAGCCTGGGTGATTCTGCGGGGTTTGGCCAAGCTGCAGCTATCGCTGTTTTTCAGGATCAACGCCGTGCTGCTGTgcgcgc                                                                                  >  1:1121265/1‑150 (MQ=255)
                                                 gcaggcCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTg                                                                                >  2:1134532/1‑150 (MQ=255)
                                                 gcaggcCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTg                                                                                >  2:435154/1‑150 (MQ=255)
                                                    ggcCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCg                                                                             <  1:813643/150‑1 (MQ=255)
                                                      ccccGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCCGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGgt                                                                           >  2:944841/1‑150 (MQ=255)
                                                       cccGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGgta                                                                          >  1:388144/1‑150 (MQ=255)
                                                             tctcATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTc                                                                    <  1:791670/147‑1 (MQ=255)
                                                                 aTCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGCTTCTGCCGAGTTCGGCCAAGCTGGGGCTAGCGCGGTTCGTCAGCAGCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCgaga                                                                >  2:791670/1‑147 (MQ=255)
                                                                  tCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCgg                                                               <  2:1123594/150‑1 (MQ=255)
                                                                   cAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCggg                                                              >  2:897128/1‑150 (MQ=255)
                                                                    aGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGc                                                             >  1:1165431/1‑150 (MQ=255)
                                                                                    tGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGtt                                             <  2:1152778/150‑1 (MQ=255)
                                                                                        gcgACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCaa                                         <  1:1181541/150‑1 (MQ=255)
                                                                                          gACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAaga                                       <  1:1065623/150‑1 (MQ=255)
                                                                                          gACGGCACTGGTGCTGTCGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAaga                                       <  1:94143/150‑1 (MQ=255)
                                                                                              ggCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGc                                    >  2:1278403/1‑150 (MQ=255)
                                                                                                cACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAAcg                                                                                                    >  2:29560/1‑84 (MQ=255)
                                                                                                cACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAAcg                                                                                                    <  1:29560/84‑1 (MQ=255)
                                                                                                 aCTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGc                                                <  1:126515/135‑1 (MQ=255)
                                                                                                 aCTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGc                                                >  2:126515/1‑135 (MQ=255)
                                                                                                  cTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGGCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGCCGgt                                <  1:1134532/150‑1 (MQ=255)
                                                                                                   tGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGCCGgtg                               >  2:1156889/1‑150 (MQ=255)
                                                                                                   tGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGCTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGCCGgtg                               >  1:1076288/1‑150 (MQ=255)
                                                                                                      tGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGCCGGTGTGc                            >  1:467508/1‑150 (MQ=255)
                                                                                                      tGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCCGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGCCGGTGTGc                            <  1:944841/150‑1 (MQ=255)
                                                                                                       gCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGCCGGTGTGCt                           <  1:897128/150‑1 (MQ=255)
                                                                                                         tGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCaa                                         >  2:1185232/1‑134 (MQ=255)
                                                                                                         tGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCaa                                         <  1:1185232/134‑1 (MQ=255)
                                                                                                               tggGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTg                                                     >  2:860280/1‑116 (MQ=255)
                                                                                                               tggGCCTGGCCTGGGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTg                                                     <  1:860280/116‑1 (MQ=255)
                                                                                                                   cctggcctggGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGCCGGTGTGCTGGGCACACGGcc               >  1:188876/1‑150 (MQ=255)
                                                                                                                   cctggcctggGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGCCGGTGTGCTGGGCACACGGcc               >  1:827868/1‑150 (MQ=255)
                                                                                                                        cctggGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGCCGGTGTGCTGGGCACACGGCCAGTgg          <  2:388144/150‑1 (MQ=255)
                                                                                                                         ctggGTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGCCGGTGTGCTGGGCACACGGCCAGTGGc         <  2:1165431/150‑1 (MQ=255)
                                                                                                                           gggTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTc                                                                              <  1:1182412/79‑1 (MQ=255)
                                                                                                                           gggTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTc                                                                              >  2:1182412/1‑79 (MQ=255)
                                                                                                                            ggTGATTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGCCGGTGTGCTGGGCACACGGCCAGTGGCGtt      >  1:871981/1‑150 (MQ=255)
                                                                                                                                aTTCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGCCGGTGTGCTGGGCACACGGCCAGTGGCGTTCttt  <  1:435154/150‑1 (MQ=255)
                                                                                                                                 ttCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGAATGGCGTGAAGGCGTTGCAAGAGGCCGGTGTGCTGGGCACACGGCCAGt            <  1:133572/139‑1 (MQ=255)
                                                                                                                                 ttCTGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGAATGGCGTGAAGGCGTTGCAAGAGGCCGGTGTGCTGGGCACACGGCCAGt            >  2:133572/1‑139 (MQ=255)
                                                                                                                                    tGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCa                                                                                                             >  1:184497/1‑39 (MQ=255)
                                                                                                                                    tGCGGGGTTCGGCCAAGCTGCCGCTATCGCTGTTCTTCa                                                                                                             <  2:184497/39‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                          
CCGCGAATTGTTCGAAGTGATCCTGTTCTACGAGACCCTGTGGCTGCAGGCAGGCCCCGCCGGGCATCAGGCGGTGCTGGCGGGTGGCGCGACGGCACTGGTGCTGTTGGTGGGCCTGGCCTGGGTGATTCTGCGGGGCTCGGCCAAGCTGCCGCTATCGCTGTTCTTCAGCATCAACGCCGCGCTGCTGTGCGCGCTGTCGGTAGTATTCGCCGGGCATGGCGTGAAGGCGTTGCAAGAGGCCGGTGTGCTGGGCACACGGCCAGTGGCGTTCTTT  >  NC_002947/233407‑233683

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: