Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,848,334 C→G 7.5% A33A (GCC→GCG PP_3397 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,848,3340CG 7.5% 214.6 / 6.9 67A33A (GCC→GCGPP_3397hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/30);  new base G (2/3);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

GTCACCAGCGCCAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCACCCGGGG  >  NC_002947/3848186‑3848482
                                                                                                                                                    |                                                                                                                                                      
gTCACCAGCGCCAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGcgcc                                                                                                                                                        >  4:95801/1‑149 (MQ=255)
     cAGCGCCAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGtt                                                                                                                                                  <  6:255622/149‑1 (MQ=255)
       gcgcCAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTcg                                                                                                                                                <  4:39553/149‑1 (MQ=255)
         gcCAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTcgcc                                                                                                                                              <  2:265806/149‑1 (MQ=255)
               ggtggtAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTg                                                                                                                                        <  7:55522/149‑1 (MQ=255)
                         aCTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTc                                                                                                                                      <  6:322661/141‑1 (MQ=255)
                         aCTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTc                                                                                                                                      >  5:322661/1‑141 (MQ=255)
                              gAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCAGTGCGCTGTTGGCCGCCGGTTTCTGTCTTACCGTTACCAGCGCCGCGTACGCCGCGCTGTCTCATGTCGATCCg                                                                                                                         >  1:34803/1‑149 (MQ=255)
                               aaGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCgg                                                                                                                        <  5:126451/149‑1 (MQ=255)
                                         gctggtgcCATGAACAGACAAAGACTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGt                                                                                                                                >  8:54162/1‑131 (MQ=255)
                                         gctggtgcCATGAACAGACAAAGACTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGt                                                                                                                                <  7:54162/131‑1 (MQ=255)
                                               gcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTggtaa                                                                                                                                    >  5:265656/1‑117 (MQ=255)
                                               gcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTggtaa                                                                                                                                    <  6:265656/121‑5 (MQ=255)
                                                ccATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGCTCCGGGGCCCTACACCTTCGc                                                                                                        >  6:187695/1‑148 (MQ=255)
                                                  aTGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTgg                                                                                                                                       >  7:111502/1‑114 (MQ=255)
                                                  aTGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTgg                                                                                                                                       <  8:111502/115‑2 (MQ=255)
                                                  aTGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCca                                                                                                      >  5:184123/1‑148 (MQ=255)
                                                    gAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTg                                                                                                                                    >  1:175388/1‑116 (MQ=255)
                                                    gAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTg                                                                                                                                    <  2:175388/116‑1 (MQ=255)
                                                        aGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTgg                                                                                               <  5:41431/149‑1 (MQ=255)
                                                          aCAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTCGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCCCAGACTGATGTCGATCCGGGGCCCTACACCTTCGACACAGGTGGc                                                                                              >  4:81491/1‑148 (MQ=255)
                                                           cAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTa                                                                                            >  7:145174/1‑149 (MQ=255)
                                                            aaaGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTgg                                                                                               <  3:261062/145‑1 (MQ=255)
                                                            aaaGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTgg                                                                                               >  4:261062/1‑145 (MQ=255)
                                                            aaaGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTAt                                                                                           >  8:90228/1‑149 (MQ=255)
                                                              acgCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCTGCGTTCGCCGCGCTGTCTGATGTCGATc                                                                                                                           <  1:18848/114‑1 (MQ=255)
                                                              aGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATc                                                                                                                           >  2:18848/1‑115 (MQ=255)
                                                                   cAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGAtgt                                                                                    >  7:6345/1‑149 (MQ=255)
                                                                    aGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGAtgt                                                                                    <  7:90228/148‑1 (MQ=255)
                                                                        ttCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTa                                                                                >  1:106530/1‑148 (MQ=255)
                                                                           aaGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGACGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCAGGGGCCCTACACCTTCGCCACAGGTTGGTATCCGATGTGGTACAAg                                                                            >  8:297154/1‑149 (MQ=255)
                                                                                   agcaTGAACACTTGGACGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCCGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGcc                                                                     >  8:96479/1‑148 (MQ=255)
                                                                                    gcaTGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTc                                                                   <  8:145174/149‑1 (MQ=255)
                                                                                     caTGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATg                                                                                                                                 <  8:75136/86‑1 (MQ=255)
                                                                                     caTGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATg                                                                                                                                 >  7:75136/1‑86 (MQ=255)
                                                                                      aTGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGa                                                                 >  5:289117/1‑149 (MQ=255)
                                                                                        gAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATc                                                                                          <  4:28643/122‑1 (MQ=255)
                                                                                        gAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATc                                                                                          >  3:28643/1‑122 (MQ=255)
                                                                                        gAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGAc                                                                >  7:150455/1‑148 (MQ=255)
                                                                                                  tCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGa                                                    >  2:60943/1‑149 (MQ=255)
                                                                                                       ccGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCCGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCtgt                                               <  7:96479/149‑1 (MQ=255)
                                                                                                         gTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGc                                             <  3:81491/149‑1 (MQ=255)
                                                                                                             gcTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGACTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGt                                         >  2:214707/1‑149 (MQ=255)
                                                                                                                   ggCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGt                                    <  8:150455/148‑1 (MQ=255)
                                                                                                                   ggCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGt                                    <  8:6345/148‑1 (MQ=255)
                                                                                                                   ggCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTg                                   >  2:180399/1‑149 (MQ=255)
                                                                                                                    gccgccGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGACTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTg                                   <  1:214707/148‑1 (MQ=255)
                                                                                                                     ccgccgGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGc                                  >  8:169469/1‑148 (MQ=255)
                                                                                                                           gTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCtgt                                               <  8:43448/129‑1 (MQ=255)
                                                                                                                           gTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCtgt                                               >  7:43448/1‑129 (MQ=255)
                                                                                                                              tctctcTTACCGTTACCAGCGCCGCGTTCGCCGCGCTTTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCAc                         <  7:297154/148‑1 (MQ=255)
                                                                                                                                tctctTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCg                      >  4:253686/1‑149 (MQ=255)
                                                                                                                                  tctTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTg                    >  8:108412/1‑149 (MQ=255)
                                                                                                                                    ttACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCTATCCGGGGCCCTACACCTTCGCca                                                                                                      <  1:48248/66‑1 (MQ=255)
                                                                                                                                    ttACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCca                                                                                                      >  2:48248/1‑66 (MQ=255)
                                                                                                                                    ttACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGcca                                <  5:292631/135‑1 (MQ=255)
                                                                                                                                    ttACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTg                    <  7:169469/147‑1 (MQ=255)
                                                                                                                                    ttACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTAGACCAATCCCTGGAGCTGTGCCAGTCGCGTGcca                                >  6:292631/1‑135 (MQ=255)
                                                                                                                                      aCCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTg                                                                                                <  6:47686/70‑1 (MQ=255)
                                                                                                                                      aCCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTg                                                                                                >  5:47686/1‑70 (MQ=255)
                                                                                                                                      aCCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTgc                   <  3:95801/146‑1 (MQ=255)
                                                                                                                                      aCCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTgc                   <  2:106530/146‑1 (MQ=255)
                                                                                                                                          ttACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTCGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAgg             <  1:60943/149‑1 (MQ=255)
                                                                                                                                          ttACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGgc            <  3:253686/149‑1 (MQ=255)
                                                                                                                                                gcgcCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCTGGCGCAcgg      <  6:289117/149‑3 (MQ=255)
                                                                                                                                                gcgcCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCAccc      >  5:149264/1‑149 (MQ=255)
                                                                                                                                                gcgcCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCAccc      >  2:50073/1‑149 (MQ=255)
                                                                                                                                                gcgcCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCAccc      >  2:263824/1‑149 (MQ=255)
                                                                                                                                                   ccgcGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCACCCggg   <  1:180399/149‑1 (MQ=255)
                                                                                                                                                    cgcgTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCACCCgggg  >  2:15657/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GTCACCAGCGCCAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCACCCGGGG  >  NC_002947/3848186‑3848482

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: