Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 881,214 C→T 6.9% Q400* (CAG→TAG)  lapB → Lipopolysaccharide assembly protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB143881,2140CT 6.9% 188.6 / 3.5 58Q400* (CAG→TAG) lapBLipopolysaccharide assembly protein B
Reads supporting (aligned to +/- strand):  ref base C (19/35);  new base T (2/2);  total (21/37)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

AAGGCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGACGGATGACCCCAGCGTGCTCAACGC  >  ACB143/881082‑881338
                                                                                                                                    |                                                                                                                            
aaGGCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAgg                                                                                                                           <  4:236709/136‑1 (MQ=255)
aaGGCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAgg                                                                                                                           <  8:163993/136‑1 (MQ=255)
 aGGCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAgg                                                                                                                           <  4:117828/135‑1 (MQ=255)
  ggCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCg                                                                                                                         <  3:236887/136‑1 (MQ=255)
     tGCCAGCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATg                                                                                                                      >  6:22636/1‑136 (MQ=255)
          gCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGTAGGCGATgggg                                                                                                                   >  4:126432/1‑134 (MQ=255)
          gCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGTAGGCGATgggg                                                                                                                   >  4:122031/1‑134 (MQ=255)
          gCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATgggg                                                                                                                   >  8:223131/1‑134 (MQ=255)
          gCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATgggg                                                                                                                   >  6:207351/1‑134 (MQ=255)
          gCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGAAgggg                                                                                                                   >  2:88579/1‑134 (MQ=255)
                agcagGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAgc                                                                                                           <  1:126925/136‑1 (MQ=255)
                        ggCGACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCa                                                                                                   >  2:187131/1‑136 (MQ=255)
                           gACCCGGAAACGGCCGGCATGTTGACATAACAGGCGTTTCAGCGTCTGGACCAGATGAAGCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCg                                                                                                <  5:130116/136‑1 (MQ=255)
                           gACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCg                                                                                                <  2:19937/136‑1 (MQ=255)
                           gACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCg                                                                                                <  4:295008/136‑1 (MQ=255)
                           gACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCg                                                                                                <  5:155477/136‑1 (MQ=255)
                           gACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCg                                                                                                <  8:36365/136‑1 (MQ=255)
                                     cggccggcAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTgaa                                                                                      <  3:170087/136‑1 (MQ=255)
                                     cggccggcAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTgaa                                                                                      <  7:66615/136‑1 (MQ=255)
                                              aGTTGACAGAACATGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACtgtgt                                                                             <  3:155981/136‑1 (MQ=255)
                                              aGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACtgtgt                                                                             <  3:107765/136‑1 (MQ=255)
                                               gTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACtgtgtg                                                                            <  1:263572/136‑1 (MQ=255)
                                                          aGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATgg                                                                 <  6:156350/136‑1 (MQ=255)
                                                          aGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATgg                                                                 <  5:30442/136‑1 (MQ=255)
                                                          aGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATgg                                                                 >  5:56356/1‑136 (MQ=255)
                                                           ggggATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCGCTGACTGTCGCTGGGGAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGc                                                                <  6:52236/133‑1 (MQ=255)
                                                           ggCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGGTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGc                                                                <  3:188368/136‑1 (MQ=255)
                                                           ggCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGc                                                                <  6:119811/136‑1 (MQ=255)
                                                           ggCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGc                                                                <  6:154157/136‑1 (MQ=255)
                                                           ggCGATTCAGCGGCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTGGGCCCTGACTGGCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGc                                                                <  6:237499/136‑1 (MQ=255)
                                                            gCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGc                                                                <  5:186245/135‑1 (MQ=255)
                                                                   aGCGTCTGGACATGTTGAACCTGTTCTTCTCTGGCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGGGACGGGGCTGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGAcc                                                         <  2:295783/135‑1 (MQ=255)
                                                                    gCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGAccc                                                        >  2:13856/1‑135 (MQ=255)
                                                                    gCGTCTGGACAAGATGAACCAGTTCTTCTCTGTAGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGACGCCATCGGCGTGCTGAAGAACTGTGTGTAAAGTTATGGCGATGAccc                                                        >  1:50839/1‑135 (MQ=255)
                                                                       tCTGGTCAAGATGATCCAGTTGTTCTCTGTCGAGTCGGCGCTGAGTGTCGCTGCGCAATTGCAGGCGATGGGGCTGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCatg                                                    <  1:141402/136‑3 (MQ=255)
                                                                               aGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGGCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGga                                            <  1:249811/136‑1 (MQ=255)
                                                                                      ccAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGTAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTg                                     <  4:179415/136‑1 (MQ=255)
                                                                                        aGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTggg                                   >  7:114183/1‑136 (MQ=255)
                                                                                             ttctCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAATTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGc                              >  4:213802/1‑136 (MQ=255)
                                                                                                  tGTCGATTCGGCCCTGTCGGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGACg                         <  4:182777/136‑1 (MQ=255)
                                                                                                  tGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGACg                         >  1:60634/1‑136 (MQ=255)
                                                                                                   gTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGACg                         <  8:114183/135‑1 (MQ=255)
                                                                                                    tCGATTCGGCCCTGACTGTCGCTGCGCTATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGACg                         <  4:159632/134‑1 (MQ=255)
                                                                                                    tCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGTAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGAGCCCAAGGTCATGGAGAAAGTGGGCAAGCTGAc                          <  6:244194/133‑1 (MQ=255)
                                                                                                    tCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGAc                          <  3:100215/133‑1 (MQ=255)
                                                                                                     cGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGAc                          >  4:229764/1‑132 (MQ=255)
                                                                                                     cGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGAc                          <  3:229764/132‑1 (MQ=255)
                                                                                                     cGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGACGGAt                      >  7:131178/1‑136 (MQ=255)
                                                                                                            gCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGACGGATGAccc                 >  4:276493/1‑134 (MQ=255)
                                                                                                            gCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGACGGATGACCCCa               <  1:65859/136‑1 (MQ=255)
                                                                                                                  aCTGTCGCTGCGCAATTGCAGGCGAGGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGACGGATGACCCCAGCGTGc         <  1:177337/136‑1 (MQ=255)
                                                                                                                   cTGTCGCTGCGCAATTGTAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGACGGATGACCCCAGCGTGCt        <  3:122031/136‑1 (MQ=255)
                                                                                                                   cTGTCGCTGCGCAATTGTAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGACGGATGACCCCAGCGTGCt        <  3:126432/136‑1 (MQ=255)
                                                                                                                     gTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGACGGATGACCCCAGCGTGCTc       >  4:107431/1‑135 (MQ=255)
                                                                                                                     gTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGACGGATGACCCCAGCGTGCTc       >  1:29602/1‑135 (MQ=255)
                                                                                                                         cTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCCTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGACGGATGCCCCCAGCGTGCTCACCGc  >  5:187029/1‑136 (MQ=255)
                                                                                                                          tGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGACGGATGACCCCAGCGTGCTCAACGc  >  2:125509/1‑135 (MQ=255)
                                                                                                                          tGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGACGGATGACCCCAGCGTGCTCAACGc  >  8:1308/1‑135 (MQ=255)
                                                                                                                                    |                                                                                                                            
AAGGCTGCCAGCCTGCAGCAGGCCGGCGACCCGGAAACGGCCGGCAAGTTGACAGAACAGGCGATTCAGCGTCTGGACAAGATGAACCAGTTCTTCTCTGTCGATTCGGCCCTGACTGTCGCTGCGCAATTGCAGGCGATGGGGCAGGAAGCCGCCGCCATCGGCGTGCTGAAGAACTGTGTGGAAAGTTATGGCGATGACCCCAAGGTCATGGAGAAAGTGGGCAAGCTGACGGATGACCCCAGCGTGCTCAACGC  >  ACB143/881082‑881338

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