Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,883,572 C→G 7.5% A197A (GCC→GCG) ‡ PP_3428 → TPR repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,883,5720CG 7.5% 215.3 / 7.3 67A197A (GCC→GCG) ‡PP_3428TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base G (2/3);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

CCATCGGCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACTCGGCCTTTCTGCAG  >  NC_002947/3883427‑3883719
                                                                                                                                                  |                                                                                                                                                   
ccATCGGCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTg                                                                                                                                                  <  7:264433/149‑1 (MQ=255)
      gCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCACGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACCAGCACGCCGCTGCAGGGCTGCagcgg                                                                                                                                            >  3:269264/1‑148 (MQ=255)
      gCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACAGGCACGCCGCTGCAGGGCTGCagcgg                                                                                                                                            <  4:269264/149‑2 (MQ=255)
      gCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCagcga                                                                                                                                            <  4:84067/149‑1 (MQ=255)
      gCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCagcga                                                                                                                                            <  3:144974/149‑1 (MQ=255)
       ccTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCagcgag                                                                                                                                           <  8:268697/149‑1 (MQ=255)
       ccTCTAATTCGGCCCCCTGCAGTCAGCCCTGTAGCAATACCAGCCAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCTAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGGGCTGCagcggc                                                                                                                                           <  7:162804/149‑3 (MQ=255)
        ctctATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCagcgagc                                                                                                                                          >  7:8889/1‑149 (MQ=255)
        ctctATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCagcgagc                                                                                                                                          <  7:352016/149‑1 (MQ=255)
         tctATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCagcgagcg                                                                                                                                         >  4:26539/1‑149 (MQ=255)
            aTTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTc                                                                                                                                      >  3:86360/1‑149 (MQ=255)
               tCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTcgcg                                                                                                                                   >  1:20809/1‑149 (MQ=255)
                  gCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGtt                                                                                                                                <  6:144592/149‑1 (MQ=255)
                   cccccTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTg                                                                                                                               >  4:135426/1‑149 (MQ=255)
                     cccTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGcc                                                                                                                             >  3:110291/1‑149 (MQ=255)
                        tGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGTCTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGccccc                                                                                                                          >  8:94177/1‑149 (MQ=255)
                         gCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCg                                                                                                                         <  6:32404/149‑1 (MQ=255)
                          cAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCgg                                                                                                                        <  8:161146/149‑1 (MQ=255)
                           aGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGc                                                                                                                       >  4:97573/1‑149 (MQ=255)
                            gTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGTCTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCa                                                                                                                      <  7:94177/149‑1 (MQ=255)
                               aGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACG‑CCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCAGTTGCCCCCGGCAACCt                                                                                                                  >  2:274191/1‑149 (MQ=255)
                                   cTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGa                                                                                                               <  8:8889/149‑1 (MQ=255)
                                       aGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGcc                                                                                                           <  3:97573/149‑1 (MQ=255)
                                        gCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGccc                                                                                                          <  2:20809/149‑1 (MQ=255)
                                          aaTACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTg                                                                                                        <  6:149235/149‑1 (MQ=255)
                                            tACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGc                                                                                                      >  4:306859/1‑149 (MQ=255)
                                               cAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTCCAGCCCTGCagc                                                                                                                                              >  1:78636/1‑107 (MQ=255)
                                                   aaGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAg                                                                                               <  6:281733/149‑1 (MQ=255)
                                                     gCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAgg                                                                                              >  3:282909/1‑148 (MQ=255)
                                                      cccTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGccc                                                                                                                            >  3:81465/1‑117 (MQ=255)
                                                      cccTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGccc                                                                                                                            <  4:81465/117‑1 (MQ=255)
                                                            aGCAACAGCCTGAAGATGCCGACG‑CCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCAGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGc                                                                                     <  1:274191/149‑1 (MQ=255)
                                                                 cAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTc                                                                                   >  3:300083/1‑147 (MQ=255)
                                                                  aGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAg                                                                                <  4:86360/149‑1 (MQ=255)
                                                                    ccTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAgcg                                                                              >  6:316832/1‑149 (MQ=255)
                                                                     cTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCGACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAgcgc                                                                             >  5:326041/1‑149 (MQ=255)
                                                                      tGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAgcgc                                                                             >  2:66386/1‑148 (MQ=255)
                                                                       gAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAgcgc                                                                             >  1:190558/1‑147 (MQ=255)
                                                                         aGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAgg                                                                         >  5:163571/1‑149 (MQ=255)
                                                                         aGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAgg                                                                         <  3:135426/149‑1 (MQ=255)
                                                                          gATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAAc                                                                                                                    <  6:91608/105‑1 (MQ=255)
                                                                          gATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAAc                                                                                                                    <  6:166886/105‑1 (MQ=255)
                                                                          gATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAAc                                                                                                                    >  5:166886/1‑105 (MQ=255)
                                                                          gATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAAc                                                                                                                    >  5:91608/1‑105 (MQ=255)
                                                                            tGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGcc                                                                      <  6:163571/149‑1 (MQ=255)
                                                                                gACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGccacca                                                                  <  5:316832/149‑1 (MQ=255)
                                                                                             cACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTg                                                                                      <  5:113922/117‑1 (MQ=255)
                                                                                             cACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTg                                                                                      >  6:113922/1‑117 (MQ=255)
                                                                                             cACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCg                                                      >  3:141942/1‑149 (MQ=255)
                                                                                              aCCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGc                                                     <  4:300083/149‑1 (MQ=255)
                                                                                                cTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCa                                                   <  3:26539/149‑1 (MQ=255)
                                                                                                  ggCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACt                                                                                                  >  2:226689/1‑100 (MQ=255)
                                                                                                        ccTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTgg                                           <  2:190558/149‑1 (MQ=255)
                                                                                                        ccTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTgg                                           <  7:144598/149‑1 (MQ=255)
                                                                                                         cTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGAt                                                                                                              <  4:140082/81‑1 (MQ=255)
                                                                                                         cTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGAt                                                                                                              >  3:140082/1‑81 (MQ=255)
                                                                                                          tGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGcatcgag                                                                 <  1:368323/125‑7 (MQ=255)
                                                                                                                   aCGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTg                                <  4:110291/149‑1 (MQ=255)
                                                                                                                                 acgcacgcCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGa                  >  6:352079/1‑149 (MQ=255)
                                                                                                                                   gcacgcCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACt                >  8:241631/1‑149 (MQ=255)
                                                                                                                                     acgcCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACTCg              <  7:316047/149‑1 (MQ=255)
                                                                                                                                        tcgcTGCAGGGCTGCAGCGAGCGCTCGCGGTTG‑CCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGa                         <  7:61329/133‑1 (MQ=255)
                                                                                                                                        tcgcTGCAGGGCTGCAGCGAGCGCTCGCGGTTG‑CCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGAAGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGa                         >  8:61329/2‑134 (MQ=255)
                                                                                                                                         cgcTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACTCGGCCt          >  4:133378/1‑149 (MQ=255)
                                                                                                                                             gCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACTCGGCCTTTCt      >  2:139352/1‑149 (MQ=255)
                                                                                                                                              cAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACTCGGCCTTTCTg     >  2:272288/1‑149 (MQ=255)
                                                                                                                                               aGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACTCGGCCTTTCTGc    <  3:133378/149‑1 (MQ=255)
                                                                                                                                                 cccTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACTCGGCCTTTCTGCAg  >  7:209400/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCATCGGCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGT‑GCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACTCGGCCTTTCTGCAG  >  NC_002947/3883427‑3883719

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: ‑