Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,196,274 Δ1 bp 7.7% coding (329/1755 nt) xcpQ → type II secretion pathway protein XcpQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,196,2740C. 7.7% 199.2 / 11.2 65coding (329/1755 nt)xcpQtype II secretion pathway protein XcpQ
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base . (3/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

CCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATCCTCAAACCGCTGATCGACCCGCGAGTTGGGG  >  NC_002947/1196125‑1196423
                                                                                                                                                      |                                                                                                                                                       
cccccGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGc                                                                                                                                                         >  2:1129588/1‑150 (MQ=255)
  cccGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCcg                                                                                                                                                      >  2:777373/1‑150 (MQ=255)
   gcGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAACATCGTGCCCAACCAGGCcgc                                                                                                                                                     >  2:1079052/2‑150 (MQ=255)
   ccGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGcgcg                                                                                                                                                      >  1:255580/1‑147 (MQ=255)
          aaGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTg                                                                                                                                             >  1:1185280/1‑150 (MQ=255)
             ggCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGAAGATATACTTCTATTGGTAAGTGATGCTCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCGCGCGGCCTggtc                                                                                                                                          <  1:1236336/150‑3 (MQ=255)
               cAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCGCGCGGCCTggtcgg                                                                                                                                        <  2:255580/150‑5 (MQ=255)
                    tCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGc                                                                                                                                   <  1:1318804/150‑1 (MQ=255)
                      aCCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCa                                                                                                                                 >  2:1313395/1‑150 (MQ=255)
                      aCCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCa                                                                                                                                 <  1:777373/150‑1 (MQ=255)
                        cGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGGAAATCGTGCCCGACCAGGCCGCGCGCCTGGAgccggtgccggt                                                                                                                               >  2:477059/1‑150 (MQ=255)
                           gcgcACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGa                                                                                                                            <  1:999277/150‑1 (MQ=255)
                            tgcACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAg                                                                                                                           <  1:417444/149‑1 (MQ=255)
                                   tCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAggcgggcg                                                                                                                    <  2:745441/150‑1 (MQ=255)
                                     gACCTGCACAACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGcgggcggg                                                                                                                  <  1:1129588/150‑1 (MQ=255)
                                      aCCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGTCAGTGGAGGCgggcgggc                                                                                                                 >  2:657285/1‑150 (MQ=255)
                                      aCCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGAGGTCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCgggcgggc                                                                                                                 >  1:969171/1‑150 (MQ=255)
                                       ccTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGgca                                                                                                                >  1:926747/1‑150 (MQ=255)
                                          gCACGACGCCGACGCAATCTACTACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGgcagca                                                                                                             <  2:969171/150‑1 (MQ=255)
                                                        gccATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGc                                                                                                <  2:767760/150‑1 (MQ=255)
                                                           aTCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGAc                                                                                             >  2:1283223/1‑150 (MQ=255)
                                                            tCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAAGCGGACGGGCAGCTGGGCGAGGGCAGCGACa                                                                                            >  2:330940/1‑150 (MQ=255)
                                                                cGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCgg                                                                                        >  2:806416/1‑149 (MQ=255)
                                                                cGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGt                                                                                        >  2:344217/1‑150 (MQ=255)
                                                                cGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGt                                                                                        <  1:1099558/150‑1 (MQ=255)
                                                                 gACATCTTCCTGGCCCACCTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGGGGGCAGCAGGGCGAGGGCAGCGACAGCGTg                                                                                       <  1:428290/150‑1 (MQ=255)
                                                                 gACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGTCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTg                                                                                       <  1:657285/150‑1 (MQ=255)
                                                                 gACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTg                                                                                       <  1:133235/150‑1 (MQ=255)
                                                                 gACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTg                                                                                       <  1:594980/150‑1 (MQ=255)
                                                                 gACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTg                                                                                       >  1:763127/1‑150 (MQ=255)
                                                                   cATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGc                                                                                     <  1:437722/150‑1 (MQ=255)
                                                                     tCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCa                                                                                   >  2:242436/1‑150 (MQ=255)
                                                                           tGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCcgcg                                                                             >  2:828566/1‑150 (MQ=255)
                                                                              cccAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTAt                                                                          >  1:367160/1‑150 (MQ=255)
                                                                                   cTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCGCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGc                                                                                          <  1:388121/129‑1 (MQ=255)
                                                                                   cTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGc                                                                                          >  2:388121/1‑129 (MQ=255)
                                                                                   cTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGc                                                                     <  1:1313395/150‑1 (MQ=255)
                                                                                     gcgcgcCCACGGCTATGCCACCCTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGt                                                                   <  1:30993/150‑1 (MQ=255)
                                                                                       gcgcCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTgc                                                                 <  1:763615/150‑1 (MQ=255)
                                                                                           ccAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCaa                                                             >  2:443608/1‑150 (MQ=255)
                                                                                                 cTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCCTCCAGCCAAGGCCGGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAAcgccgc                                                       <  2:397842/150‑1 (MQ=255)
                                                                                                  tATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACgccgcc                                                      >  1:195101/1‑150 (MQ=255)
                                                                                                    tGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAg                                                    >  2:74396/1‑150 (MQ=255)
                                                                                                         ccGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGAACAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAAc                                               >  1:1356056/1‑150 (MQ=255)
                                                                                                          cGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACa                                              <  2:809229/150‑1 (MQ=255)
                                                                                                           gTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAg                                             <  2:195101/150‑1 (MQ=255)
                                                                                                           gTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAg                                             <  2:367160/150‑1 (MQ=255)
                                                                                                                ccTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACt                                        >  1:1220271/1‑150 (MQ=255)
                                                                                                                 cTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGCTCTAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTg                                       <  1:14855/150‑1 (MQ=255)
                                                                                                                  tGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTgg                                      >  2:1279664/1‑150 (MQ=255)
                                                                                                                   gCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTggg                                     <  2:926747/150‑1 (MQ=255)
                                                                                                                     ccAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCa                                   >  1:1309421/1‑150 (MQ=255)
                                                                                                                       aaCGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATc                                 <  2:1309421/150‑1 (MQ=255)
                                                                                                                         cGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATCCt                               <  2:1185280/150‑1 (MQ=255)
                                                                                                                             aGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATCCTCaaa                           >  2:1003399/1‑150 (MQ=255)
                                                                                                                              gCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATCCTCAAAc                          >  1:612193/1‑150 (MQ=255)
                                                                                                                                    aaaTCGTGCCCGACCAGGCCGCGCGCCTGGAGCGGGTCCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATCCTCAAACCGCTGa                    <  1:74396/150‑1 (MQ=255)
                                                                                                                                    aaaTCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATCCTCAAACCGCTGa                    <  1:344217/150‑1 (MQ=255)
                                                                                                                                     aaTCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATCCTCAAACCGCTGAt                   >  1:364198/1‑150 (MQ=255)
                                                                                                                                         gTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATCCTCAAACCGCTGATCGAc               >  1:1158284/1‑150 (MQ=255)
                                                                                                                                             ggctccAGGCGCGCGGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATCCTCAAACCGCTGATCGACCcgc           >  1:701586/5‑150 (MQ=255)
                                                                                                                                               gACCAGGCCGCGCGCCGGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATCCTCAAACCGCTGATCGACCCGCGa         <  1:477059/150‑1 (MQ=255)
                                                                                                                                                tccAGGCGCGCGGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATCCTCAAACCGCTGATCGACCCGCGAg        >  2:1126676/2‑150 (MQ=255)
                                                                                                                                                tccAGGCGCGCGGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATCCTCAAACCGCTGATCGACCCGCGAg        >  1:980192/2‑150 (MQ=255)
                                                                                                                                                   aGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATCCTCAAACCGCTGATCGACCCGCGAGTTg     <  2:1158284/150‑1 (MQ=255)
                                                                                                                                                   aGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGAGCAGCAGGGGGATGGCAGCTACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATCCTCAAACCGCTGATCGACCCGCGAGTTg     <  1:806416/150‑1 (MQ=255)
                                                                                                                                                    ggCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATCCTCAAACCGCTGATCGACCCGCGAGTTgg    <  1:1042685/150‑1 (MQ=255)
                                                                                                                                                      ccGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATCCTCAAACCGCTGATCGACCCGCGAGTTgggg  <  1:443608/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
CCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGGCTATGCCACCGTCGACCTGCCCAACGGCAGCGTGAAAATCGTGCCCGACCAGGCCGCGCGCCTGGAGCCGGTGCCAGTGGAGGCGGGCGGGCAGCAGGGCGAGGGCAGCGACAGCGTGGCCACCCGCGTATTCAGCGTGCGCAACGCCGCCAGCGAACAGGTACTGGGCATCCTCAAACCGCTGATCGACCCGCGAGTTGGGG  >  NC_002947/1196125‑1196423

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: