Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,539,407 G→A 5.5% A32T (GCC→ACC)  PP_3127 → exopolysaccharide transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,539,4070GA 5.5% 222.8 / 2.5 73A32T (GCC→ACC) PP_3127exopolysaccharide transport protein
Reads supporting (aligned to +/- strand):  ref base G (37/32);  new base A (2/2);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

TTTGACGTTGACGTACAGGTACCCCGACATCTGAACTTCAGGCAGGTGGCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACCAACAACACCAGCTTTGCCGATGAGCGCATCGAGGTGATCCGTCAA  >  NC_002947/3539261‑3539553
                                                                                                                                                  |                                                                                                                                                  
tttGCCCTTGACGTACAGGTACCCCGACATCTTAACTTCAGGCAGGTGGCCCTATGAAGTCTTACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCTAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTcgcc                                                                                                                                                  <  1:77398/149‑1 (MQ=255)
tttGACGTTGACGTACAGGTACCCCGACATCTGAACTTCAGGCAGGTGGCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTcgcc                                                                                                                                                  >  3:231548/1‑149 (MQ=255)
  tgacgttgacgtACAGGTACCCCGACATCTGAACTTCAGGCAGGTGGCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTcgccgc                                                                                                                                                >  7:181531/1‑149 (MQ=255)
              aCAGGTACCCCGACATCTGAACTTCAGGCAGGTGGCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTcg                                                                                                                                                    <  4:188601/133‑1 (MQ=255)
              aCAGGTACCCCGACATCTGAACTTCAGGCAGGTGGCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTcg                                                                                                                                                    >  3:188601/1‑133 (MQ=255)
                 ggTACCCCGACATCTGAACTTCAGGCAGGTGGCCCTATGAAGTCTGACTACTAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTcgc                                                                                                                                 <  2:156732/149‑1 (MQ=255)
                  gTACCCCGACATCTGAACTTCAGGCAGGTGGCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTcgcg                                                                                                                                >  7:42344/1‑149 (MQ=255)
                        cGACATCTGAACTTCAGGCAGGTGGCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATg                                                                                                                          >  8:37262/1‑149 (MQ=255)
                          aCATCTGAACTTCAGGCAGGTGGCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGt                                                                                                                        <  4:318431/149‑1 (MQ=255)
                           cATCTGAACTTCAGGCAGGTGGCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGt                                                                                                                        <  4:231548/148‑1 (MQ=255)
                               tGAACTTCAGGCAGGTGGCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTgccgc                                                                                                                   >  6:168115/1‑149 (MQ=255)
                                 aaCTTCAGGCAGGTGGCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGccgccg                                                                                                                 <  5:168115/149‑1 (MQ=255)
                                  aCTTCAGGCAGGTGGCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGa                                                                                                                <  7:37262/149‑1 (MQ=255)
                                    ttCAGGCAGGTGGCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATc                                                                                                              >  2:56631/1‑149 (MQ=255)
                                       aggcaggTGGCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTAc                                                                                                           <  8:42344/149‑1 (MQ=255)
                                       aggcaggTGGCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGAAAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTAc                                                                                                           >  4:75957/1‑149 (MQ=255)
                                             gTGGCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTc                                                                                                      >  6:276097/1‑148 (MQ=255)
                                               ggCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGAc                                                                                                   >  5:146112/1‑149 (MQ=255)
                                                 cccTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCg                                                                                                 >  5:314080/1‑149 (MQ=255)
                                                   cTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCgg                                                                                                >  1:74557/1‑148 (MQ=255)
                                                          gTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTc                                                                                        >  6:225401/1‑149 (MQ=255)
                                                          gTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCACCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTc                                                                                        >  6:214432/1‑149 (MQ=255)
                                                          gTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCACAGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGTCACCATTc                                                                                        >  6:209225/1‑149 (MQ=255)
                                                               aCTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGccgccg                                                                                                                 <  8:57320/119‑1 (MQ=255)
                                                               aCTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCtggtg                                                                                   <  5:225401/149‑1 (MQ=255)
                                                               aCTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCtggtg                                                                                   <  5:276097/149‑1 (MQ=255)
                                                                   cGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGt                                                                               >  6:65024/1‑149 (MQ=255)
                                                                   cGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGt                                                                               >  2:176701/1‑149 (MQ=255)
                                                                    gAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTc                                                                              <  6:146112/149‑1 (MQ=255)
                                                                        tctcCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCgcag                                                                          >  7:82132/1‑149 (MQ=255)
                                                                          tcCCTCATAGATTACATCGCCTTTATCAAGGACAGGGCCGTGGTGTTAGGGGTGAGTGCCTCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCtggtg                                                                                   <  8:49515/138‑1 (MQ=255)
                                                                          tcCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTAtggtg                                                                                   >  7:49515/1‑138 (MQ=255)
                                                                          tcCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGACGCCGATCTACCAGTCGACCGGCACCATTATGGTGGACTCgcagca                                                                        >  1:327644/1‑149 (MQ=255)
                                                                            ccTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCgcagca                                                                        <  8:181531/147‑1 (MQ=255)
                                                                            ccTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGa                                                                      <  3:75957/149‑1 (MQ=255)
                                                                              tCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGAtt                                                                    >  3:27768/1‑149 (MQ=255)
                                                                              tCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGAtt                                                                    >  1:143814/1‑149 (MQ=255)
                                                                                agagATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATg                                                                                                                          >  6:8984/1‑93 (MQ=255)
                                                                                 gagaTTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCg                                                                 >  5:157345/1‑149 (MQ=255)
                                                                                  agaTTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGc                                                                <  1:176701/149‑1 (MQ=255)
                                                                                        cATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAg                                                           <  1:56631/148‑1 (MQ=255)
                                                                                             ccATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTgc                                                                                                                      >  4:79056/1‑84 (MQ=255)
                                                                                             ccATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTgc                                                                                                                      <  3:79056/84‑1 (MQ=255)
                                                                                             ccATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGgtg                                                     >  5:245910/1‑149 (MQ=255)
                                                                                              cATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCtggt                                                                                    >  8:262444/1‑117 (MQ=255)
                                                                                              cATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGc                                                          <  6:172618/143‑1 (MQ=255)
                                                                                              cATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGc                                                          >  5:172618/1‑143 (MQ=255)
                                                                                              cATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGgtgt                                                    <  4:27768/149‑1 (MQ=255)
                                                                                                ttATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTc                                                   >  7:129936/1‑148 (MQ=255)
                                                                                                  aTCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCACCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGAc                                                <  5:209225/149‑1 (MQ=255)
                                                                                                   tCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCACCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGAc                                                <  5:214432/148‑1 (MQ=255)
                                                                                                      aGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACcaac                                            <  6:245910/149‑1 (MQ=255)
                                                                                                      aGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACcaac                                            <  8:129936/149‑1 (MQ=255)
                                                                                                        gACAGGGCGCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCAGCGAACGTGGGCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACcaacaa                                          <  2:327644/149‑1 (MQ=255)
                                                                                                        gACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACcaacaa                                          >  4:180253/1‑149 (MQ=255)
                                                                                                        gACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACcaacaa                                          <  3:180253/149‑1 (MQ=255)
                                                                                                         aCAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACcaacaa                                          >  5:5322/1‑148 (MQ=255)
                                                                                                         aCAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACCaacaac                                         <  6:5322/149‑1 (MQ=255)
                                                                                                         aCAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGAGTTCGCCTGAGCTGGTGTCGACcaacaa                                          >  7:73467/1‑148 (MQ=255)
                                                                                                          aatggCACTGGTGTTAGGGGTTAGTGCCGCAGGAATTCTCGCCGCGACCGTGTCCGGCGCGTTGATGGTGCCTCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACCaacaac                                         <  8:73467/145‑1 (MQ=255)
                                                                                                             ggCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACCAacaaca                                        <  1:263331/146‑1 (MQ=255)
                                                                                                             ggCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACCAacaaca                                        >  2:263331/1‑146 (MQ=255)
                                                                                                              gCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCg                                                                 >  8:196212/1‑120 (MQ=255)
                                                                                                              gCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCg                                                                 <  7:196212/120‑1 (MQ=255)
                                                                                                                ccTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACCAACAACACCAGCt                                  >  8:223154/1‑149 (MQ=255)
                                                                                                                ccTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACCAACAACACCAGCt                                  >  5:242444/1‑149 (MQ=255)
                                                                                                                    gtgtTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACCAACAACACCAGCTTTGc                              <  8:82132/149‑1 (MQ=255)
                                                                                                                    gtgtTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACCAACAACACCAGCTTTGc                              <  6:157345/149‑1 (MQ=255)
                                                                                                                                  tGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACCAACAACACCAGCTTTGCCGATGAGCGCATc                 <  2:143814/148‑1 (MQ=255)
                                                                                                                                    ccGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCtg                                                                                      <  6:319334/77‑1 (MQ=255)
                                                                                                                                    ccGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCtg                                                                                      >  5:319334/1‑77 (MQ=255)
                                                                                                                                     cGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACCAACAACACCAGCTTTGCCGATGAGCGCATCGAgg             >  6:295139/1‑149 (MQ=255)
                                                                                                                                      gCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACCAACAACACCAGCTTTGCCGATGAGCGCATCGAGGt            <  5:65024/149‑1 (MQ=255)
                                                                                                                                            aTTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACCAACAACACCAGCTTTGCCGATGAGCGCATCGAGGTGATCCg      >  6:116334/1‑149 (MQ=255)
                                                                                                                                            aTTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACCAACAACACCAGCTTTGCCGATGAGCGCATCGAGGTGATCCg      >  5:213057/1‑149 (MQ=255)
                                                                                                                                             ttCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACcaac                                            <  6:7171/110‑1 (MQ=255)
                                                                                                                                             ttCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACcaac                                            >  5:7171/1‑110 (MQ=255)
                                                                                                                                                tcGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCAGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACCAACAACACCAGCTTTGCCGATGAGCGCATAGAGGTGATCCGTCaa  >  2:223302/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTTGACGTTGACGTACAGGTACCCCGACATCTGAACTTCAGGCAGGTGGCCCTATGAAGTCTGACTACGAGCTCTCCCTCAGAGATTACATCGCCATTATCAAGGACAGGGCCCTGGTGTTAGGGGTGAGTGCCGCAGTCATTCTCGCCGCGACCGTGGCCGTCGCGTTGATGGTGCCGCCGATCTACCAGTCGACCGGCACCATTCTGGTGGAGTCGCAGCAGATTTCGCCTGAGCTGGTGTCGACCAACAACACCAGCTTTGCCGATGAGCGCATCGAGGTGATCCGTCAA  >  NC_002947/3539261‑3539553

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