Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,973,302 Δ1 bp 6.1% coding (615/1713 nt) PP_2600 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,973,3020A. 6.1% 416.8 / 19.5 131coding (615/1713 nt)PP_2600hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (60/62);  new base . (4/4);  total (64/67)
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.43e-01

Only 100 of 131 total aligned reads displayed.

GTCGCGACCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTC  >  NC_002947/2973156‑2973449
                                                                                                                                                  |                                                                                                                                                    
gTCGCGACCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCa                                                                                                                                                   >  2:1038607/1‑150 (MQ=255)
    cgACCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGc                                                                                                                                              <  2:672876/150‑1 (MQ=255)
     gACCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCg                                                                                                                                             >  2:445663/1‑150 (MQ=255)
        ccACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCGGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGc                                                                                                                                          >  1:2198018/1‑150 (MQ=255)
        ccACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGCCATGGTCGGGc                                                                                                                                          <  2:187339/150‑1 (MQ=255)
         cACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCATGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCCGGTGAAAACCCCGACCATGGCGGGcc                                                                                                                                         >  1:126063/1‑150 (MQ=255)
         cACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGcc                                                                                                                                         <  2:1189106/150‑1 (MQ=255)
         cACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGcc                                                                                                                                         >  1:1889334/1‑150 (MQ=255)
          aCGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGAc                                                                                                                                                     <  2:1374600/138‑1 (MQ=255)
          aCGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGAc                                                                                                                                                     >  1:1374600/1‑138 (MQ=255)
          aCGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCa                                                                                                                                        >  2:1914612/1‑150 (MQ=255)
           cGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAg                                                                                                                                       <  1:1555603/150‑1 (MQ=255)
               ccGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAggtgg                                                                                                                                   <  2:1889334/150‑1 (MQ=255)
               ccGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAg                                                                                                                                       >  1:1112521/1‑146 (MQ=255)
               ccGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAg                                                                                                                                       <  2:1113079/146‑1 (MQ=255)
               ccGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAg                                                                                                                                       <  2:1112521/146‑1 (MQ=255)
               ccGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAg                                                                                                                                       >  1:1113079/1‑146 (MQ=255)
                 gcgcACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCATTTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTa                                                                                                                                 >  2:1899259/1‑150 (MQ=255)
                 gcgcACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTa                                                                                                                                 <  2:2169639/150‑1 (MQ=255)
                   gcACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACa                                                                                                                               <  1:222259/150‑1 (MQ=255)
                     aCCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAgg                                                                                                                             >  2:845438/1‑150 (MQ=255)
                        tGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTgg                                                                                                                          <  1:1409686/150‑1 (MQ=255)
                          gCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAgg                                                                                                                                      >  2:736437/1‑136 (MQ=255)
                            cACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTggcggc                                                                                                                      >  1:323355/1‑150 (MQ=255)
                               ccTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCaa                                                                                                                   >  2:481854/1‑150 (MQ=255)
                                cTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCaaa                                                                                                                  <  2:467802/150‑1 (MQ=255)
                                    aTCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAgcgc                                                                                                              <  2:1302181/150‑1 (MQ=255)
                                    aTCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAgcgc                                                                                                              <  2:2039147/150‑1 (MQ=255)
                                      cagcagGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCtt                                                                                                            <  2:2351480/150‑1 (MQ=255)
                                       agcagGCGTTTGGCCGGCGCCGACGCTTGCAGTTGGACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTg                                                                                                           >  2:789856/1‑150 (MQ=255)
                                        gcagGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCgg                                                                                                                                            >  2:1288235/1‑116 (MQ=255)
                                        gcagGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCgg                                                                                                                                            <  1:1288235/116‑1 (MQ=255)
                                                  ggCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGc                                                                                                                                              >  2:805180/1‑104 (MQ=255)
                                                  ggCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGc                                                                                                                                              <  1:805180/104‑1 (MQ=255)
                                                  ggCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTgcag                                                                                                >  1:708984/1‑150 (MQ=255)
                                                     cGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCaaa                                                                                                                  >  1:231456/1‑129 (MQ=255)
                                                     cGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCaaa                                                                                                                  >  1:231039/1‑129 (MQ=255)
                                                     cGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCaaa                                                                                                                  <  2:231456/129‑1 (MQ=255)
                                                     cGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCaaa                                                                                                                  <  2:231039/129‑1 (MQ=255)
                                                     cGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTCCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGcagcag                                                                                             <  1:548332/150‑1 (MQ=255)
                                                           cGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCacac                                                                                                      >  2:632499/1‑135 (MQ=255)
                                                           cGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCacac                                                                                                      <  1:632499/135‑1 (MQ=255)
                                                               gCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGc                                                                                   <  1:393744/150‑1 (MQ=255)
                                                                 ttGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGa                                                                                 >  2:833205/1‑150 (MQ=255)
                                                                     aGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGctaccag                                                                             <  2:2400289/150‑7 (MQ=255)
                                                                     aGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCgg                                                                             <  1:445663/150‑1 (MQ=255)
                                                                       gTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTggcg                                                                                                                        >  1:386023/1‑105 (MQ=255)
                                                                       gTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTggcg                                                                                                                        <  2:386023/105‑1 (MQ=255)
                                                                         gCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCa                                                                         <  2:2106026/150‑1 (MQ=255)
                                                                               ccAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGa                                                                   <  1:789856/150‑1 (MQ=255)
                                                                               ccAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCAGGTTGa                                                                   >  1:238753/1‑150 (MQ=255)
                                                                                cAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAg                                                                  <  2:1292477/150‑1 (MQ=255)
                                                                                   gTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGCGCATTGCCCGCCCCACGCGGCTCTTGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTAAGGaa                                                               <  2:573562/150‑1 (MQ=255)
                                                                                         tgCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCgg                                                                                                                                            >  2:1663068/1‑67 (MQ=255)
                                                                                         tgCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCgg                                                                                                                                            <  1:1663068/67‑1 (MQ=255)
                                                                                           ccGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGc                                                                                                                                          >  2:1445643/1‑67 (MQ=255)
                                                                                           ccGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGc                                                                                                                                          <  1:1445643/67‑1 (MQ=255)
                                                                                           ccGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCaa                                                       >  1:1390533/1‑150 (MQ=255)
                                                                                           ccGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCaa                                                       >  1:1155594/1‑150 (MQ=255)
                                                                                             gATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTc                                                                          <  2:324910/129‑1 (MQ=255)
                                                                                               tttAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGAcc                                                                                                                                                    <  1:1726077/54‑1 (MQ=255)
                                                                                                   gCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTAAAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGCGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCa                                               >  2:470447/1‑150 (MQ=255)
                                                                                                     cAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGtt                                                                     >  1:1526155/1‑126 (MQ=255)
                                                                                                     cAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGtt                                                                     <  2:1526155/126‑1 (MQ=255)
                                                                                                       ggCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAggc                                                            <  1:1693833/133‑1 (MQ=255)
                                                                                                       ggCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAggc                                                            >  2:1693833/1‑133 (MQ=255)
                                                                                                          gcgcCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCAt                                                                        >  1:1187498/1‑118 (MQ=255)
                                                                                                             ccaccaGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGa                                     >  2:560220/1‑150 (MQ=255)
                                                                                                                 caGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTc                                                                          >  1:2167487/1‑109 (MQ=255)
                                                                                                                  aGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGc                                <  2:198411/150‑1 (MQ=255)
                                                                                                                  aGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGGGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAACGGCGTCACCGGCCATGGTCGGc                                >  2:1829318/1‑150 (MQ=255)
                                                                                                                      ccACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGc                            <  1:2337205/150‑1 (MQ=255)
                                                                                                                         cacCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCa                         >  1:1985397/1‑150 (MQ=255)
                                                                                                                             gTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGa                                                                                 >  1:1376972/1‑90 (MQ=255)
                                                                                                                             gTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGa                                                                                 <  2:1376972/90‑1 (MQ=255)
                                                                                                                             gTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCa                                               <  1:1946054/124‑1 (MQ=255)
                                                                                                                              ttGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATg                                                                       >  2:1222298/1‑99 (MQ=255)
                                                                                                                              ttGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGa                                                                   >  2:796004/1‑103 (MQ=255)
                                                                                                                              ttGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGa                                                                   <  1:796004/103‑1 (MQ=255)
                                                                                                                                gACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGa                  <  2:1985397/150‑1 (MQ=255)
                                                                                                                                gACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGa                  <  1:1899259/150‑1 (MQ=255)
                                                                                                                                   cAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGTCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCACGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGt               >  1:290989/1‑150 (MQ=255)
                                                                                                                                    aGGTACATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGtt              >  2:1296569/1‑150 (MQ=255)
                                                                                                                                      gTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTc            <  1:1047216/150‑1 (MQ=255)
                                                                                                                                        acctGGCCCGCCATGGTCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTcgc          >  2:863851/4‑150 (MQ=255)
                                                                                                                                         cctGGCCCGCCATGGTCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTcgcg         >  1:209142/3‑150 (MQ=255)
                                                                                                                                         cATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGc                                                                               >  2:2291839/1‑80 (MQ=255)
                                                                                                                                          ctGGCCCGCCATGGTCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGa                  <  1:2266659/139‑1 (MQ=255)
                                                                                                                                          ctGGCCCGCCATGGTCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGa                  >  2:2266659/2‑140 (MQ=255)
                                                                                                                                           tGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCa       <  1:833205/150‑1 (MQ=255)
                                                                                                                                           tGGCCCGCCATGGTCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCa       <  1:863851/150‑1 (MQ=255)
                                                                                                                                            ggCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCAt      <  1:2244606/150‑1 (MQ=255)
                                                                                                                                              cccGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCa       >  1:1158256/1‑147 (MQ=255)
                                                                                                                                              cccGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCa       <  2:1158256/147‑1 (MQ=255)
                                                                                                                                              cccGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATgg    >  1:1119716/1‑150 (MQ=255)
                                                                                                                                               ccGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCtt                                                                                                            >  2:1786080/1‑45 (MQ=255)
                                                                                                                                               ccGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAg                                                                                <  1:77401/73‑1 (MQ=255)
                                                                                                                                               ccGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGa                                                                   >  1:1227931/1‑86 (MQ=255)
                                                                                                                                               ccGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGa                                                                   <  2:1227931/86‑1 (MQ=255)
                                                                                                                                                cGACCATGGCGGGCCAGGTGGTACAGGTGGCGTCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCACGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTc  <  2:290989/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTCGCGACCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGGCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTC  >  NC_002947/2973156‑2973449

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