Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,573,164 T→G 6.6% F168C (TTC→TGC)  pagL‑II → lipid A 3‑O‑deacylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,573,1640TG 6.6% 187.5 / 4.0 61F168C (TTC→TGC) pagL‑IIlipid A 3‑O‑deacylase
Reads supporting (aligned to +/- strand):  ref base T (25/32);  new base G (2/2);  total (27/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 = 8.93e-01

CAGCGACTGGGGTCGGCGGTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGA  >  NC_002947/3573016‑3573306
                                                                                                                                                    |                                                                                                                                              
cAGCGACTGGGGTCGGCGGTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGtt                                                                                                                                                >  2:197263/1‑149 (MQ=255)
 aGCGACTGGGGTCGGCGGTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTc                                                                                                                                               <  7:25567/149‑1 (MQ=255)
 aGCGACTGGGGTCGGCGGTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTc                                                                                                                                               <  5:216054/149‑1 (MQ=255)
    gACTGGGGTCGGCGGTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTac                                                                                                                                            <  3:108134/149‑1 (MQ=255)
    gACTGGGGTCGGCGGTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTa                                                                                                                                             >  6:250798/1‑148 (MQ=255)
    gACTGGGGTCGGCGGTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTa                                                                                                                                             <  5:250798/148‑1 (MQ=255)
      cTGGGGTCGGCGGTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTacag                                                                                                                                          >  8:347546/1‑149 (MQ=255)
      cTGGGGTCGGCGGTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTacag                                                                                                                                          <  3:334731/149‑1 (MQ=255)
          ggTCGGCGGTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTTCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTacagcaca                                                                                                                                      <  1:257235/149‑1 (MQ=255)
                cggTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTc                                                                                                                                <  1:197263/149‑1 (MQ=255)
                  gTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAACGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCt                                                                                                                               <  1:289570/148‑1 (MQ=255)
                  gTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCt                                                                                                                               >  2:289570/1‑148 (MQ=255)
                          tGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTAc                                                                                                                       <  6:336428/148‑1 (MQ=255)
                            aaGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACgg                                                                                                                     <  7:281668/148‑1 (MQ=255)
                            aaGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACgg                                                                                                                     <  1:351549/148‑1 (MQ=255)
                                   tATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCAAAGAGTCGTACTCGCTGTTCTacagcaca                                                                                                                                      <  6:16870/124‑1 (MQ=255)
                                   tATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCAAAGAGTCGTACTCGCTGTTCTacagcaca                                                                                                                                      >  5:16870/1‑124 (MQ=255)
                                         cTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGc                                                                                                       <  7:93323/149‑1 (MQ=255)
                                              gCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCgg                                                                                                  <  6:212131/149‑1 (MQ=255)
                                               ccTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCgg                                                                                                  >  1:224628/1‑148 (MQ=255)
                                                cTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTc                                                                                                                                               <  5:100540/102‑1 (MQ=255)
                                                cTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTc                                                                                                                                               >  6:100540/1‑102 (MQ=255)
                                                     gTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAAtt                                                                                             <  3:6326/147‑1 (MQ=255)
                                                     gTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAAtt                                                                                             >  4:6326/1‑147 (MQ=255)
                                                     gTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCa                                                                                           <  3:357626/149‑1 (MQ=255)
                                                         ccGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGCGGCCACTGCGCGGCCCCGGCAATTCAGaa                                                                                        <  5:240487/148‑1 (MQ=255)
                                                          cGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAgg                                                                                      <  8:46665/149‑1 (MQ=255)
                                                                  cAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGa                                                                               >  6:135652/1‑148 (MQ=255)
                                                                    gAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAAc                                                                            >  3:368526/1‑149 (MQ=255)
                                                                        gTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCg                                                                        <  7:179819/149‑1 (MQ=255)
                                                                           ggCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGc                                                                      <  1:232666/148‑1 (MQ=255)
                                                                              aTCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTgcgc                                                                                                          <  2:272836/109‑1 (MQ=255)
                                                                              aTCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTgcgc                                                                                                          >  1:272836/1‑109 (MQ=255)
                                                                              aTCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTg                                                                   >  4:372106/1‑148 (MQ=255)
                                                                                ccGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCAc                                                                >  2:91098/1‑149 (MQ=255)
                                                                                      cATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTgcg                                                          >  4:215550/1‑149 (MQ=255)
                                                                                        tGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGcgcg                                                        >  3:287213/1‑149 (MQ=255)
                                                                                         gCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCgg                                                       >  1:186463/1‑149 (MQ=255)
                                                                                           actactCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGc                                                                                                 >  4:30984/1‑105 (MQ=255)
                                                                                           actactCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGc                                                                                                 <  3:30984/105‑1 (MQ=255)
                                                                                           actactCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGt                                                      <  2:186463/148‑1 (MQ=255)
                                                                                                  cAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCggg                                              >  5:47674/1‑149 (MQ=255)
                                                                                                  cAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGC‑‑GCACCGCGCACGGTgcaggccg                                              <  1:330946/147‑7 (MQ=255)
                                                                                                      gCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGtat                                          <  7:347546/149‑1 (MQ=255)
                                                                                                      gCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGta                                           >  1:101239/1‑148 (MQ=255)
                                                                                                           gATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATa                                     <  5:135652/149‑1 (MQ=255)
                                                                                                               aaGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCt                                  <  3:372106/148‑1 (MQ=255)
                                                                                                                 gcagcCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGtt                               >  1:308792/1‑149 (MQ=255)
                                                                                                                 gcagcCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGtt                               <  2:101239/149‑1 (MQ=255)
                                                                                                                    gcCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGaa                            >  8:701/1‑149 (MQ=255)
                                                                                                                     ccAAACGGCGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAAc                           >  4:142318/1‑149 (MQ=255)
                                                                                                                               gCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGTGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGt                 <  3:215550/149‑1 (MQ=255)
                                                                                                                                cATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGt                 >  3:65501/1‑148 (MQ=255)
                                                                                                                                cATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGt                 >  5:49991/1‑148 (MQ=255)
                                                                                                                                 aTAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCa               >  6:157906/1‑149 (MQ=255)
                                                                                                                                       tCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATAcccc         <  2:224628/149‑1 (MQ=255)
                                                                                                                                          agaacGCTGTGCTGTAGAACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTc                   <  6:220699/132‑1 (MQ=255)
                                                                                                                                          agaacGCTGTGCTGTAGAACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTc                   >  5:220699/5‑136 (MQ=255)
                                                                                                                                            cTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttgtt    >  6:219476/1‑149 (MQ=255)
                                                                                                                                            cTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttgt     <  4:65501/148‑1 (MQ=255)
                                                                                                                                              cGCTGTGCTGTAGAACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttgttg   >  1:201171/1‑148 (MQ=255)
                                                                                                                                              cGCTGTGCTGTAGAACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGa  <  2:201171/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CAGCGACTGGGGTCGGCGGTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGA  >  NC_002947/3573016‑3573306

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