Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,358,599 T→G 5.7% I453L (ATC→CTC)  PP_3834 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,358,5990TG 5.7% 223.8 / 4.3 70I453L (ATC→CTC) PP_3834hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base G (2/2);  total (36/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 = 6.06e-01

GATATCGGCGATCAGGTATTCGTCATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCGGGGGCGCTGATT  >  NC_002947/4358451‑4358747
                                                                                                                                                    |                                                                                                                                                     
gATATCGGCGATCAGGTATTCGTCATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGAt                                                                                                                                                       <  4:234875/149‑1 (MQ=255)
            cAGGTATTCGTCATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTg                                                                                                                                           >  6:370883/1‑149 (MQ=255)
              ggTATTCGTCATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGc                                                                                                                                          >  7:216130/1‑148 (MQ=255)
               gTATTCGTCATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGc                                                                                                                                          >  1:80779/1‑147 (MQ=255)
                 aTTCGTCATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAgg                                                                                                                                      >  3:37756/1‑149 (MQ=255)
                    cGTCATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCg                                                                                                                                   >  2:110388/1‑149 (MQ=255)
                     gTCATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCgg                                                                                                                                  <  3:181218/149‑1 (MQ=255)
                       cATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCgg                                                                                                                                  <  3:179103/147‑1 (MQ=255)
                           gACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCTGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCGTGCGCAAGGTAATGCCTAGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATg                                                                                                                            >  6:244991/1‑149 (MQ=255)
                             cAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGt                                                                                                                          >  6:279289/1‑149 (MQ=255)
                             cAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGt                                                                                                                          <  5:145128/149‑1 (MQ=255)
                                tCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGAc                                                                                                                       >  2:191346/1‑149 (MQ=255)
                                tCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCACAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGAc                                                                                                                       >  3:21757/1‑149 (MQ=255)
                                    gATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTc                                                                                                                   <  5:35406/149‑1 (MQ=255)
                                        gtgtTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACg                                                                                                                >  2:248869/1‑148 (MQ=255)
                                         tgtTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCa                                                                                                              >  2:325990/1‑149 (MQ=255)
                                           ttCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCAcc                                                                                                            >  5:280956/1‑149 (MQ=255)
                                           ttCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCAc                                                                                                             <  5:370883/148‑1 (MQ=255)
                                            tCAGCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGGCGGCGTTCCCTAAAGCACACGGCAGGAACCTTCTAAGGCGTTGCGATCGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCAcc                                                                                                            <  1:110388/148‑1 (MQ=255)
                                              aTCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCtt                                                                                                          <  7:211693/148‑1 (MQ=255)
                                                   cTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGAc                                                                                                    <  6:139139/149‑1 (MQ=255)
                                                       tccgtcACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGa                                                                                                 >  4:327924/1‑148 (MQ=255)
                                                        ccgtcACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCaa                                                                                                                                        >  6:313525/1‑108 (MQ=255)
                                                        ccgtcACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCaa                                                                                                                                        <  5:313525/108‑1 (MQ=255)
                                                         cgtcACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCgg                                                                                                                                  <  2:169369/113‑1 (MQ=255)
                                                          gtcACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAg                                                                                              >  8:338976/1‑148 (MQ=255)
                                                             aCCCGCTCTCCTGCGTAGTAGTCGTTCCCGAAAGCACACCGCAGGCACATTCGAACGCCTTGCGCAATGTAATGCCTTTAGTCACGATCAGCATCTCTTTAAAGGCAGGGTCATGGTTACTGTCACGCACCTTTATTACCGACATCATg                                                                                          <  5:76085/149‑1 (MQ=255)
                                                               ccGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATcagca                                                                                                                                                  <  5:226065/91‑1 (MQ=255)
                                                               ccGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATcagca                                                                                                                                                  >  6:226065/1‑91 (MQ=255)
                                                                cGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATgcgc                                                                                       <  1:325990/149‑1 (MQ=255)
                                                                  ctctCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCg                                                                                     <  8:216130/149‑1 (MQ=255)
                                                                       cTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGtt                                                                                <  1:248869/149‑1 (MQ=255)
                                                                           ggcgtcggcgtTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGa                                                                            <  6:280956/149‑1 (MQ=255)
                                                                              gtcggcgtTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTAGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTc                                                                         <  5:244991/149‑1 (MQ=255)
                                                                                  gcgtTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGt                                                                     >  8:244870/1‑149 (MQ=255)
                                                                                   cgtTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGtt                                                                    <  5:279289/149‑1 (MQ=255)
                                                                                        ccGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTc                                                                         <  2:52381/139‑1 (MQ=255)
                                                                                        ccGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTc                                                                         >  1:52381/1‑139 (MQ=255)
                                                                                        ccGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGAc                                                                >  6:43034/1‑148 (MQ=255)
                                                                                          gAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACa                                                               >  3:28061/1‑147 (MQ=255)
                                                                                          gAAAGCACACCGCAGACACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACaaa                                                             <  1:11913/149‑1 (MQ=255)
                                                                                                acacCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCAt                                                       >  3:204668/1‑149 (MQ=255)
                                                                                                      cAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCa                                                  >  5:311328/1‑148 (MQ=255)
                                                                                                       aGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATa                                                <  3:327924/149‑1 (MQ=255)
                                                                                                         gCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCAt                                                       >  7:339565/1‑140 (MQ=255)
                                                                                                         gCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCAt                                                       <  8:339565/140‑1 (MQ=255)
                                                                                                           aCCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGt                                            >  6:193651/1‑149 (MQ=255)
                                                                                                             cTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGc                                                                              <  6:355878/113‑1 (MQ=255)
                                                                                                             cTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGc                                                                              >  5:355878/1‑113 (MQ=255)
                                                                                                              ttCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAATGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGt                                         >  5:347624/1‑149 (MQ=255)
                                                                                                                    cGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑cgacga                                  <  5:2663/149‑1 (MQ=255)
                                                                                                                     gCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑cgacga                                  >  5:24764/1‑148 (MQ=255)
                                                                                                                       cTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATAc                               >  3:68414/1‑149 (MQ=255)
                                                                                                                        ttGCGGAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTAACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCTCATCAGCATAATGTTGT‑cgacgcttct                              >  3:191391/1‑144 (MQ=255)
                                                                                                                          gcgcAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAACGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACTATACTc                             <  6:347624/148‑1 (MQ=255)
                                                                                                                           cgcAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGa                           <  4:68414/149‑1 (MQ=255)
                                                                                                                                     tGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGgc                 >  3:330120/1‑149 (MQ=255)
                                                                                                                                        cTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATTCGCCGTAGTTGCGAGTCAAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGagcg              >  2:174/1‑149 (MQ=255)
                                                                                                                                        cTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGcgcg              >  3:68311/1‑149 (MQ=255)
                                                                                                                                        cTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGcgcg              >  4:297340/1‑149 (MQ=255)
                                                                                                                                           gCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGTCCGACGATACTCGACAGGTCCGGCGCgg             <  1:191346/148‑1 (MQ=255)
                                                                                                                                           gCGGCACGATCACCAGCTCGTGCAAGGCCGGGTCATGGGGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCggg            <  4:37756/148‑1 (MQ=255)
                                                                                                                                            ttgCACGAGCTGCTGATCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCgggg           >  7:59157/3‑148 (MQ=255)
                                                                                                                                            cGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCgggg           <  5:43034/148‑1 (MQ=255)
                                                                                                                                            cGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCgggg           <  6:193778/148‑1 (MQ=255)
                                                                                                                                             ggCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCGGGGgc         <  2:80779/149‑1 (MQ=255)
                                                                                                                                              gCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCTGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCGGGGgcg        >  5:189434/1‑149 (MQ=255)
                                                                                                                                              gCACGAGCTGCTGATCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑cga                                     >  6:286927/1‑120 (MQ=255)
                                                                                                                                              gCACGAGCTGCTGATCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑cga                                     <  5:286927/120‑1 (MQ=255)
                                                                                                                                              gCACGAGCTGCTGATCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCgggg           <  8:59157/146‑1 (MQ=255)
                                                                                                                                                 cGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCTGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCGGGGGCGCTg     <  6:189434/149‑1 (MQ=255)
                                                                                                                                                    tCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCGGGGGCGCTGAtt  >  4:351755/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATATCGGCGATCAGGTATTCGTCATCGACAATCAGGATGGTGTTCATCGGCTCGTCCGTCACCCGCTCTCCTGCGGCGTCGGCGTTCCCGAAAGCACACCGCAGGCACCTTCGAACGCCTTGCGCAAGGTAATGCCTTGCGGCACGATCAGCAGCTCGTGCAAGGCCGGGTCATGGTGACTGTCACGCACCTTGAGGACCGACAGCATGCGCCGTAGTTGCGAGTCCAGTTCGACAAAGCGCATCAGCATAAGGTTGT‑CGACGATACTCGACAGGTCCGGCGCGGGGGCGCTGATT  >  NC_002947/4358451‑4358747

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: ‑