Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,259,686 A→G 12.1% intergenic (‑390/+33) psuK ← / ← fruA pseudouridine kinase/fused fructose‑specific PTS enzymes: IIBcomponent/IIC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,259,6860AG12.1% 154.0 / 12.4 59intergenic (‑390/+33)psuK/fruApseudouridine kinase/fused fructose‑specific PTS enzymes: IIBcomponent/IIC components
Reads supporting (aligned to +/- strand):  ref base A (28/23);  new base G (3/4);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

AAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAAT  >  NC_000913/2259538‑2259832
                                                                                                                                                    |                                                                                                                                                   
aaaaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacagga                                                                                                                                                     <  1:139399/149‑1 (MQ=255)
aaaaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacagga                                                                                                                                                     <  1:540452/149‑1 (MQ=255)
aaaaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGCTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacagga                                                                                                                                                     >  3:348523/1‑149 (MQ=255)
         gAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAAtt                                                                                                                                            >  1:224020/1‑149 (MQ=255)
              taatCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGc                                                                                                                                      >  4:381388/1‑149 (MQ=255)
                 tCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCt                                                                                                                                   <  1:206488/149‑1 (MQ=255)
                   tGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGt                                                                                                                                 >  1:86167/1‑149 (MQ=255)
                          gAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAacacac                                                                                                                          >  4:79611/1‑149 (MQ=255)
                                 tAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTttat                                                                                                                   <  2:446906/149‑1 (MQ=255)
                                     cTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACg                                                                                                               >  3:64158/1‑149 (MQ=255)
                                                atgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTC‑CTGTCCTGTAac                                                                                                                              >  2:15079/1‑123 (MQ=38)
                                                atgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCA‑GAAATTTCCTGTCCTGTAac                                                                                                                              <  1:15079/123‑1 (MQ=38)
                                                   aaCGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGTCAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACt                                                                                                 <  4:197969/149‑1 (MQ=255)
                                                   aaCGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACt                                                                                                 >  1:428226/1‑149 (MQ=255)
                                                   aaCGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACt                                                                                                 <  1:494811/149‑1 (MQ=255)
                                                    aCGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTg                                                                                                <  2:329117/149‑1 (MQ=255)
                                                     cGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGc                                                                                               <  2:431787/149‑1 (MQ=255)
                                                      gTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCg                                                                                              >  1:315405/1‑149 (MQ=255)
                                                       tttGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGt                                                                                             <  2:24854/149‑1 (MQ=255)
                                                           aaGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCAc                                                                                         <  4:47915/149‑1 (MQ=255)
                                                             gTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTttatt                                                                                                                  >  3:204524/1‑122 (MQ=255)
                                                             gTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTttatt                                                                                                                  <  4:204524/122‑1 (MQ=255)
                                                             ctttCTACTGAAAGATAAACCCATTTAAATTCATAAAATTACAATCACAAACCGGTCATTCTCTTTGCGCATCTTTAATATTACATTACCGGAAATTT‑CTGCCCTGTCACACATATTTTATTACGCTCCTTCCGCTACTGTGTCCACtt                                                                                       >  4:513571/2‑149 (MQ=255)
                                                                     tGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGt                                                                               <  4:73035/149‑1 (MQ=255)
                                                                     tGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGt                                                                               >  4:55926/1‑149 (MQ=255)
                                                                     tGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGt                                                                               >  4:109353/1‑149 (MQ=255)
                                                                      gAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGtt                                                                              <  3:505925/149‑1 (MQ=255)
                                                                         aGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCa                                                                           <  4:338001/149‑1 (MQ=255)
                                                                              aaaCCATTTAAATTCATAAAATTACATTTACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAg                                                                      >  4:203656/1‑149 (MQ=255)
                                                                               aaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAgg                                                                     >  3:260109/1‑149 (MQ=255)
                                                                                aCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACATGACAGGAAATTT‑CTTCCCTGTAACACACATTTTATTATGCTGCATTCGCTACTGCGTCCCCTTCCGGACGTTTCAGGAAGGc                                                                    >  4:45231/1‑149 (MQ=255)
                                                                                aCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGc                                                                    <  1:546082/149‑1 (MQ=255)
                                                                                   aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTcc                                                                                           <  4:248711/123‑1 (MQ=255)
                                                                                   aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTcc                                                                                           >  3:248711/1‑123 (MQ=255)
                                                                                   aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATa                                                                 >  2:430444/1‑149 (MQ=255)
                                                                                    tttAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAg                                                                <  4:6466/149‑1 (MQ=255)
                                                                                     ttAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAgg                                                               <  3:357603/149‑1 (MQ=255)
                                                                                         aTTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaa                                                           >  3:148388/1‑149 (MQ=255)
                                                                                          ttCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaaa                                                          <  1:219262/149‑1 (MQ=255)
                                                                                            cATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAcc                                                        >  3:62119/1‑149 (MQ=255)
                                                                                             aTAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAccc                                                       >  2:438811/1‑149 (MQ=255)
                                                                                                  aTTACATTCACAAACCGGTCTTTCTCTTTGCGCTCATCAAACGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGccac                                                  >  3:370236/1‑149 (MQ=255)
                                                                                                       aTTCACAACCCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCg                                             <  2:280244/149‑1 (MQ=255)
                                                                                                       aTTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCg                                             >  1:12859/1‑149 (MQ=255)
                                                                                                           acaAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAg                                                                      >  1:571690/1‑120 (MQ=255)
                                                                                                           acaAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAg                                                                      <  2:571690/120‑1 (MQ=255)
                                                                                                            caAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCg                                        >  1:81967/1‑149 (MQ=255)
                                                                                                              aaCCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGc                                      >  2:556891/1‑149 (MQ=255)
                                                                                                                   gTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCaataa                                 >  1:239597/1‑149 (MQ=255)
                                                                                                                     cATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAAtt                               <  3:63396/149‑1 (MQ=255)
                                                                                                                      aTTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTg                              >  3:186960/1‑149 (MQ=255)
                                                                                                                            tttGCGCTCATCAAATGTTACAGGGCAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCa                        <  2:153387/149‑1 (MQ=255)
                                                                                                                            tttGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCa                        <  2:275898/149‑1 (MQ=255)
                                                                                                                            tttGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCa                        >  1:575239/1‑149 (MQ=255)
                                                                                                                               gcgcTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGt                     <  4:348523/149‑1 (MQ=255)
                                                                                                                                 gcTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTaa                   >  2:278380/1‑149 (MQ=255)
                                                                                                                                    catcaAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAAccc                <  1:4321/149‑1 (MQ=255)
                                                                                                                                           tGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCg         <  2:224020/149‑1 (MQ=255)
                                                                                                                                             ttACAGGGCAGAAATTTC‑CTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGc       <  1:360185/149‑1 (MQ=255)
                                                                                                                                             ttACAGGGCAGAAATTTC‑CTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGc       >  4:504131/1‑149 (MQ=255)
                                                                                                                                             ttACAGGGCAGAAATTTC‑CTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGc       >  2:360185/1‑149 (MQ=255)
                                                                                                                                             ttACAGGGCAGAAATTTC‑CTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGc       <  3:504131/149‑1 (MQ=255)
                                                                                                                                                  ggacaggaAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCGTAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAAt  >  4:320742/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT‑CTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAAT  >  NC_000913/2259538‑2259832

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