Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,149,184 G→A 5.4% A361V (GCG→GTG)  eptC ← LPS heptose I phosphoethanolamine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,149,1840GA 5.4% 237.7 / 2.2 74A361V (GCG→GTG) eptCLPS heptose I phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base G (38/32);  new base A (2/2);  total (40/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 = 2.98e-01

CATCAAACTTGCCCTGATTTTCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCA  >  NC_000913/4149041‑4149328
                                                                                                                                               |                                                                                                                                                
cATCAAACTTGCCCTGATTTTCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGCTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCActc                                                                                                                                             >  5:335535/1‑149 (MQ=255)
   cAAACTAGCCCTGATTTTCCGGGAAGCGGTATTAGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCAATCACCACATCCAGGAACGGCTTCAGCACGTTGGGTTCGTATTCACGCGCACTCAGc                                                                                                                                          >  4:81345/1‑149 (MQ=255)
          gCCCTGATTTTCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGt                                                                                                                                    <  3:184256/148‑1 (MQ=255)
            ccTGATTTTCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTTGTATTCACGCGCACTCTGCGTACGTTGc                                                                                                                                 >  8:46854/1‑149 (MQ=255)
            ccTGATTTTCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGc                                                                                                                                 <  7:60456/149‑1 (MQ=255)
             cTGATTTTCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCt                                                                                                                                >  5:211740/1‑149 (MQ=255)
             cTGATTTTCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCt                                                                                                                                >  1:53298/1‑149 (MQ=255)
                aTTTTCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTg                                                                                                                               <  2:331253/147‑1 (MQ=255)
                aTTTTCCGGGGAGCGGTATTTGTATTTTATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTg                                                                                                                               <  8:148900/147‑1 (MQ=255)
                  tttCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCACACTCTGCGTACGTTGCTGGtt                                                                                                                            >  3:286666/1‑148 (MQ=255)
                   ttCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTc                                                                                                                           <  1:87495/148‑1 (MQ=255)
                            gCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTAc                                                                                                                  <  6:211740/148‑1 (MQ=255)
                              ggTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGc                                                                                                               <  2:15190/149‑1 (MQ=255)
                              ggTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCACACTCTGCGTACGTTGCTGGTTCATGTAGTACTGc                                                                                                               <  4:286666/149‑1 (MQ=255)
                                tatttgtatttgATATGCGTACCGAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTGCTGGAAGGGCTTCAGCACGTTGGTGTCGTATTGACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCtt                                                                                                             <  3:272786/149‑1 (MQ=255)
                                    tgtatttgATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCg                                                                                                         >  5:211709/1‑149 (MQ=255)
                                     gtatttgATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCgg                                                                                                        >  1:229999/1‑149 (MQ=255)
                                        tttgATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACg                                                                                                                                     <  2:347684/117‑1 (MQ=255)
                                        tttgATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACg                                                                                                                                     >  1:347684/1‑117 (MQ=255)
                                        tttgATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCt                                                                                                     >  1:386345/1‑149 (MQ=255)
                                        tttgATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCt                                                                                                     >  4:59162/1‑149 (MQ=255)
                                           gATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACg                                                                                                                                     <  4:95054/114‑1 (MQ=255)
                                           gATATGCGTAACCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCAATTCCTGGAACGGCTTCAGCAATTTGGTGTCGTATTCACGCGCACTCTTCTTACGTTGCTGGTTCATGTAGTACTGCTTGTCGTTCTgtc                                                                                                  >  2:309465/1‑147 (MQ=255)
                                              atGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGgcgc                                                                                                >  5:277372/1‑148 (MQ=255)
                                                 cGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTCCTGGTTCATgtag                                                                                                                     <  7:130119/124‑1 (MQ=255)
                                                 cGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTACTGGTTCATgtag                                                                                                                     >  8:130119/1‑124 (MQ=255)
                                                   tACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGaaaa                                                                                           >  7:207167/1‑148 (MQ=255)
                                                    aCCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAACTa                                                                                         >  3:67852/1‑149 (MQ=255)
                                                     cccAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATa                                                                                         <  2:229999/148‑1 (MQ=255)
                                                      ccAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGc                                                                                                               <  7:64025/125‑1 (MQ=255)
                                                      ccAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGc                                                                                                               >  8:64025/1‑125 (MQ=255)
                                                        agcagATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGt                                                                                     <  8:367155/149‑1 (MQ=255)
                                                             aTGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCa                                                                                >  6:110972/1‑149 (MQ=255)
                                                              tGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCAt                                                                               >  5:272879/1‑149 (MQ=255)
                                                                aaCAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATgg                                                                             <  3:59162/149‑1 (MQ=255)
                                                                aaCAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATg                                                                              <  2:386345/148‑1 (MQ=255)
                                                                   aaTGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTAt                                                                          >  7:318756/1‑149 (MQ=255)
                                                                    aTGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTAtt                                                                         <  8:207167/149‑1 (MQ=255)
                                                                     tGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCACACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCa                                                                                <  4:102233/141‑1 (MQ=255)
                                                                     tGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCACACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCAtagatcgg                                                                        >  3:102233/1‑142 (MQ=255)
                                                                      taTCAGTTTCTTCGTCGGTGTGTCATTCAGCACTTCAGGGTACGGCTTGATCTCGTTGGTGTCTTATTCTCGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCTGTCTGGCGGGTAATTAGGTTCAGGATGGTATTg                                                                        <  4:105798/147‑1 (MQ=255)
                                                                          aGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTgcgggc                                                                   >  8:165665/1‑149 (MQ=255)
                                                                           gTTTGTTCGGGGCAGGGTCATTCAGCACTTCCTGGATCGGCTGCAGCACGTTGGTGTCGTATTCTGGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTTCTGCTTGTCGGTCTGGCGCGATAATACCGTCAGCATGGTATTgcgggcg                                                                  <  3:81345/149‑1 (MQ=255)
                                                                           gTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTgcgggcg                                                                  <  4:67852/149‑1 (MQ=255)
                                                                           gTTTCTTCGGCGCAGGGGCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATCCCGTCAGCATGGTATTgcgggcg                                                                  <  2:114305/149‑1 (MQ=255)
                                                                            tttCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTgcgggcgg                                                                 <  1:320238/149‑1 (MQ=255)
                                                                                        aGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTAATGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTctgct                                                     >  7:85539/1‑149 (MQ=255)
                                                                                         gggTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTctgctg                                                    <  2:53298/149‑1 (MQ=255)
                                                                                           gTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCtggt                                                  >  6:117772/1‑149 (MQ=255)
                                                                                              aTTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGttggttgg                                               <  6:211709/149‑1 (MQ=255)
                                                                                                  aGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACTCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGAt                                           >  1:2243/1‑149 (MQ=255)
                                                                                                     aCTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATc                                          >  7:402253/1‑147 (MQ=255)
                                                                                                       ttCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAg                                       >  5:222084/1‑148 (MQ=255)
                                                                                                           tGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGt                                  >  1:373931/1‑149 (MQ=255)
                                                                                                             gAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGttt                                >  8:4781/1‑149 (MQ=255)
                                                                                                             gAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGATCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTAAGCATGTTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGgtt                                >  5:121421/1‑149 (MQ=255)
                                                                                                                         cACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCa                                        >  2:58395/1‑129 (MQ=255)
                                                                                                                         cACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCa                                        <  1:58395/129‑1 (MQ=255)
                                                                                                                         cACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCt                    <  4:279637/149‑1 (MQ=255)
                                                                                                                         cACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGGCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCt                    >  5:24956/1‑149 (MQ=255)
                                                                                                                            gTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGtt                 <  6:277372/149‑1 (MQ=255)
                                                                                                                                gtgtCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTtcat             <  7:165665/149‑1 (MQ=255)
                                                                                                                                gtgtCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTtcat             >  3:147235/1‑149 (MQ=255)
                                                                                                                                 tgtCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAAt                                                                                          >  3:220553/1‑71 (MQ=255)
                                                                                                                                 tgtCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAAt                                                                                          <  4:220553/71‑1 (MQ=255)
                                                                                                                                   tCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTtc               >  6:350085/1‑144 (MQ=255)
                                                                                                                                   tCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTtc               <  5:350085/144‑1 (MQ=255)
                                                                                                                                   tCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTcatcat          >  1:214580/1‑149 (MQ=255)
                                                                                                                                     gTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGt        >  7:231067/1‑149 (MQ=255)
                                                                                                                                      tATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGGCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGtt       <  5:110972/149‑1 (MQ=255)
                                                                                                                                       aTTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGtt       >  3:251045/1‑148 (MQ=255)
                                                                                                                                         tCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTtcat    >  3:306171/1‑149 (MQ=255)
                                                                                                                                          cACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTtcat    >  1:366399/1‑148 (MQ=255)
                                                                                                                                           aCGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGTTTTTATAACCCGCCTGTTTCATCATGTtcatca  <  2:2243/149‑1 (MQ=255)
                                                                                                                                            cgcgCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTtca     >  7:214903/1‑145 (MQ=255)
                                                                                                                                            cgcgCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTtca     <  8:214903/145‑1 (MQ=255)
                                                                                                                                             gcgcACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGt        <  1:275568/141‑1 (MQ=255)
                                                                                                                                             gcgcACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGt        >  2:275568/1‑141 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CATCAAACTTGCCCTGATTTTCCGGGTAGCGGTATTTGTATTTGATATGCGTACCCAGCAGATGAACAATGATCAGTTTCTTCGGCGCAGGGTCATTCAGCACTTCCTGGAACGGCTTCAGCACGTTGGTGTCGTATTCACGCGCACTCTGCGTACGTTGCTGGTTCATGTAGTACTGCTTGTCGGTCTGGCGCGAAAATACCGTCAGCATGGTATTGCGGGCGGTCATCGTCTGCTGGTTGGTGATCCAGAAGGTTTTATAACCCGCCTGTTTCATCATGTTCATCA  >  NC_000913/4149041‑4149328

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: