Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,284,695 T→C 10.9% intergenic (‑122/+107) 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,284,6950TC10.9% 254.8 / 16.8 94intergenic (‑122/+107)PP_2010/PP_2011cytochrome b561/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (6/77);  new base C (3/7);  total (9/85)
Fisher's exact test for biased strand distribution p-value = 5.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

AATGGGTGGTGGGGGGGGATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAG  >  NC_002947/2284546‑2284826
                                                                                                                                                     |                                                                                                                                   
actGGGTGGTGGGGGGGGATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCct                                                                                                                                     <  1:3662135/148‑1 (MQ=255)
aaTGGGTGGTGGGGGGGGATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCct                                                                                                                                     <  2:1029698/150‑1 (MQ=255)
aaTGGGTGGTGGGGGGGGATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCct                                                                                                                                     >  2:3200433/1‑150 (MQ=255)
  tGGGTGGTGGGGGGGGATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCctct                                                                                                                                   >  2:2023948/1‑150 (MQ=255)
   gggTGGTGGGGGGGGATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCtt                                                                                                                                  <  2:3095213/150‑1 (MQ=255)
       ggtggGGGGGGATGAAACCATGAGTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCt                                                                                                                              <  2:4761948/150‑1 (MQ=255)
        gtggGGGGGGATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTg                                                                                                                             <  2:3903490/150‑1 (MQ=255)
        gtggGGGGGGATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTg                                                                                                                             <  2:893170/150‑1 (MQ=255)
        gtggGGGGGGATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTg                                                                                                                             <  1:2165379/150‑1 (MQ=255)
               gggATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGcc                                                                                                                      <  2:4800870/150‑1 (MQ=255)
                    gAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTc                                                                                                                 <  2:357332/150‑1 (MQ=255)
                    gAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTc                                                                                                                 <  1:3908227/150‑1 (MQ=255)
                     aaaCCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTcc                                                                                                                <  1:4047594/150‑1 (MQ=255)
                       aCCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCca                                                                                                              <  1:286016/150‑1 (MQ=255)
                       aCCATGAGTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCca                                                                                                              <  1:4299112/150‑1 (MQ=255)
                         cctGAGTGGGCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTGTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCcaca                                                                                                            <  1:4901982/148‑1 (MQ=255)
                         cATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTCGCCTGCTCCGGCCCCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                           >  2:4259710/1‑150 (MQ=38)
                          aTGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAg                                                                                                           <  1:1464749/150‑1 (MQ=255)
                          aTGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTCGCCTGCTCCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGt                                                                                                          <  1:512335/150‑1 (MQ=38)
                          aTGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTCGCCTGCTCCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGt                                                                                                          <  1:512340/150‑1 (MQ=38)
                            gATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCTCTCTGGAGTCATTTTAATGGCTACCTGATAGAGATCTTTGCGGCAGGTCTGTGAAAATTGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTc                                                                                                         <  2:4152335/150‑1 (MQ=255)
                            gATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCCCTTCGCTGGCACGCCCGCTCCCACAGTc                                                                                                         <  2:3376958/150‑1 (MQ=255)
                             aTTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTcc                                                                                                        <  2:5045936/150‑1 (MQ=255)
                                    cTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGca                                                                                                 <  2:856790/150‑1 (MQ=255)
                                    cTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGca                                                                                                 <  2:339078/150‑1 (MQ=255)
                                      cGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGcaca                                                                                               <  2:1555840/150‑1 (MQ=255)
                                         catcaGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTc                                                                                            <  1:2606243/150‑1 (MQ=255)
                                           tcaGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCtt                                                                                          <  2:2411769/150‑1 (MQ=255)
                                           tcaGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCtt                                                                                          <  1:3829131/150‑1 (MQ=255)
                                           tcaGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCtt                                                                                          <  1:1248736/150‑1 (MQ=255)
                                           tcaGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCtt                                                                                          <  1:4735116/150‑1 (MQ=255)
                                           tcaGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCtt                                                                                          <  1:4733303/150‑1 (MQ=255)
                                           tcaGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACACTCCAACTGCACAGTcgt                                                                                          >  2:3847167/1‑148 (MQ=255)
                                           tcaGGTCTCAATCGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCtt                                                                                          <  2:3993512/150‑1 (MQ=255)
                                            caGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTg                                                                                         <  2:2105341/150‑1 (MQ=255)
                                            caGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTg                                                                                         <  1:3200433/150‑1 (MQ=255)
                                            caGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTg                                                                                         <  2:3631744/150‑1 (MQ=255)
                                            caGGTCTGAATGGCGCTCTGGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTg                                                                                         <  1:3124879/150‑1 (MQ=255)
                                             cggcttGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc                                                                                        <  2:1436477/145‑3 (MQ=255)
                                             aGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCttcc                                                                                        <  2:2153092/150‑3 (MQ=255)
                                             aGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTgg                                                                                        <  1:1118313/150‑1 (MQ=255)
                                             aGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc                                                                                        <  2:1849468/150‑3 (MQ=255)
                                             aGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc                                                                                        <  2:3075580/150‑3 (MQ=255)
                                             aGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc                                                                                        <  2:3005811/150‑3 (MQ=255)
                                             aGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc                                                                                        <  1:3215924/150‑3 (MQ=255)
                                             aGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc                                                                                        <  2:4721877/150‑3 (MQ=255)
                                             aGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc                                                                                        <  1:4891332/150‑3 (MQ=255)
                                             aGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc                                                                                        <  1:59129/150‑3 (MQ=255)
                                             aGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCGGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc                                                                                        <  1:5018359/150‑3 (MQ=255)
                                             aGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTGTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTgg                                                                                        <  1:505058/150‑1 (MQ=255)
                                               gTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGaa                                                                                      <  1:2666276/150‑1 (MQ=255)
                                               gTCTGAATGGCGCTCTTGCTGGCACTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGaa                                                                                      <  1:1954932/150‑1 (MQ=255)
                                                tCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAg                                                                                     >  1:3434417/1‑150 (MQ=255)
                                                      tGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgtg                                                                               <  1:2218799/150‑1 (MQ=255)
                                                      tGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTCGCCTGCTCCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGc                                                                              <  1:4259710/150‑1 (MQ=37)
                                                      gggCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACtgtg                                                                               <  2:4002085/149‑1 (MQ=255)
                                                       ggCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGc                                                                              <  2:2164939/150‑1 (MQ=255)
                                                       ggCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTCGCCTGCTCCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa                                                                             <  1:2245052/150‑1 (MQ=18)
                                                        gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCa                                                                             <  2:3755361/150‑1 (MQ=255)
                                                        gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa                                                                             <  2:781494/150‑1 (MQ=255)
                                                        gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa                                                                             <  2:1686706/150‑1 (MQ=37)
                                                        gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa                                                                             <  1:399964/150‑1 (MQ=37)
                                                        gcgcTCGTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCa                                                                             <  1:494478/150‑1 (MQ=255)
                                                         cgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAg                                                                            <  1:1898935/150‑1 (MQ=255)
                                                         cgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAg                                                                            <  1:1242111/150‑1 (MQ=255)
                                                          gcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGt                                                                           <  1:1613242/150‑1 (MQ=255)
                                                          gcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGt                                                                           <  2:3675467/150‑1 (MQ=255)
                                                           ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGTAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                          <  1:1981390/150‑1 (MQ=255)
                                                           ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                          <  2:2631138/150‑1 (MQ=255)
                                                           ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                          <  2:4765325/150‑1 (MQ=255)
                                                           ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                          <  1:4876153/150‑1 (MQ=255)
                                                           ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                          <  2:3837264/150‑1 (MQ=255)
                                                           ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                          <  2:2038781/150‑1 (MQ=255)
                                                           ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGtt                                                                          <  2:3448649/150‑1 (MQ=255)
                                                           ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGtt                                                                          <  2:524996/150‑1 (MQ=37)
                                                           ctcgTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                          <  2:1387012/150‑1 (MQ=255)
                                                            tctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTg                                                                         <  1:757568/150‑1 (MQ=255)
                                                                  tGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACtgt                                                                   <  1:4100986/150‑1 (MQ=25)
                                                                   ggCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACtgtg                                                                  <  1:3760294/150‑1 (MQ=255)
                                                                    gcgcTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTgg                                                                 <  2:2391235/150‑1 (MQ=255)
                                                                      gcTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACACTCCAACAGCACAGTCGTGGAAGCCGGTGCAGGTGGACTGTgggc                                                               >  2:610975/1‑149 (MQ=255)
                                                                       cTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAg                                                                            >  2:635069/1‑136 (MQ=25)
                                                                       cTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAg                                                                            <  1:635069/136‑1 (MQ=25)
                                                                        tCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGc                                                             <  2:3434417/150‑1 (MQ=255)
                                                                                    ttAATTGCTACCGTAAACACATCTTTGCGCCAGGTCTGTCAAAAATGTGTTGCCACCCACCGGCCTCTTCGCTGGCACGCCCGACCCAACAGTCAAACTGCACAGTATTGGAAGCATGTGCAGTTGGACTGTGGTAGCGGGCGTGCCcgc                                                 <  2:925234/150‑1 (MQ=255)
                                                                                         ggCTCCCGGAAACAGATCTTGGCGCCAGGTCTGTGAAAAATGTGTGGCCAGCCACCGGCCACTTCGCGGGCACGCCCGCTCCCACCGCCCACCTGCCCCGCCTTGGAAGCCTGGGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAgg                                            <  1:1889614/150‑1 (MQ=21)
                                                                                         ggCTACCGTAAACATATCTTTGCGCCAGGTCTGTTAAAAATGTGTTGCCAGCCCCCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAgg                                            <  1:1889559/150‑1 (MQ=255)
                                                                                         gccTACCGTAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAgg                                            <  1:3699419/148‑1 (MQ=255)
                                                                                            tACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGTAGCCTGTGCAGTTGGACTGGGGGTGCGGGCGTGCCCGCGAAGGGGc                                         >  2:4340313/1‑150 (MQ=255)
                                                                                                 gAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCGTCGCTGGCACGCCCGCTCCCACAGTCCACCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAg                                    <  2:1682701/150‑1 (MQ=255)
                                                                                                      aGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCCTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCAGAGTGTTGGAAGCCTGTGGAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGGCGGAGCAGGc                               >  2:3915857/1‑150 (MQ=37)
                                                                                                           tttGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACACTCTTGGAAGCCTGTGCAGTTGGCCTGTGGGAGCGGGCGTGCCCGCGGAGGGGCCGGAGCAGGCGACag                          >  2:295977/1‑149 (MQ=255)
                                                                                                             tGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCa                        <  2:3367501/150‑1 (MQ=255)
                                                                                                              gcgcCAGTTCTGTCAAAAATGTGTCGCCAGCCACCGGCCTCTTCGCTGCCACGCCCGCTCCCACAGTCAAACGGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAg                       <  1:4340313/150‑1 (MQ=25)
                                                                                                                gcCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCCCTTCGCGGGCACGccc                                                                                                                     >  1:3439898/1‑54 (MQ=37)
                                                                                                                gcCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCCCTTCGCGGGCACGccc                                                                                                                     <  2:3439898/54‑1 (MQ=37)
                                                                                                                                  gtgtTGCCCGCCACGGGCCTCTTCGCTGGCACGCCCGCCCCCACCGTCCACCTGAACAGTCTTGTATGCCTGTGCTGTGGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCa   <  1:3847167/150‑1 (MQ=25)
                                                                                                                                   tgtTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTATTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAg  <  1:2023948/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                   
AATGGGTGGTGGGGGGGGATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAG  >  NC_002947/2284546‑2284826

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