Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,886,602 C→T 7.0% D423D (GAC→GAT rmuC → DNA recombination limiting inversions protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,886,6020CT 7.0% 214.1 / 3.9 71D423D (GAC→GATrmuCDNA recombination limiting inversions protein
Reads supporting (aligned to +/- strand):  ref base C (35/31);  new base T (3/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

TGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGACCTG  >  NC_002947/1886454‑1886749
                                                                                                                                                    |                                                                                                                                                   
tGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGAc                                                                                                                                                     >  3:224777/1‑149 (MQ=255)
      tCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCGTCCAGGTCCTGGATGAAca                                                                                                                                                <  5:274798/148‑2 (MQ=255)
        aGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCGTCCAGGTCCTGGATGAAcagc                                                                                                                                              <  5:154851/148‑4 (MQ=255)
        aGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTggg                                                                                                                                             >  4:233796/1‑149 (MQ=255)
           ccTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAg                                                                                                                                          <  8:186609/149‑1 (MQ=255)
               ccACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGcc                                                                                                                                        <  2:273592/147‑1 (MQ=255)
                     tGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGc                                                                                                                                      <  6:248672/143‑1 (MQ=255)
                     tGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGc                                                                                                                                      >  5:248672/1‑143 (MQ=255)
                      gCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTgcag                                                                                                                                <  4:224777/148‑1 (MQ=255)
                       cTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTgcagc                                                                                                                               <  5:281037/148‑1 (MQ=255)
                              ccTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGAc                                                                                                                       <  5:298110/149‑1 (MQ=255)
                                ttGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCACTTGCAGCAGGTGGAc                                                                                                                       >  7:123928/1‑147 (MQ=255)
                                 tGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGc                                                                                                                                            >  4:75910/1‑125 (MQ=255)
                                 tGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGc                                                                                                                                            <  3:75910/125‑1 (MQ=255)
                                   cGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAgg                                                                                                                   <  4:22629/148‑1 (MQ=255)
                                   cGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGc                                                                                                                  <  7:131963/149‑1 (MQ=255)
                                     ggTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCTCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCtt                                                                                                                <  6:285349/149‑1 (MQ=255)
                                     ggTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCtt                                                                                                                <  7:105210/149‑1 (MQ=255)
                                      gTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTa                                                                                                               >  7:38564/1‑149 (MQ=255)
                                         aTCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTgca                                                                                                                                 <  2:266751/128‑1 (MQ=255)
                                         aTCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTgca                                                                                                                                 >  1:266751/1‑128 (MQ=255)
                                             aCAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGGCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCgg                                                                                                          >  2:260244/1‑147 (MQ=255)
                                                 ccTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAgcgc                                                                                                    >  7:22195/1‑149 (MQ=255)
                                                       gAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGc                                                                                                                                            <  1:224378/103‑1 (MQ=255)
                                                       gAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGc                                                                                                                                            >  2:224378/1‑103 (MQ=255)
                                                          gCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGc                                                                                           >  2:295294/1‑149 (MQ=255)
                                                          gCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGc                                                                                           >  7:142846/1‑149 (MQ=255)
                                                              gAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCtgtg                                                                                       <  8:8458/149‑1 (MQ=255)
                                                                   tCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAg                                                                                  >  7:38719/1‑149 (MQ=255)
                                                                     aaGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAggg                                                                                <  8:22195/149‑1 (MQ=255)
                                                                        ggCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGt                                                                             <  1:295294/149‑1 (MQ=255)
                                                                          cAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAATCTGTGCGAAGGGCGTgg                                                                           >  8:198990/1‑149 (MQ=255)
                                                                           aGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTgg                                                                           >  6:273165/1‑148 (MQ=255)
                                                                                  ccGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTgg                                                                   >  5:152136/1‑149 (MQ=255)
                                                                                   cgcgAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCtgt                                                                                        >  6:71585/1‑127 (MQ=255)
                                                                                   cgcgAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGt                                                                  >  2:66501/1‑149 (MQ=255)
                                                                                        aaTTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGcc                                                             <  7:13999/149‑1 (MQ=255)
                                                                                        aaTTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGcc                                                             <  8:142846/149‑1 (MQ=255)
                                                                                            gCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTa                                                          <  1:261651/148‑1 (MQ=255)
                                                                                              cGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAgcgcg                                                                                                   >  3:115200/1‑105 (MQ=255)
                                                                                              cGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAgcgcg                                                                                                   <  4:115200/105‑1 (MQ=255)
                                                                                               gAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGc                                                        >  7:17966/1‑147 (MQ=255)
                                                                                                     gCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGt                                                                  <  7:263451/131‑1 (MQ=255)
                                                                                                     gCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGt                                                                  >  8:263451/1‑131 (MQ=255)
                                                                                                     gCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGAGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTAGCAACCTGGTTAGACGTAGCGAGCAGCt                                                >  7:131849/1‑149 (MQ=255)
                                                                                                      ccGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCtt                                               <  3:302000/149‑1 (MQ=255)
                                                                                                            ggctCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTTGGCAGCCGCTTGCAGCAGGTGGACAATGCTTATGAGTCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTATCCGTAGCGagcagc                                                 >  7:155877/1‑141 (MQ=255)
                                                                                                                ctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGagc                                                    <  3:286138/134‑1 (MQ=255)
                                                                                                                ctACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGagc                                                    >  4:286138/1‑134 (MQ=255)
                                                                                                                 tACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGa                                                      <  4:135131/131‑1 (MQ=255)
                                                                                                                 tACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGa                                                      >  3:135131/1‑131 (MQ=255)
                                                                                                                 tACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAgctgctc                                      <  8:38719/147‑2 (MQ=255)
                                                                                                                  aCGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGgc                                   >  5:195166/1‑149 (MQ=255)
                                                                                                                     aCAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGcc                                >  2:28237/1‑149 (MQ=255)
                                                                                                                        aGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGt                             <  3:233796/149‑1 (MQ=255)
                                                                                                                          ttCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGc                           <  6:195166/149‑1 (MQ=255)
                                                                                                                            cGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCACCCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGc                           >  1:168605/1‑147 (MQ=255)
                                                                                                                             gtgcgtttCAGCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCATGTGGACAAGGCTTATGCGTCAGCGCGCAACAAGCTGTGCGAAGGGCGTAGCAACCTGTTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAg                        <  8:131849/143‑1 (MQ=255)
                                                                                                                              tGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTa                                                                >  7:104251/1‑108 (MQ=255)
                                                                                                                              tGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTa                                                                <  8:104251/108‑1 (MQ=255)
                                                                                                                                cagTTCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGc                       >  3:204165/3‑147 (MQ=255)
                                                                                                                                 agTTCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGc                       >  7:160951/2‑146 (MQ=255)
                                                                                                                                 agTTCGTCCAGGTCCTGGATGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGc                       >  6:112319/2‑146 (MQ=255)
                                                                                                                                   ttCATCCAGGACCTGGACGAACTGGGCAGCCTCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAaga                   <  8:123928/148‑1 (MQ=255)
                                                                                                                                     cATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCt                >  7:202050/1‑149 (MQ=255)
                                                                                                                                        ccAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGc             >  7:60472/1‑149 (MQ=255)
                                                                                                                                         cAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCAGCTGGGCGCCCGTGCCAGCAAGAGCTTGcc            <  1:28237/149‑1 (MQ=255)
                                                                                                                                               cTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGa      >  2:230436/1‑149 (MQ=255)
                                                                                                                                               cTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGa      >  5:254223/1‑149 (MQ=255)
                                                                                                                                               cTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGa      >  2:178398/1‑149 (MQ=255)
                                                                                                                                                  gACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGAcat   >  7:273283/1‑147 (MQ=255)
                                                                                                                                                   aCGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGACctg  >  3:87317/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGATTGTCAGCCCTACCACGTTGCTGGCCACCTTGCGGGTTATCGACAGCCTGTGGAAGCAGGAGCGTCAAGGGCAGAACGCCCGCGAAATTGCCGAGCGGGCCGGCTGGCTCTACGACAAGTTCGTGCTGTTCATCCAGGACCTGGACGAACTGGGCAGCCGCTTGCAGCAGGTGGACAAGGCTTATGCGGCAGCGCGCAACAAGCTGTGCGAAGGGCGTGGCAACCTGGTTAGCCGTAGCGAGCAGCTTAAGCTGCTGGGCGCCCGTGCCAGCAAGAGCTTGCCGGCAGACCTG  >  NC_002947/1886454‑1886749

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