Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,885,304 C→T 6.8% P397S (CCA→TCA)  PP_3429 → two component system response regulator/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,885,3040CT 6.8% 172.7 / 2.4 59P397S (CCA→TCA) PP_3429two component system response regulator/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (25/30);  new base T (2/2);  total (27/32)
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.17e-01

CGCCTTCGAACCCTTCTTCAGCACCAAGGCGACCGGCCAAGGCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGCCAGGCGCTGGGCGAAAACGGCTTCCCCTG  >  NC_002947/3885156‑3885449
                                                                                                                                                       |                                                                                                                                                 
cGCCTTCGAACCCTTCTTCAGCACCAAGGCGACCGGCCAAGGCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGGGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGc                                                                                                                                                   >  8:278717/1‑149 (MQ=255)
 gCCTTCGAACCCTTCTTCAGCACCAAGGCGACCGGCCAAGGCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGc                                                                                                                                                   >  6:3368/1‑148 (MQ=255)
      cGAACCCTTCTTCAGCACCAAGGCGACCGGCCAAGGCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGc                                                                                                                                              >  2:274795/1‑148 (MQ=255)
          cccTTCTTCAGCACCAAGGCGACCGGCCAAGGCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCTTGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGt                                                                                                                                         <  5:305894/149‑1 (MQ=255)
            cttcttCAGCACCAAGGCGACCGGCCAAGGCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTc                                                                                                                                        >  3:135007/1‑148 (MQ=255)
              tcttcAGCACCAAGGCGACCGGCCAAGGCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGAt                                                                                                                                     >  7:302661/1‑149 (MQ=255)
                  cAGCACCAAGGCGACCGGCCAAGGCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATc                                                                                                                                    <  7:278717/146‑1 (MQ=255)
                  cAGCACCAAGGCGACCGGCCAAGGCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAg                                                                                                                                  <  7:81042/148‑1 (MQ=255)
                      aCCAAGGCGACCGGCCAAGGCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCt                                                                                                                             <  4:135007/149‑1 (MQ=255)
                      aCCAAGGCGACCGGCCAAGGCATTGGCCTGGGGTTGTCGATGCTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCt                                                                                                                             <  6:305136/149‑1 (MQ=255)
                        cAAGGCGACCGGCCAAGGCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGt                                                                                                                           >  5:248444/1‑149 (MQ=255)
                           ggCGACCGGCCAAGGCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCAAGG‑‑GAACCCAGGTGGAGCAGTATTTGCCACGCCATGTCGATCAGGCCCTGTCa                                                                                                                         >  7:18067/1‑148 (MQ=255)
                            gCGACCGGCCAAGGCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATa                                                                                                                       <  2:306744/149‑1 (MQ=255)
                            gCGACCGGCCAAGGCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATa                                                                                                                       <  8:302661/149‑1 (MQ=255)
                            gCGACCGGCCAAGGCATTGGCCTGGGGGTGGCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATa                                                                                                                       <  7:219941/149‑1 (MQ=255)
                              cccaggCCAATGCCTTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATa                                                                                                                       >  6:259516/5‑147 (MQ=255)
                               ccaggCCAATGCCTTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAg                                                                                                                    >  5:323963/4‑149 (MQ=255)
                                ccGGCCAAGGCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGc                                                                                                                   <  6:233827/149‑1 (MQ=255)
                                      aaGGCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCg                                                                                                             >  5:83004/1‑149 (MQ=255)
                                        ggCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGt                                                                                                           <  3:291563/149‑1 (MQ=255)
                                           aTTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTgc                                                                                                         >  4:1776/1‑148 (MQ=255)
                                            ttGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGc                                                                                                                                <  5:226462/124‑1 (MQ=255)
                                            ttGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGc                                                                                                                                >  6:226462/1‑124 (MQ=255)
                                                ccTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACTCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTgccagcca                                                                                                   <  6:310866/149‑1 (MQ=255)
                                                 cTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGccagccag                                                                                                  <  5:3368/149‑1 (MQ=255)
                                                  tGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCCTGTCGATCAGGCCCGGTCATACAGCGAGCCGGTGCCAGCCAgg                                                                                                 >  4:309331/1‑149 (MQ=255)
                                                         gTCGATGGTGTACGGTTTCAGCAAGCAGTCCCAAGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACg                                                                                          >  2:189291/1‑149 (MQ=255)
                                                            gATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCt                                                                                       <  3:1776/149‑1 (MQ=255)
                                                               ggTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGc                                                                                     <  6:83004/148‑1 (MQ=255)
                                                                  gtACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGc                                                                                  <  1:129223/148‑1 (MQ=255)
                                                                   tACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGcc                                                                                 >  4:20420/1‑148 (MQ=255)
                                                                                       ccATTGGGATATCACCCTCAGCAGCGAGATTGGGCACGGCCGCTCCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTc                                                                <  1:274795/148‑1 (MQ=255)
                                                                                       ccATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACAGGTGCTGGTTGTCGAg                                                             >  8:280445/1‑148 (MQ=255)
                                                                                       ccATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAg                                                             >  2:63761/1‑148 (MQ=255)
                                                                                       ccATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAg                                                             >  2:179427/1‑148 (MQ=255)
                                                                                        cATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGa                                                           <  7:280445/149‑1 (MQ=255)
                                                                                             gCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGAt                                                       <  6:323963/148‑1 (MQ=255)
                                                                                                      ccTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTgc                                             >  4:2494/1‑149 (MQ=255)
                                                                                                      ccTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTgc                                             <  5:259516/149‑1 (MQ=255)
                                                                                                      ccTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCACGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTgc                                             <  5:59999/149‑1 (MQ=255)
                                                                                                          agcagcGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTgcgc                                           <  1:179427/147‑1 (MQ=255)
                                                                                                                  gaTTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCttt                                 <  3:20420/149‑1 (MQ=255)
                                                                                                                    ttGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGc                               >  6:199332/1‑149 (MQ=255)
                                                                                                                     tGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGc                               >  5:256907/1‑148 (MQ=255)
                                                                                                                     tGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGc                               >  5:252533/1‑148 (MQ=255)
                                                                                                                        ccAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGCCAgg                           >  3:25157/1‑149 (MQ=255)
                                                                                                                         cAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGCCAGgc                          <  1:189291/149‑1 (MQ=255)
                                                                                                                          aGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGTCACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCtt                                  <  4:190801/140‑1 (MQ=255)
                                                                                                                          aGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGTCACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCtt                                  >  3:190801/1‑140 (MQ=255)
                                                                                                                          aGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGGCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGCCAGgcg                         >  1:221747/1‑149 (MQ=255)
                                                                                                                                 aaCCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGCCAGGCGCTgg                     <  3:309331/149‑1 (MQ=255)
                                                                                                                                    ccAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGCCAGGCGCTGGGc                   >  8:271760/1‑148 (MQ=255)
                                                                                                                                      aGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGCCAGGCGCTGGGCGa                 >  5:185739/1‑148 (MQ=255)
                                                                                                                                       ggTGGAGCTGTATTTGTCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGCCAGGCGCTGGGc                   <  5:25496/145‑1 (MQ=255)
                                                                                                                                       ggTGGAGCTGTATTTGTCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGCCAGGCGCTGGGc                   <  5:20727/145‑1 (MQ=255)
                                                                                                                                       ggTGGAGCTGTATTTGTCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGCCAGGCGCTGGGc                   >  6:25496/1‑145 (MQ=255)
                                                                                                                                       ggTGGAGCTGTATTTGTCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGCCAGGCGCTGGGc                   >  6:20727/1‑145 (MQ=255)
                                                                                                                                              cTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGCCAGGCGCTGGGCGAAAACGGCtt        <  5:199332/149‑1 (MQ=255)
                                                                                                                                               tggaTTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGCACCGGTGCCAGCCAGGCAAGACGGCGGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGCCAGGCGCTGGGCGAAAACGGCTTc       <  6:248444/146‑1 (MQ=255)
                                                                                                                                                 tATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGCCAGGCGCTGGGCGAAAACGGCTTcc      <  2:221747/148‑1 (MQ=255)
                                                                                                                                                     tGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGCCAGGCGCTGGGCGAAAACGGCTTCCCCTg  <  4:25157/148‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
CGCCTTCGAACCCTTCTTCAGCACCAAGGCGACCGGCCAAGGCATTGGCCTGGGGTTGTCGATGGTGTACGGTTTCAGCAAGCAGTCCCATGGGCATATCACCCTCAGCAGCGAGATTGGCCAGG‑‑GAACCCAGGTGGAGCTGTATTTGCCACGCCATGTCGATCAGGCCCTGTCATACAGCGAACCGGTGCCAGCCAGGCAAGACGGCTGCAGCCGACACGTGCTGGTTGTCGAGGACGATCCCCATGTGCGCCAACTGCTTTGCCAGGCGCTGGGCGAAAACGGCTTCCCCTG  >  NC_002947/3885156‑3885449

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: