Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,848,338 T→A 7.3% F35I (TTC→ATC) ‡ PP_3397 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,848,3380TA 7.3% 201.1 / 5.1 68F35I (TTC→ATC) ‡PP_3397hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/30);  new base A (2/3);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CAGCGCCAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCACCCGGGG  >  NC_002947/3848191‑3848482
                                                                                                                                                    |                                                                                                                                                 
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                                                                                                                                      >  5:322661/1‑141 (MQ=255)
                    aCTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTc                                                                                                                                      <  6:322661/141‑1 (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                                                                                                                                       <  8:111502/115‑2 (MQ=255)
                                             aTGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTgg                                                                                                                                       >  7:111502/1‑114 (MQ=255)
                                             aTGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCca                                                                                                      >  5:184123/1‑148 (MQ=255)
                                               gAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTg                                                                                                                                    <  2:175388/116‑1 (MQ=255)
                                               gAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTg                                                                                                                                    >  1:175388/1‑116 (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)
                                                                                                                                                    |                                                                                                                                                 
CAGCGCCAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCACCCGGGG  >  NC_002947/3848191‑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: