Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,742,960 C→T 6.1% H5Y (CAC→TAC) ‡ PP_2397 → EF hand domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,742,9600CT 6.1% 289.6 / 6.5 98H5Y (CAC→TAC) ‡PP_2397EF hand domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (44/48);  new base T (2/4);  total (46/52)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

CCTGCCAGCGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAACCCCAGTGGAGCT  >  NC_002947/2742813‑2743107
                                                                                                                                                   |                                                                                                                                                    
cctgccAGCGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCAc                                                                                                                                                    >  2:312131/1‑150 (MQ=255)
  tgccAGCGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCt                                                                                                                                                  >  2:1128937/1‑150 (MQ=255)
    ccAGCGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTg                                                                                                                                                >  2:919878/1‑150 (MQ=255)
    ccAGCGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTg                                                                                                                                                >  1:788313/1‑150 (MQ=255)
    ccAGCGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTg                                                                                                                                                >  2:600298/1‑150 (MQ=255)
     cAGCGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTgg                                                                                                                                               >  1:2054439/1‑150 (MQ=255)
      aGCGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTggg                                                                                                                                              <  1:1941026/150‑1 (MQ=255)
      aGCGAAAACGCCTTGGCGGCACCCGAGCCCATAGCGGCGGCCCCGCAAGCCAACGAAACGTTAAGCAACCCGGCCAAGCACACCAAGCCATCGGCCCCGCCCAACCCGCTGTAAGACTTATCACGCCTCATGGCCCCAAGGCACCTTggg                                                                                                                                              <  1:165205/150‑1 (MQ=255)
       gggAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTgggg                                                                                                                                             <  1:2215716/148‑1 (MQ=255)
        cGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGc                                                                                                                                            <  2:1490932/150‑1 (MQ=255)
        cGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGc                                                                                                                                            <  2:124617/150‑1 (MQ=255)
        cGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAACCCCGCCAAGCCAGCGGCCCCGCCGAATCCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGTGCCTTGGGGc                                                                                                                                            <  2:5844/150‑1 (MQ=255)
          aaaaGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGTGCCTTGGGGcca                                                                                                                                          <  1:212038/150‑2 (MQ=255)
             aGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGtt                                                                                                                                       <  2:1749437/150‑1 (MQ=255)
              ggCGGTGGCGGCACCCGAGCCCGTAGCGGCGGAACAACAAGCATATGATAATTTGAGCAAGCCGGCCAAGGCCCTCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGccatgag                                                                                                                                      <  1:2271364/150‑6 (MQ=255)
                 ggTTGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGAAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGAGTACGCCTCAGGGCCCCAAGGCACCTTGGGGCCCGTTGCCg                                                                                                                                   <  1:2159344/150‑1 (MQ=255)
                  gTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGTGCCTTGGGGcca                                                                                                                                          >  2:1089209/1‑141 (MQ=255)
                  gTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGTGCCTTGGGGcca                                                                                                                                          <  1:1089209/142‑2 (MQ=255)
                         cACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGgtgt                                                                                                                           <  2:484011/150‑1 (MQ=255)
                         cACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGgtgt                                                                                                                           <  1:380130/150‑1 (MQ=255)
                           cccGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGgtgtgt                                                                                                                         >  1:526605/1‑150 (MQ=255)
                            ccGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTa                                                                                                                        <  1:855282/150‑1 (MQ=255)
                            ccGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTa                                                                                                                        >  1:148269/1‑150 (MQ=255)
                                 cccGTAGCGGCGTCCCCACAAGACAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCACCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAACGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCa                                                                                                                   <  1:548711/150‑1 (MQ=255)
                                  ccGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCAc                                                                                                                  >  2:789494/1‑150 (MQ=255)
                                  ccGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAACGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCAc                                                                                                                  >  1:1950929/1‑150 (MQ=255)
                                    gTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTg                                                                                                                >  2:1337139/1‑150 (MQ=255)
                                         ggcggcCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCt                                                                                                           <  1:786851/150‑1 (MQ=255)
                                          gcggcCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGt                                                                                                               <  2:469267/145‑1 (MQ=255)
                                          gcggcCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGt                                                                                                               >  1:469267/1‑145 (MQ=255)
                                                ccGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCtgtg                                                                                                   >  1:285389/1‑150 (MQ=255)
                                                  gCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGtgct                                                                                                 <  1:505348/150‑1 (MQ=255)
                                                    aaGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGGGGTGTGGATCGCACTGTATTTGCCTGtgctgc                                                                                               >  1:386749/1‑150 (MQ=255)
                                                     aGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCa                                                                                              <  1:600298/150‑1 (MQ=255)
                                                     aGCCAAGGAAAAGAACAGCACGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCa                                                                                              <  1:1970574/150‑1 (MQ=255)
                                                       ccAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGa                                                                                            <  1:919878/150‑1 (MQ=255)
                                                        cAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGAc                                                                                           <  2:2128877/150‑1 (MQ=255)
                                                         aaGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGAcc                                                                                          >  1:1622323/1‑150 (MQ=255)
                                                          aGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGgtgt                                                                                                                           <  2:2166472/117‑1 (MQ=255)
                                                          aGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGgtgt                                                                                                                           >  1:2166472/1‑117 (MQ=255)
                                                            gAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTAc                                                                                       >  2:1942482/1‑150 (MQ=255)
                                                              aaaGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTg                                                                                     <  1:1787057/150‑1 (MQ=255)
                                                               aaGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGc                                                                                    <  2:285389/150‑1 (MQ=255)
                                                                 ggACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCa                                                                                  <  1:1128937/150‑1 (MQ=255)
                                                                       cAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCa                                                                                                                                                     <  2:1056655/78‑1 (MQ=255)
                                                                       cAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCa                                                                                                                                                     >  1:1056655/1‑78 (MQ=255)
                                                                       cAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACaaa                                                                            >  1:2183689/1‑150 (MQ=255)
                                                                         aGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAg                                                                          >  1:1735786/1‑150 (MQ=255)
                                                                         aGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAg                                                                          <  2:135165/150‑1 (MQ=255)
                                                                            cgtccAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCtgtg                                                                                                   <  2:1663101/119‑1 (MQ=255)
                                                                            cGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCtgtg                                                                                                   >  1:1663101/1‑122 (MQ=255)
                                                                              gCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAgcgcgc                                                                     <  1:1337139/150‑1 (MQ=255)
                                                                               ccAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGggtggt                                                                                                                             >  2:452958/1‑94 (MQ=255)
                                                                               ccAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGggtggt                                                                                                                             <  1:452958/94‑1 (MQ=255)
                                                                                 aaGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTg                                                                  <  2:517775/150‑1 (MQ=255)
                                                                                 aaGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTg                                                                  >  1:1189064/1‑150 (MQ=255)
                                                                                 aaGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTg                                                                  >  2:1937687/1‑150 (MQ=255)
                                                                                     ccGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCAc                                                              >  1:1831234/1‑150 (MQ=255)
                                                                                     ccGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACACCCAACCGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCAc                                                              >  1:1591293/1‑150 (MQ=255)
                                                                                       gCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGa                                                            >  2:1566712/1‑150 (MQ=255)
                                                                                         cAAGCCAGCGGCCCCGCCCAAACAGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTTTGTATGGCAGTGTATCTGCATGTGCTGGAGACCTACTGCCATAGAAAAGAGGGCCTGGGACGATa                                                          >  1:848462/1‑150 (MQ=255)
                                                                                              cAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCACGCCTGGCACGACACATcc                                                     >  2:523334/1‑150 (MQ=255)
                                                                                               aGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTTCCATACAAAAGCGCGCCTGGCACGATACATccc                                                    <  2:1622323/150‑1 (MQ=255)
                                                                                                gCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATTGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATcccc                                                   <  2:256653/150‑1 (MQ=255)
                                                                                                gCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATcccc                                                   <  2:788313/150‑1 (MQ=255)
                                                                                                  gggcccGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGa                                                 <  1:1731812/147‑1 (MQ=255)
                                                                                                  ggCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGa                                                 <  2:2054439/150‑1 (MQ=255)
                                                                                                     cccGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAgcgc                                                                       >  2:2196656/1‑125 (MQ=255)
                                                                                                     cccGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAgcgc                                                                       <  1:2196656/125‑1 (MQ=255)
                                                                                                       cGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAAc                                            <  1:1937687/150‑1 (MQ=255)
                                                                                                         cccAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCg                                          <  2:2183689/150‑1 (MQ=255)
                                                                                                                  cTGGAAGACTGATCACGCCTCATCGCCCCAAGGCACCTTGGGGCACGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCt                                 <  1:312131/150‑1 (MQ=255)
                                                                                                                     gAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTaag                              >  1:1026638/1‑149 (MQ=255)
                                                                                                                     gAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAAt                              >  2:497732/1‑150 (MQ=255)
                                                                                                                     gAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAAt                              >  2:1577286/1‑150 (MQ=255)
                                                                                                                     gAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAAt                              >  1:1125636/1‑150 (MQ=255)
                                                                                                                       aGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATgg                            <  1:723185/150‑1 (MQ=255)
                                                                                                                       aGACTGATAACGTGTACTTGCCCCTAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATgg                            <  2:1591293/150‑1 (MQ=255)
                                                                                                                            gATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTccc                       >  2:2228588/1‑149 (MQ=255)
                                                                                                                              tCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACt                     <  2:526605/150‑1 (MQ=255)
                                                                                                                                  gCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTcc                 >  1:307368/1‑150 (MQ=255)
                                                                                                                                    cTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCaa               >  2:1403515/1‑150 (MQ=255)
                                                                                                                                     tCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAAc              >  1:1577640/1‑150 (MQ=255)
                                                                                                                                       aTGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACACACCGCAACGGCCTAATGGCTCCACTTTCCAAccc            >  1:693045/1‑150 (MQ=255)
                                                                                                                                       aTGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAAccc            <  1:1566712/150‑1 (MQ=255)
                                                                                                                                       aTGGCCCCAAGGCACCTTGGGGCCCGTTGCCCGGTGGCATGTGGAGCAGGGTAGTGCTCTGTGCTGCAGAAAACCGGCTATACAAAAGCGCGCCTGGCACTATACAGCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAAccc            <  1:523334/150‑1 (MQ=255)
                                                                                                                                        tGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAAcccc           >  1:2285830/1‑150 (MQ=255)
                                                                                                                                        gggCCCCAAGGTGCCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGtgctgc                                                                                               <  2:1801881/65‑1 (MQ=38)
                                                                                                                                        gggCCCCAAGGTGCCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGtgctgc                                                                                               >  1:1801881/2‑66 (MQ=38)
                                                                                                                                         ggCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAACCCCa          <  2:1831234/150‑1 (MQ=255)
                                                                                                                                         ggCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAACCCCa          <  1:1577286/150‑1 (MQ=255)
                                                                                                                                            cccAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTcc                        <  1:952799/133‑1 (MQ=255)
                                                                                                                                            cccAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTcc                        >  2:952799/1‑133 (MQ=255)
                                                                                                                                              cAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATAc                                                         >  1:1906021/1‑98 (MQ=255)
                                                                                                                                              cAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATAc                                                         <  2:1906021/98‑1 (MQ=255)
                                                                                                                                                 ggCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAAccc            <  1:697106/140‑1 (MQ=255)
                                                                                                                                                 ggCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAAccc            >  2:697106/1‑140 (MQ=255)
                                                                                                                                                 ggCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAACCCCAGTGGAGCt  >  1:1239632/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCTGCCAGCGAAAAGGCGGTGGCGGCACCCGAGCCCGTAGCGGCGGCCCCGCAAGCCAAGGAAAAGGACAGCAAGCCGGCCAAGCCCGCCAAGCCAGCGGCCCCGCCCAAACCGCTGGAAGACTGATCACGCCTCATGGCCCCAAGGCACCTTGGGGCCCGTTGCCGGGTGGTGTGTATCGCACTGTATCTGCCTGTGCTGCAGACCTACTGCCATACAAAAGCGCGCCTGGCACGATACATCCCCGACAAACCGCAACGGCCTAATGGCTCCACTTTCCAACCCCAGTGGAGCT  >  NC_002947/2742813‑2743107

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