Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,856,783 C→A 5.2% P237H (CCC→CAC)  PP_5133 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,856,7830CA 5.2% 315.5 / 5.6 96P237H (CCC→CAC) PP_5133membrane protein
Reads supporting (aligned to +/- strand):  ref base C (53/38);  new base A (3/2);  total (56/40)
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.80e-01

CTCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGTGGCGCTG  >  NC_002947/5856638‑5856930
                                                                                                                                                 |                                                                                                                                                   
cTCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCcga                                                                                                                                                  >  5:321280/1‑149 (MQ=255)
 tCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCcgac                                                                                                                                                 >  3:191895/1‑149 (MQ=255)
            ggCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGAt                                                                                                                                      >  5:58930/1‑149 (MQ=255)
              cATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATcc                                                                                                                                    >  7:184803/1‑149 (MQ=255)
                   tACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCgcgc                                                                                                                                >  3:275903/1‑148 (MQ=255)
                   tACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCgcgc                                                                                                                                >  6:282781/1‑148 (MQ=255)
                     cTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAg                                                                                                                             >  6:209353/1‑149 (MQ=255)
                      tGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgc                                                                                                                            >  6:110023/1‑149 (MQ=255)
                        gTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAg                                                                                                                                             <  6:36865/130‑1 (MQ=255)
                        gTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAg                                                                                                                                             >  5:36865/1‑130 (MQ=255)
                        gTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcgg                                                                                                                          >  1:173694/1‑149 (MQ=255)
                         tCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcggc                                                                                                                         <  1:286965/149‑1 (MQ=255)
                           ttGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGa                                                                                                                       >  6:198795/1‑149 (MQ=255)
                           ttGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGa                                                                                                                       >  6:285586/1‑149 (MQ=255)
                           ttGCTGGCGACCCGGCACTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGa                                                                                                                       >  1:136064/1‑149 (MQ=255)
                                   gACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGc                                                                                                               <  5:285586/149‑1 (MQ=255)
                                      ccGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAg                                                                                                            >  1:176706/1‑149 (MQ=255)
                                      ccGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAg                                                                                                            >  8:175963/1‑149 (MQ=255)
                                      ccGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAg                                                                                                            >  5:91489/1‑149 (MQ=255)
                                         gcgcTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCCCCAGCGGCGACAGCCTGCCAGccc                                                                                                         <  5:110023/149‑1 (MQ=255)
                                           gcTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTg                                                                                                       >  3:145838/1‑149 (MQ=255)
                                           gcTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTg                                                                                                       <  8:184803/149‑1 (MQ=255)
                                             tATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGa                                                                                                     >  5:202547/1‑149 (MQ=255)
                                              aTGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCcgc                                                                                                                                  <  3:145478/119‑1 (MQ=255)
                                              aTGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCcgc                                                                                                                                  >  4:145478/1‑119 (MQ=255)
                                              aTGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGAc                                                                                                    >  8:261152/1‑149 (MQ=255)
                                                       tCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGc                                                                                           >  4:14541/1‑149 (MQ=255)
                                                         tGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTgg                                                                                                      <  6:167221/136‑1 (MQ=255)
                                                         tGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTgg                                                                                                      >  5:167221/1‑136 (MQ=255)
                                                         tGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGc                                                                                           >  7:243194/1‑147 (MQ=255)
                                                           ggAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGcc                                                                                          <  4:159945/146‑1 (MQ=255)
                                                           ggAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAgg                                                                                       >  4:144690/1‑149 (MQ=255)
                                                            gAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGc                                                                                      >  1:118177/1‑149 (MQ=255)
                                                             aaaCCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGcc                                                                                                          <  3:79169/128‑1 (MQ=255)
                                                             aaaCCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGcc                                                                                                          >  4:79169/1‑128 (MQ=255)
                                                               accaccTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTg                                                                                   <  5:198795/149‑1 (MQ=255)
                                                               accaccTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTg                                                                                   <  3:144690/149‑1 (MQ=255)
                                                                ccaccTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCTAGCACTGGGCGCCCTCCCGTCGTTGCTGTGCTTCGCCTTGCCCGACGAAGCAATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTgg                                                                                  <  8:209939/149‑1 (MQ=255)
                                                                ccaccTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCGCGCCAGGCCTgg                                                                                  <  7:245563/149‑1 (MQ=255)
                                                                ccaccTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTgg                                                                                  <  4:275903/149‑1 (MQ=255)
                                                                  accTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGc                                                                                                                   >  5:184885/1‑114 (MQ=255)
                                                                  accTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGc                                                                                                                   <  6:184885/114‑1 (MQ=255)
                                                                   ccTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGc                                                                                >  5:37166/1‑148 (MQ=255)
                                                                    cTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCa                                                                              <  6:91489/149‑1 (MQ=255)
                                                                     ttGCTCGCCGCCGACCCGTTCATTAACCTGACCCAAGCACTGGGCTCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAg                                                                             >  7:266291/1‑149 (MQ=255)
                                                                            ccgccgACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAgctggctg                                                                      >  3:54739/1‑149 (MQ=255)
                                                                                cgACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCg                                                                  >  1:31511/1‑149 (MQ=255)
                                                                                 gACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCgg                                                                 >  3:86371/1‑149 (MQ=255)
                                                                                    ccGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgt                                                              >  2:65700/1‑149 (MQ=255)
                                                                                     cGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtg                                                             <  7:175963/149‑1 (MQ=255)
                                                                                      gTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtgg                                                            >  6:271427/1‑149 (MQ=255)
                                                                                      gTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtgg                                                            >  1:85954/1‑149 (MQ=255)
                                                                                        tcattcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGCCGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtggtg                                                          >  2:306274/1‑149 (MQ=255)
                                                                                          attcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCt                                                        >  2:140926/1‑149 (MQ=255)
                                                                                             caCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTgg                                                      <  2:176706/148‑1 (MQ=255)
                                                                                             caCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTgg                                                      <  5:271427/148‑1 (MQ=255)
                                                                                              aCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcggc                                                                                                                         <  4:306976/80‑1 (MQ=255)
                                                                                              aCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcggc                                                                                                                         >  3:306976/1‑80 (MQ=255)
                                                                                              aCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGt                                                     <  4:145838/148‑1 (MQ=255)
                                                                                               gcTGACCCAAGCCCTCGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCGCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGtat                                                   <  8:243194/148‑1 (MQ=255)
                                                                                                cTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGtata                                                  >  3:170716/1‑149 (MQ=255)
                                                                                                     ccAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTgg                                                      >  4:18828/1‑140 (MQ=255)
                                                                                                     ccAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTgg                                                      <  3:18828/140‑1 (MQ=255)
                                                                                                     ccAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGcc                                             >  1:148353/1‑149 (MQ=255)
                                                                                                     ccAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGcc                                             >  6:4033/1‑149 (MQ=255)
                                                                                                            cTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGcc                                      >  5:310879/1‑149 (MQ=255)
                                                                                                             tGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGTCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGccc                                     >  7:26138/1‑149 (MQ=255)
                                                                                                               ggCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGcccc                                    <  7:235479/148‑1 (MQ=255)
                                                                                                               ggCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCg                                   >  1:105406/1‑149 (MQ=255)
                                                                                                               ggCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCg                                   >  1:101544/1‑149 (MQ=255)
                                                                                                                gcgcCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGc                                  <  2:148353/149‑1 (MQ=255)
                                                                                                                gcgcCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGc                                  <  2:31511/149‑1 (MQ=255)
                                                                                                                  gcCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtg                                                             >  4:178839/1‑120 (MQ=255)
                                                                                                                  gcCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtg                                                             <  3:178839/120‑1 (MQ=255)
                                                                                                                        ccGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCctgctgct                          >  3:23552/1‑149 (MQ=255)
                                                                                                                         cgtcgtTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCctgctgctg                         >  7:139880/1‑149 (MQ=255)
                                                                                                                             gtTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGcc                      <  2:101544/148‑1 (MQ=255)
                                                                                                                             gtTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGcc                      <  2:105406/148‑1 (MQ=255)
                                                                                                                              ttGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCgg                    <  4:86371/149‑1 (MQ=255)
                                                                                                                              ttGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCgg                    <  2:118177/149‑1 (MQ=255)
                                                                                                                              ttGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCgg                    <  4:170716/149‑1 (MQ=255)
                                                                                                                                gCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGt                  >  2:91859/1‑149 (MQ=255)
                                                                                                                                  tGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGCCCTGCTGCCCCGCCTGCTGCTGGCCggg                   <  2:66171/146‑1 (MQ=255)
                                                                                                                                  tGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGCCCTGCTGCCCCGCCTGCTGCTGGCCggg                   >  1:66171/1‑146 (MQ=255)
                                                                                                                                   tggCTTCGTCGGGCACGGCGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCt                                                                                    <  1:230847/79‑1 (MQ=38)
                                                                                                                                   tggCTTCGTCGGGCACGGCGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCt                                                                                    >  2:230847/2‑80 (MQ=38)
                                                                                                                                    ggCTTCGTCGGGCACGGCGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATg              >  4:32157/1‑149 (MQ=255)
                                                                                                                                    ggCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTa                <  7:261152/147‑1 (MQ=255)
                                                                                                                                      cTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGc             <  7:176999/148‑1 (MQ=255)
                                                                                                                                      cTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGc             >  8:176999/1‑148 (MQ=255)
                                                                                                                                       ttCGTCGGGCACGGCGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGc             >  6:315940/1‑147 (MQ=255)
                                                                                                                                       ttCGTCGGGCACGGCGAAGCCATGATCCGCGCCACCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGc             <  5:315940/147‑1 (MQ=255)
                                                                                                                                          gCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCtgt         <  1:91859/148‑1 (MQ=255)
                                                                                                                                          gCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCtgt         <  2:173694/148‑1 (MQ=255)
                                                                                                                                            cGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGtggc      <  2:85954/149‑1 (MQ=255)
                                                                                                                                               gCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGtggcgct   >  1:129514/1‑149 (MQ=255)
                                                                                                                                                cccGCCGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTTGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCtgtg        <  2:45991/143‑1 (MQ=255)
                                                                                                                                                cccGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGtggcgcca  <  3:14541/149‑3 (MQ=255)
                                                                                                                                                cccGACGAAGCCATGATCCGCGCCAGCGGCGACAACCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGCTTATGCCtgtg        >  1:45991/1‑143 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGTGGCGCTG  >  NC_002947/5856638‑5856930

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