Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 491,089 G→A 6.9% N444N (AAC→AAT) ‡ lptD ← LPS‑assembly protein LptD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947491,0890GA 6.9% 258.5 / 5.6 86N444N (AAC→AAT) ‡lptDLPS‑assembly protein LptD
Reads supporting (aligned to +/- strand):  ref base G (36/44);  new base A (2/4);  total (38/48)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

CGCCAGGCTGGAACGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCATCCAGCGCTGAT  >  NC_002947/490941‑491237
                                                                                                                                                    |                                                                                                                                                    
cGCCAGGCTGGAACGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTgg                                                                                                                                                      <  4:390305/149‑1 (MQ=255)
   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                                                                                                                              <  5:208203/149‑1 (MQ=255)
                        cGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACg                                                                                                                              >  6:102371/1‑149 (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                                                                                                                                         >  4:187049/1‑129 (MQ=255)
                                 tCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTcgc                                                                                                                                         <  3:187049/129‑1 (MQ=255)
                                 tCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTgg                                                                                                                     >  6:134293/1‑149 (MQ=255)
                                   aGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTc                                                                                                                   <  2:324451/149‑1 (MQ=255)
                                   aGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTc                                                                                                                   >  7:266194/1‑149 (MQ=255)
                                   aGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTc                                                                                                                   <  2:359247/149‑1 (MQ=255)
                                   aGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTc                                                                                                                   <  2:202747/149‑1 (MQ=255)
                                           cAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCa                                                                                                           <  5:405939/149‑1 (MQ=255)
                                            aGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTc                                                                                                             >  5:307248/1‑146 (MQ=255)
                                            aGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCa                                                                                                           <  6:307248/148‑1 (MQ=255)
                                            aGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAg                                                                                                          <  5:155119/149‑1 (MQ=255)
                                            aGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAg                                                                                                          >  2:77979/1‑149 (MQ=255)
                                              gCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTc                                                                                                                   >  7:67109/1‑138 (MQ=255)
                                              gCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTc                                                                                                                   <  8:67109/138‑1 (MQ=255)
                                              gCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGt                                                                                                        >  3:355990/1‑149 (MQ=255)
                                              gCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGt                                                                                                        <  6:144631/149‑1 (MQ=255)
                                              gCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGt                                                                                                        <  6:349151/149‑1 (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)
                                                                                                                                                    |                                                                                                                                                    
CGCCAGGCTGGAACGGCAGCGTACCGTTCAGGGTCAGCTGCGGCAGGCGATCATACGGCGTGATCTGCGAGATGGTCGCCATCTCGTAGGCATGCACATTCAACCGCGCGGTGTAGCTGTCGCCTCGGTAGGTCAACGCACCCTGCTGGTTCAGCAGGTCGCGGCTTTCCACGCCAATCTGGTCGGATTCCAGGTCCTGGAAGTAGAACGGGTCGCTGATGTCGGTGTAGTCCACCTCGGTCATCAGGCGTTCATCCAGGCCACCCTTGTGCTGCCAGTTGACCATCCAGCGCTGAT  >  NC_002947/490941‑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: ‑