Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 491,091 T→C 6.6% N444D (AAC→GAC) ‡ lptD ← LPS‑assembly protein LptD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947491,0910TC 6.6% 258.2 / 6.6 90N444D (AAC→GAC) ‡lptDLPS‑assembly protein LptD
Reads supporting (aligned to +/- strand):  ref base T (39/45);  new base C (2/4);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01

CAGGCTGGAACGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCATCCAGCGCTGAT  >  NC_002947/490944‑491237
                                                                                                                                                   |                                                                                                                                                  
cAGGCTGGAACGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGtt                                                                                                                                                    >  4:89765/1‑148 (MQ=255)
  ggCTGGAACGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTcag                                                                                                                                                 >  1:104409/1‑149 (MQ=255)
    cTGGAACGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTcag                                                                                                                                                 >  6:208203/1‑147 (MQ=255)
        aaCGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGAACCAGCAgggt                                                                                                                                           <  7:81954/149‑3 (MQ=255)
        aaCGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGAACCAGCAgggt                                                                                                                                           <  8:32569/149‑3 (MQ=255)
        aaCGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGAACCAGCAgggt                                                                                                                                           <  2:375470/149‑3 (MQ=255)
        aaCGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGAACCAGCAgggt                                                                                                                                           <  2:329409/149‑3 (MQ=255)
           ggCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTcgcg                                                                                                                                        <  3:183841/149‑1 (MQ=255)
            gCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCgg                                                                                                                                       <  6:20075/149‑1 (MQ=255)
             cAGCGTACCGTTCAGGATCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGc                                                                                                                                      >  1:381887/1‑149 (MQ=255)
                     cGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACg                                                                                                                              >  6:102371/1‑149 (MQ=255)
                     cGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACg                                                                                                                              <  5:208203/149‑1 (MQ=255)
                      gTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGc                                                                                                                             <  1:335181/149‑1 (MQ=255)
                       ttCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGc                                                                                                                             >  3:16757/1‑148 (MQ=255)
                           gggTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGTTTTCCACGCCAATc                                                                                                                        <  5:397400/149‑1 (MQ=255)
                            ggTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCt                                                                                                                       >  5:352338/1‑149 (MQ=255)
                              tCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTcgc                                                                                                                                         <  3:187049/129‑1 (MQ=255)
                              tCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTcgc                                                                                                                                         >  4:187049/1‑129 (MQ=255)
                              tCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTgg                                                                                                                     >  6:134293/1‑149 (MQ=255)
                                aGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTc                                                                                                                   <  2:202747/149‑1 (MQ=255)
                                aGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTc                                                                                                                   <  2:324451/149‑1 (MQ=255)
                                aGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTc                                                                                                                   <  2:359247/149‑1 (MQ=255)
                                aGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTc                                                                                                                   >  7:266194/1‑149 (MQ=255)
                                        cAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCa                                                                                                           <  5:405939/149‑1 (MQ=255)
                                         aGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTc                                                                                                             >  5:307248/1‑146 (MQ=255)
                                         aGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCa                                                                                                           <  6:307248/148‑1 (MQ=255)
                                         aGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAg                                                                                                          >  2:77979/1‑149 (MQ=255)
                                         aGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAg                                                                                                          <  5:155119/149‑1 (MQ=255)
                                           gCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTc                                                                                                                   >  7:67109/1‑138 (MQ=255)
                                           gCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTc                                                                                                                   <  8:67109/138‑1 (MQ=255)
                                           gCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGt                                                                                                        <  6:349151/149‑1 (MQ=255)
                                           gCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGt                                                                                                        <  6:144631/149‑1 (MQ=255)
                                           gCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGt                                                                                                        >  3:355990/1‑149 (MQ=255)
                                           gCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCGCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCa                                                                                                           >  5:250647/1‑146 (MQ=255)
                                           gCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCCCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCa                                                                                                           <  6:250647/146‑1 (MQ=255)
                                             gATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTcc                                                                                                      >  5:383422/1‑149 (MQ=255)
                                                 aTACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTgg                                                                                                   <  2:67636/148‑1 (MQ=255)
                                                    cGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAg                                                                                                <  8:139656/148‑1 (MQ=255)
                                                      gCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCATGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAg                                                                                             <  5:87839/149‑1 (MQ=255)
                                                           aTCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACg                                                                                         <  4:16757/148‑1 (MQ=255)
                                                           aTCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACg                                                                                         <  5:139720/148‑1 (MQ=255)
                                                             ctacGAGATGGCGACCATCTCGCAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGt                                                                                      >  8:369199/4‑149 (MQ=255)
                                                                cGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGc                                                                                   >  7:119668/1‑149 (MQ=255)
                                                                   gaTGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGACCTGGAAGTAGAACGGGTCGCTGa                                                                                >  3:167412/1‑149 (MQ=255)
                                                                       gTCGCCATCTCGTAGGCATGCAGATTCAACCCACCAGAGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCCGCTTTCCACGCCAATCTCGTCGGATTCCAGGTCCTGGAAGTAGCACGGGTCGCTGATGTc                                                                            >  7:363203/1‑149 (MQ=255)
                                                                        tCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCg                                                                           >  3:245228/1‑149 (MQ=255)
                                                                        tCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCg                                                                           <  2:381887/149‑1 (MQ=255)
                                                                          gCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGgt                                                                         <  5:102371/149‑1 (MQ=255)
                                                                          gCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGgt                                                                         <  4:355990/149‑1 (MQ=255)
                                                                           ccATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGgtg                                                                        <  8:266194/149‑1 (MQ=255)
                                                                             aTCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCATCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTa                                                                      <  8:119668/149‑1 (MQ=255)
                                                                             aTCTCGTAGGCATGCAACTTCAACCGCGCGGTGTAGCTGGCGCCTCTGTAGGGCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTa                                                                      <  4:97911/149‑1 (MQ=255)
                                                                                  gTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAg                                                                     <  2:104409/145‑1 (MQ=255)
                                                                                  gTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAg                                                                     <  6:289935/145‑1 (MQ=255)
                                                                                  gTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCa                                                                 >  6:218127/1‑149 (MQ=255)
                                                                                  gTAGGCACGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAg                                                                                                <  7:234593/118‑1 (MQ=255)
                                                                                  gTAGGCACGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAg                                                                                                >  8:234593/1‑118 (MQ=255)
                                                                                            acaTTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGt                                                                             <  7:92939/127‑1 (MQ=255)
                                                                                            acaTTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGt                                                                             >  8:92939/1‑127 (MQ=255)
                                                                                               ttCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAg                                                    >  4:199942/1‑149 (MQ=255)
                                                                                                  aaCCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTATAACGGGTCGCTGATGTCGGTGTAGTCAACCTCGGTCATCAGGct                                                 >  6:317792/1‑148 (MQ=255)
                                                                                                  aaCCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGGCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCg                                                 <  6:352338/149‑1 (MQ=255)
                                                                                                    ccGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGtt                                               >  4:193086/1‑149 (MQ=255)
                                                                                                    ccGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGCAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGtt                                               <  1:316949/149‑1 (MQ=255)
                                                                                                       cgcgGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCCGCAGGTCGCGGCTTTCCACGCCAATCTGGGCGGATTCCAGGTCCGGGACGGAGACCGGGTCGCTGATGTCGGTGTAGTGCAGCTCGGTCATCAGGCGtggcg                                            >  6:321095/1‑145 (MQ=255)
                                                                                                         cgGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATcc                                          <  3:89765/149‑1 (MQ=255)
                                                                                                           gtgtAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAg                                        <  5:218127/149‑1 (MQ=255)
                                                                                                             gtAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGc                                      <  6:383422/149‑1 (MQ=255)
                                                                                                               aGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCa                                    >  7:75511/1‑149 (MQ=255)
                                                                                                                 cTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCAcc                                  >  4:188508/1‑149 (MQ=255)
                                                                                                                    tCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCtt                               >  3:30603/1‑149 (MQ=255)
                                                                                                                    tCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCtt                               >  5:245385/1‑149 (MQ=255)
                                                                                                                      gCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTtg                              <  4:337657/148‑1 (MQ=255)
                                                                                                                       ccTCGGTAGGTCAACGCACCCTGCTGGTTGAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTtgtg                            <  5:134293/149‑1 (MQ=255)
                                                                                                                       ccTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTtgtg                            <  7:148497/149‑1 (MQ=255)
                                                                                                                       ccTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTtgtg                            <  7:369199/149‑1 (MQ=255)
                                                                                                                               ggTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAg                     <  5:264636/148‑1 (MQ=255)
                                                                                                                                gTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCTGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGtt                   <  5:321095/149‑1 (MQ=255)
                                                                                                                                  cAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGa                 >  5:35894/1‑149 (MQ=255)
                                                                                                                                     cGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCa              >  3:383088/1‑149 (MQ=255)
                                                                                                                                      gcgaccTGCTGAACCAGCAGGGTGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTTGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCAt             >  2:239315/5‑149 (MQ=255)
                                                                                                                                      gcgaccTGCTGAACCAGCAGGGTGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGATCAt             >  4:403202/5‑149 (MQ=255)
                                                                                                                                        aCCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGAATCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCATcc           <  3:193086/149‑1 (MQ=255)
                                                                                                                                           ctgctgGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTtg                              <  2:395287/127‑1 (MQ=255)
                                                                                                                                           ctgctgGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTtg                              >  1:395287/1‑127 (MQ=255)
                                                                                                                                            tgctgGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCATCAAgcg       >  7:77676/1‑149 (MQ=255)
                                                                                                                                                ggTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGtt                   >  5:214432/1‑133 (MQ=255)
                                                                                                                                                ggTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGtt                   <  6:214432/133‑1 (MQ=255)
                                                                                                                                                ggTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCATCCAGCGCTGa   >  2:406940/1‑149 (MQ=255)
                                                                                                                                                 gTTCAGCAGGTCGCGGCTTTCCACGGCAATCTGGTCGGATTCCAGGTACTGGAAGTAGAACGGGTCGCTGATGTGGGTGTAGTCGACCGCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTCCGAGTTGACCATCCAGCGCTGAt  >  4:15331/1‑149 (MQ=255)
                                                                                                                                                 gTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCATCCAGCGCTGAt  >  2:304903/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGGCTGGAACGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCATCCAGCGCTGAT  >  NC_002947/490944‑491237

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