Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,282,780 T→G 9.7% intergenic (‑129/+100) PP_2010 ← / ← PP_2011 cytochrome b561/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,282,7800TG 9.7% 189.6 / 13.1 72intergenic (‑129/+100)PP_2010/PP_2011cytochrome b561/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (26/39);  new base G (3/4);  total (29/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

TGGGGGGGGATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTG  >  NC_002947/2282633‑2282840
                                                                                                                                                   |                                                            
gggggggggATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCttg                                                            <  2:2204263/149‑3 (MQ=255)
 ggggggggATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGc                                                           <  1:488741/150‑1 (MQ=255)
 ggggggggATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGc                                                           <  2:1118029/150‑1 (MQ=255)
        gATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCg                                                    >  2:1460668/1‑150 (MQ=255)
             aaCCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTccc                                               >  2:883482/1‑150 (MQ=255)
             aaCCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTccc                                               >  1:890327/1‑150 (MQ=255)
              aCCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCca                                              >  1:2129749/1‑150 (MQ=255)
              aCCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCca                                              >  1:2129604/1‑150 (MQ=255)
                  tGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGt                                          <  2:1457787/150‑1 (MQ=255)
                   gATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTc                                         >  2:2528361/1‑150 (MQ=255)
                   gATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTc                                         <  2:313887/150‑1 (MQ=255)
                       ggTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAAc                                     <  1:2528361/150‑1 (MQ=255)
                           cTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGca                                 >  2:941613/1‑150 (MQ=255)
                            tCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGcac                                >  1:2883623/1‑150 (MQ=255)
                             cGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGcaca                               <  1:158273/150‑1 (MQ=255)
                             cGTAATCAGGTCTGAATGGCGCTCTTGCTGGCGGTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGcaca                               <  1:790810/150‑1 (MQ=255)
                                   caGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTg                         <  1:1095201/150‑1 (MQ=255)
                                   caGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTg                         <  1:2060557/150‑1 (MQ=255)
                                     ggTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGa                       <  2:2883623/150‑1 (MQ=255)
                                     ggTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTCGCCTGCTCCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttccaa                      <  2:1553897/150‑5 (MQ=21)
                                         tGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttccaagac                   >  1:92806/1‑143 (MQ=255)
                                           aaTGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCt                                                              <  2:2035241/105‑1 (MQ=255)
                                           aaTGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCt                                                              >  1:2035241/1‑105 (MQ=255)
                                               gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCa             <  2:253235/150‑1 (MQ=255)
                                               gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCa             <  2:2129749/150‑1 (MQ=255)
                                               gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCa             <  2:2129604/150‑1 (MQ=255)
                                               gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa             <  1:883482/150‑1 (MQ=38)
                                                cgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAg            <  2:92806/150‑1 (MQ=37)
                                                cgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAg            <  1:941613/150‑1 (MQ=37)
                                                 gcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGt           >  1:1061562/1‑150 (MQ=255)
                                                 gcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGt           <  2:2781881/150‑1 (MQ=37)
                                                  ctctTGCTGGCGCTTAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt          <  2:1755782/150‑1 (MQ=255)
                                                   tctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTg         <  2:1061562/150‑1 (MQ=255)
                                                     ttGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGa       <  2:890327/150‑1 (MQ=255)
                                                     ttGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCAAAGTCTTGGAAGCCTGTGCAGTTGGa       <  1:1628883/150‑1 (MQ=255)
                                                       gCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa             >  1:745494/1‑142 (MQ=255)
                                                       gCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa             <  2:745494/142‑1 (MQ=255)
                                                          ggCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACtgtg  <  1:1937329/150‑1 (MQ=255)
                                                          ggCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACtgtg  <  1:1460668/150‑1 (MQ=255)
                                                           gcgcTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAg            <  1:721233/139‑1 (MQ=21)
                                                           gcgcTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAg            >  2:721233/1‑139 (MQ=21)
                                                             gcTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGc                    >  1:900961/1‑129 (MQ=255)
                                                             gcTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGc                    <  2:900961/129‑1 (MQ=255)
                                                                    aGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTcc                                                <  1:1568568/94‑1 (MQ=255)
                                                                    aGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTcc                                                >  2:1568568/1‑94 (MQ=255)
                                                                         ttttAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGc              >  2:2467701/1‑123 (MQ=255)
                                                                         ttttAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGc              <  1:2467701/123‑1 (MQ=255)
                                                                             aaTGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGca                                 >  1:2278784/1‑100 (MQ=255)
                                                                             aaTGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGca                                 <  2:2278784/100‑1 (MQ=255)
                                                                                 gCTACCGGAAACAGATCTTTGCGCCAGGTCTATGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgt                <  2:2071075/113‑1 (MQ=255)
                                                                                 gCTACCGGAAACAGATCTTTGCGCCAGGTCTATGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgt                >  1:2071075/1‑113 (MQ=255)
                                                                                      cGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCaa                                      <  2:762647/86‑1 (MQ=255)
                                                                                      cGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCaa                                      >  1:762647/1‑86 (MQ=255)
                                                                                      cGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACtgt   >  2:1824899/1‑121 (MQ=255)
                                                                                      cGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACtgt   <  1:1824899/121‑1 (MQ=255)
                                                                                           aCAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgtg               >  1:440437/1‑104 (MQ=255)
                                                                                           aCAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgtg               <  2:440437/104‑1 (MQ=255)
                                                                                             aGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc                        <  1:54263/93‑3 (MQ=255)
                                                                                             aGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc                        >  2:54263/1‑91 (MQ=255)
                                                                                                         cAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt          >  1:2423187/1‑95 (MQ=255)
                                                                                                         cAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt          <  2:2423187/95‑1 (MQ=255)
                                                                                                               tgtgAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGc                    <  2:204289/79‑1 (MQ=255)
                                                                                                               tgtgAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGc                    >  1:204289/1‑79 (MQ=255)
                                                                                                                 tgAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgtg               >  2:2927520/1‑82 (MQ=255)
                                                                                                                 tgAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgtg               <  1:2927520/82‑1 (MQ=255)
                                                                                                                        tgtgTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa             <  2:2852409/77‑1 (MQ=14)
                                                                                                                        tgtgTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa             >  1:2852409/1‑77 (MQ=14)
                                                                                                                                    ctccGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACtgtg  >  1:187419/3‑76 (MQ=14)
                                                                                                                                    ctccGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACtgtg  <  2:187419/74‑1 (MQ=14)
                                                                                                                                    caCCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc                        >  1:329844/1‑52 (MQ=255)
                                                                                                                                    caCCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc                        <  2:329844/54‑3 (MQ=11)
                                                                                                                                       cGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCtt                          >  2:1297409/1‑49 (MQ=255)
                                                                                                                                       cGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCtt                          <  1:1297409/49‑1 (MQ=255)
                                                                                                                                              ttCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgt                <  2:2687802/52‑1 (MQ=255)
                                                                                                                                              ttCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgt                >  1:2687802/1‑52 (MQ=255)
                                                                                                                                                   tGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTgg        <  1:191957/55‑1 (MQ=255)
                                                                                                                                                   tGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTgg        >  2:191957/1‑55 (MQ=255)
                                                                                                                                                   |                                                            
TGGGGGGGGATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTG  >  NC_002947/2282633‑2282840

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