Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,289,588 G→C 5.8% L246L (CTC→CTG) ‡ PP_4663 ← transmembrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,289,5880GC 5.8% 237.7 / 6.4 86L246L (CTC→CTG) ‡PP_4663transmembrane protein
Reads supporting (aligned to +/- strand):  ref base G (42/37);  new base C (3/2);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-01

AGGCCGCGCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGC  >  NC_002947/5289441‑5289735
                                                                                                                                                   |                                                                                                                                                   
aGGCCGCGCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGa                                                                                                                                                    <  1:37526/149‑1 (MQ=255)
aGGCCGCGCTGGACGTGCACTCAGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGa                                                                                                                                                    <  1:284026/149‑1 (MQ=255)
  gCCGCGCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAgg                                                                                                                                                  <  1:374787/149‑1 (MQ=255)
    cgcgCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGggc                                                                                                                                                >  4:258117/1‑149 (MQ=255)
     gcgcTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAggtac                                                                                                                                               >  5:126306/1‑146 (MQ=255)
     gcgcTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGggcg                                                                                                                                               <  6:25063/149‑1 (MQ=255)
      cgcTGGACGTGCACTCGGGCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGGTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGggcgg                                                                                                                                              <  4:286380/149‑1 (MQ=255)
           gACGTGCTCTCGGCCCTTCCGCACATTGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCTCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCCTCGGCGGCCa                                                                                                                                           <  7:350023/147‑1 (MQ=255)
            aCGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGc                                                                                                                                        <  5:106212/149‑1 (MQ=255)
               tGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCCTCGGCGGCCacctgct                                                                                                                                     <  8:272939/149‑7 (MQ=255)
                 cACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGAt                                                                                                                                    <  4:262559/148‑1 (MQ=255)
                  aCTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGGCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATg                                                                                                                                   >  3:295323/1‑148 (MQ=255)
                   cTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGa                                                                                                                                  >  4:238768/1‑148 (MQ=255)
                    tcggcccTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGaa                                                                                                                                 >  7:260328‑M2/8‑148 (MQ=255)
                    tcggcccTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTAACTGGAGCGCGATTTCGAAATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCAGAGGGCGGCCAGCCCGATGaa                                                                                                                                 >  1:88831‑M2/8‑148 (MQ=255)
                      ggcccTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGaaaaa                                                                                                                              <  1:212881‑M2/144‑1 (MQ=255)
                        cccTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAg                                                                                                                             >  4:173968‑M2/4‑148 (MQ=255)
                         ccTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGCGCGCGGTTTCGATATCCCGCGGGTGGACCTGCTCAGCGCTTCCTCCGCGAACCGTTGATCAGCAGGTGGCCGCCGCGGGCGGGCAGGCCGGTGAAAAAGc                                                                                                                            >  1:84075‑M2/3‑148 (MQ=255)
                         ccTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGc                                                                                                                            >  4:283014‑M2/3‑148 (MQ=255)
                          cTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCaa                                                                                                                          >  3:23626‑M2/2‑149 (MQ=255)
                           ttCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGc                                                                                                                                        <  4:116265/134‑1 (MQ=255)
                           ttCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGc                                                                                                                                        <  4:196477/134‑1 (MQ=255)
                           ttCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGc                                                                                                                                        >  3:196477/1‑134 (MQ=255)
                           ttCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGc                                                                                                                                        >  3:116265/1‑134 (MQ=255)
                             ccGCACTTGGCCGCCTTCTTCTAAAAGCGATTGTAATACCTGGAGCGCTTTTTCTATATCCAGGGGGTCTACTTGCTCAGCGCTTCGTCCGAAAACCTTTTATCAGCTGGTGCCCGCCGTGGGCGGCCAGGCCGATGAAAAAGCAACg                                                                                                                        <  4:370852/148‑1 (MQ=255)
                               gCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCtt                                                                                                                     <  3:283014/149‑1 (MQ=255)
                               gCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCtt                                                                                                                     >  4:337350/1‑149 (MQ=255)
                               gCACATGCCCGCCTTCTTCGAAAACCGCATGCACTACCGGGAGCGCGATTGCGATACCCAGCGGGTCGACATGCACCGCGCGCCGCCCGCTAACCTTTTCACACCAGGTGGCCCCCGCGGGCGGCCAGGCCGATGAAAAAGCAACGCtt                                                                                                                     <  4:312781/149‑1 (MQ=255)
                                  caTGGCCGCCTTCTTCTAAACGCGCAGGGACTACCTGGAGCGCGATTTCTATATCCAGCGGGTCGACATGCTCAGCGCTTCGCCCGCAAACCGTTGATCAGCAGGGGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGaa                                                                                                                  <  8:774/149‑1 (MQ=255)
                                    tGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAAca                                                                                                                <  4:37697/149‑1 (MQ=255)
                                      gccgccTTCTTCTAAACGCGCTTGCACTAGCTGGCGCGGTATTTCTATATCCGGCGGGTCGACATGCTCCGCGCTTCGCCCGCTAACCTTTTATCCGCCGGGGGCCGCCGCGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAAcag                                                                                                               <  8:50434/148‑1 (MQ=255)
                                      gccgccTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAAcag                                                                                                               >  4:174727/1‑148 (MQ=255)
                                           cttcttCGAAAAGCGCATGGACTACCTGCAGCGCGATTTCTATATCCACCGGGTCGACATGCTCAGCGCTGCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAgcgc                                                                                                         <  2:84075/149‑1 (MQ=255)
                                            ttcttcGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAgcgcg                                                                                                        >  7:8374/1‑149 (MQ=255)
                                              gttcTAAAAGCCCTTGGACTACGTGTAGCGCTTTTTCTATATCCAGGGGGTCGACATGCTCAGCGCTTCGTCCGCTAACGTTTTATAAGCAGGTGGCCGCCTTGGGGGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAgcgcgc                                                                                                       <  4:175018/147‑1 (MQ=255)
                                                  gAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATc                                                                                                  >  5:10121/1‑149 (MQ=255)
                                                  gAAAAGCGCATGGACAACCTGGAGCGCGATTTCGAAATCCAGCGGGTCGACATGCTCACCGCTTCGTACGCTAACCGTTGAACAGCAGGTGGCCGCCGAGGCCGGCCAGGCCGAAGAAAAACAAAAGCATCAACAGCAGCGCCCTGATc                                                                                                  >  7:319083/1‑149 (MQ=255)
                                                            tgcaCTACCTGGAGCGCGATTTCTATATCACGCGGGTCGACATGCGCAGCGCTTCGTCCGCTAACCGTTGATCAGCCGGTGGCCGCCTAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCa                                                                                        <  7:49639/146‑1 (MQ=255)
                                                                cTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAgg                                                                                                                                                  <  1:47537/87‑1 (MQ=255)
                                                                cTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAgg                                                                                                                                                  >  2:47537/1‑87 (MQ=255)
                                                                cTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCa                                                                                    >  3:89110/1‑149 (MQ=255)
                                                                   ccaggAGCGCGATTTCGATATCCAGCGGGACCACATCCACAGCGCATCGTCAGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACTCAAGAACAGCAGCGCGCTGATACGCTCACGCAGCCACTg                                                                                 >  4:361754/4‑149 (MQ=255)
                                                                   ccTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTg                                                                                 >  6:325933/1‑149 (MQ=255)
                                                                    cTGGAGCGCTATTTCTATATCCTGCGGGTCGACATGCTCCGCGCTTGGTCCGCTAACCGTTGAGCAGCTGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTgg                                                                                <  4:295323/149‑1 (MQ=255)
                                                                     tGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGc                                                                               >  2:134720/1‑149 (MQ=255)
                                                                       gAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGc                                                                               <  5:325933/147‑1 (MQ=255)
                                                                         gcgcGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCAGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGc                                                                           >  5:156896/1‑149 (MQ=255)
                                                                           gcgATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCa                                                                         <  7:355052/149‑1 (MQ=255)
                                                                            cgATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAg                                                                        >  8:340877/1‑149 (MQ=255)
                                                                               tttCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCAt                                                                     >  4:275034/1‑149 (MQ=255)
                                                                                ttCGATATCCAGCGGGTCGAAATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCAt                                                                     <  3:213898/148‑1 (MQ=255)
                                                                                   gATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTg                                                                                             >  2:71904/1‑121 (MQ=255)
                                                                                         cAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCaa                                                            <  8:8374/148‑1 (MQ=255)
                                                                                         cAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACAGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAg                                                           >  5:46270/1‑149 (MQ=255)
                                                                                          aGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAgg                                                          >  4:178151/1‑149 (MQ=255)
                                                                                              ggTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCgg                                                      >  7:332258/1‑149 (MQ=255)
                                                                                               gTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGgc                                                     <  7:1482/149‑1 (MQ=255)
                                                                                                    cATGCTCAGCGCTTCGGCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAg                                                <  8:260328/149‑1 (MQ=255)
                                                                                                             cgcTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATcagcag                                       <  3:238768/149‑1 (MQ=255)
                                                                                                              gcTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGTGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCagcagc                                      <  4:89110/149‑1 (MQ=255)
                                                                                                              gcTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCagcagc                                      <  3:258117/149‑1 (MQ=255)
                                                                                                                ttCGTGCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGa                                    >  7:340055/1‑149 (MQ=255)
                                                                                                                  cGTCCGCTAACCGTTTATCCGCTGGTGGCCGCCGTGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGAc                                   <  4:23626/148‑1 (MQ=255)
                                                                                                                     ccGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGTCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGccc                               >  4:138696/1‑149 (MQ=255)
                                                                                                                      cGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGcccc                              >  2:269670/1‑149 (MQ=255)
                                                                                                                       gCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGcccc                              >  2:339327/1‑148 (MQ=255)
                                                                                                                        cTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCa                            >  2:316065/1‑149 (MQ=255)
                                                                                                                         tAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAg                           >  1:77076/1‑149 (MQ=255)
                                                                                                                             cGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGcc                       >  5:80428/1‑149 (MQ=255)
                                                                                                                              gTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGccccctaggcct                      >  5:94284/1‑142 (MQ=255)
                                                                                                                              aacggttagcagGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAGCGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCt                      >  6:17908/8‑149 (MQ=255)
                                                                                                                               ttGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCa                                                 <  6:276188/121‑1 (MQ=255)
                                                                                                                               ttGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCa                                                 >  5:276188/1‑121 (MQ=255)
                                                                                                                                  aTCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCATCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGc                   <  6:80428/148‑1 (MQ=255)
                                                                                                                                   tCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCa                 <  3:173968/149‑1 (MQ=255)
                                                                                                                                   tCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCCGGCCa                 <  7:340877/149‑1 (MQ=255)
                                                                                                                                     tcggcctGGCCGCCCTCGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCAcc               >  7:74000/7‑149 (MQ=255)
                                                                                                                                       cagGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCa                                                 >  5:708/1‑113 (MQ=255)
                                                                                                                                       cagGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCa                                                 <  6:708/113‑1 (MQ=255)
                                                                                                                                       cagGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCAccc              >  2:115542/1‑148 (MQ=255)
                                                                                                                                        agGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAGCGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCaa            <  5:17908/149‑1 (MQ=255)
                                                                                                                                         gggggCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCTACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAg           <  2:88831/146‑1 (MQ=255)
                                                                                                                                          gTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAg                                                           >  3:334159/1‑100 (MQ=255)
                                                                                                                                          gTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAg                                                           <  4:334159/100‑1 (MQ=255)
                                                                                                                                          ctGGCCGCCCTCGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAgg          >  3:349900/2‑149 (MQ=255)
                                                                                                                                           tGGCCGCCCTCGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGCCGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAgg          >  6:206626/1‑148 (MQ=255)
                                                                                                                                               cgccgAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCCCCCAAGGCATcc     >  7:1648/1‑149 (MQ=255)
                                                                                                                                                 ccgAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTg   >  4:140005/1‑149 (MQ=255)
                                                                                                                                                  cgAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGACCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGc  >  3:41230/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGGCCGCGCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGC  >  NC_002947/5289441‑5289735

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