Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,205,609 A→T 5.9% L275Q (CTG→CAG) ‡ PP_4586 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,205,6090AT 5.9% 254.2 / 4.3 84L275Q (CTG→CAG) ‡PP_4586hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (42/37);  new base T (3/2);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

CTGCTTGCTCAGACGCGCCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTGGGCGGGT  >  NC_002947/5205467‑5205753
                                                                                                                                              |                                                                                                                                                
cTGCTTGCTCAGACGCGCCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCa                                                                                                                                            >  1:15405/1‑149 (MQ=255)
    ttgctCAGACGCGCCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCaa                                                                                                                                        >  8:87909/1‑149 (MQ=255)
        tCAGACGCGCCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTcg                                                                                                                                    >  3:269818/1‑149 (MQ=255)
          aGACGCGCCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTcgc                                                                                                                                   <  2:126699/148‑1 (MQ=255)
            aCGCGCCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCTAACAGTTGCACCAATTCGCgca                                                                                                                                <  6:372234/149‑1 (MQ=255)
             cgcgCCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCgcag                                                                                                                               <  6:244281/149‑1 (MQ=255)
              gcgcCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCgcag                                                                                                                               >  4:324071/1‑148 (MQ=255)
               cgcCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAATTGGTGCAACTGTTCgcgc                                                                                                                                 >  4:218142/1‑145 (MQ=255)
               cgcCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAATTGGTGCAACTGTTCgcgc                                                                                                                                 <  3:218142/145‑1 (MQ=255)
                       aGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTa                                                                                                                     <  3:125443/149‑1 (MQ=255)
                       aGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGACGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTa                                                                                                                     >  5:374680/1‑149 (MQ=255)
                        gACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTAt                                                                                                                    >  5:79019/1‑149 (MQ=255)
                          ccTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTAtt                                                                                                                   <  5:263431/148‑1 (MQ=255)
                                   ggTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCg                                                                                                         >  1:11006/1‑149 (MQ=255)
                                    gTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGc                                                                                                        >  8:243481/1‑149 (MQ=255)
                                     tcatcaCTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCt                                                                                                       >  7:68483/1‑149 (MQ=255)
                                       atcaCTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGcc                                                                                                                >  6:362231/1‑138 (MQ=255)
                                       atcaCTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGcc                                                                                                                <  5:362231/138‑1 (MQ=255)
                                       atcaCTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGa                                                                                                     >  5:284253/1‑149 (MQ=255)
                                        tcaCTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTGACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGAc                                                                                                    >  5:384206/1‑149 (MQ=255)
                                           cTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCAt                                                                                                 >  1:180617/1‑149 (MQ=255)
                                            tCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACGAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCGATTGCCGTTTTCGCTGACCATg                                                                                                >  7:239514/1‑149 (MQ=255)
                                            tCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCgcag                                                                                                                               >  8:306929/1‑118 (MQ=255)
                                            tCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCgcag                                                                                                                               <  7:306929/118‑1 (MQ=255)
                                            tCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATg                                                                                                >  8:292018/1‑149 (MQ=255)
                                               tCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGc                                                                                             <  3:283496/149‑1 (MQ=255)
                                                cccgcccgAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGcc                                                                                            >  3:228779/1‑149 (MQ=255)
                                                  cgcccgAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCgg                                                                                          <  5:258091/149‑1 (MQ=255)
                                                     ccgAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCa                                                                                       <  6:284253/149‑1 (MQ=255)
                                                          cggcgTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGtt                                                                                  <  7:87909/149‑1 (MQ=255)
                                                           ggcgTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGtt                                                                                  >  1:319348/1‑148 (MQ=255)
                                                            gcgTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCa                                                                                <  1:371495/149‑1 (MQ=255)
                                                             cgTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAg                                                                               <  2:11006/149‑1 (MQ=255)
                                                               tAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGTCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAgg                                                                              >  5:362714/1‑148 (MQ=255)
                                                                     cTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAAc                                                                        >  5:355950/1‑148 (MQ=255)
                                                                      tCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAAc                                                                        >  4:318487/1‑147 (MQ=255)
                                                                      tCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAAc                                                                        >  2:236800/1‑147 (MQ=255)
                                                                      tCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGa                                                                      >  6:8126/1‑149 (MQ=255)
                                                                        gTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAg                                                                     <  1:124927/148‑1 (MQ=255)
                                                                           ggCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGccc                                                                  <  4:269818/148‑1 (MQ=255)
                                                                            gCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCt                                                                 <  6:374680/148‑1 (MQ=255)
                                                                             cGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTgc                                                               <  6:384206/149‑1 (MQ=255)
                                                                              gTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTgcg                                                              >  7:221427/1‑149 (MQ=255)
                                                                               taataaTCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTgcg                                                              >  7:144303/1‑148 (MQ=255)
                                                                                aataatCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTATTCAGGTACAACGAGCCCTgcg                                                              >  8:59256/1‑147 (MQ=255)
                                                                                 ataatCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACGAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCa                                                            <  7:243481/148‑1 (MQ=255)
                                                                                 ataatCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCa                                                            <  3:318487/148‑1 (MQ=255)
                                                                                    atCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACACCGCGCCCTGCGCAGCCa                                                        >  6:969/1‑149 (MQ=255)
                                                                                    atCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCa                                                        >  4:352159/1‑149 (MQ=255)
                                                                                       ggTCAGGCTCTCTACCAGTCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATa                                                      <  6:362714/148‑1 (MQ=255)
                                                                                       ggTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATa                                                      <  5:337344/148‑1 (MQ=255)
                                                                                       ggTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATa                                                      <  2:15405/148‑1 (MQ=255)
                                                                                       ggTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAg                                                     <  6:285560/149‑1 (MQ=255)
                                                                                           aGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGccc                                                  <  6:252970/148‑1 (MQ=255)
                                                                                             gCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAg                                                     >  2:375199/1‑143 (MQ=255)
                                                                                             gCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAg                                                     <  1:375199/143‑1 (MQ=255)
                                                                                             gCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTc                                               <  6:79019/149‑1 (MQ=255)
                                                                                               tctctACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGc                                             >  2:145974/1‑149 (MQ=255)
                                                                                                   tACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCcaca                                         >  4:62384/1‑149 (MQ=255)
                                                                                                     ccagccagTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATg                                       <  8:68483/149‑1 (MQ=255)
                                                                                                     ccagccagTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCGATTGCCGTTGTCGCGGCCCAGGGGCCGGGCACGGTTCCGGGACAACGCGCCCTGCGCAGCCATCGCCCCTCGcccccgg                                       >  8:327342/1‑144 (MQ=255)
                                                                                                        gccagTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCg                                    <  1:145974/149‑1 (MQ=255)
                                                                                                         ccagTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGt                                   >  7:282411/1‑149 (MQ=255)
                                                                                                               aCTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTAATGCCGTTTTCGCTGACCATGGGCCGGGGACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCACTCGCCACATGCCGTCAcgcg                             >  8:4638/1‑149 (MQ=255)
                                                                                                                cTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTCTTCAGGTACAACGAGCCATCCGCAGCCATAGCCCCTCGCCACATGCCGTCAcgctc                            >  8:160948/1‑147 (MQ=255)
                                                                                                                  gCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACAGGCCGTCCCGCGccc                          >  6:260931/1‑149 (MQ=255)
                                                                                                                     atCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGa                       >  3:990/2‑149 (MQ=255)
                                                                                                                        aCGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGAcc                     <  2:180617/148‑1 (MQ=255)
                                                                                                                           ccccTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCaa                   >  3:148378/1‑147 (MQ=255)
                                                                                                                              cTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAagag                >  8:355434/1‑147 (MQ=255)
                                                                                                                              cTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGagca              <  7:355434/149‑1 (MQ=255)
                                                                                                                              cTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCGATTGCCGTTTTCTCTGACCATGGGCCGGGCACTGTTCGGGTACAACGAGCCCTGCGCAGCAATAGCCCCTCGCCACAAGCCGTCACGCGCCCCGACCAagagac              >  4:63731/1‑147 (MQ=255)
                                                                                                                                  aTCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCtt          <  3:62384/149‑1 (MQ=255)
                                                                                                                                   tCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTg         >  1:40938/1‑149 (MQ=255)
                                                                                                                                     tGCGCGAATTGGTGCAACTGTTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTggg       >  1:198424/1‑149 (MQ=255)
                                                                                                                                      gcgcGAATTGGTGCAACTGTTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCg                                    <  2:370775/119‑1 (MQ=255)
                                                                                                                                      gcgcGAATTGGTGCAACTGTTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCg                                    >  1:370775/1‑119 (MQ=255)
                                                                                                                                       cgcgAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACATCGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTgg        <  5:260931/146‑1 (MQ=255)
                                                                                                                                       cgcgAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTg         <  7:160948/145‑1 (MQ=255)
                                                                                                                                        gcgAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTATTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTGGGCgg    <  7:59256/149‑1 (MQ=255)
                                                                                                                                         cgAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTGGGCggg   <  5:88394/149‑1 (MQ=255)
                                                                                                                                          gAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCcaca                                         <  2:219265/110‑1 (MQ=255)
                                                                                                                                          gAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCcaca                                         >  1:219265/1‑110 (MQ=255)
                                                                                                                                          gAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTGGGCgtgt  <  2:40938/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CTGCTTGCTCAGACGCGCCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTGGGCGGGT  >  NC_002947/5205467‑5205753

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: