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

AAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAAT  >  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)
              taatCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGc                                                                                                                                      >  4:381388/1‑149 (MQ=255)
                 tCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCt                                                                                                                                   <  1:206488/149‑1 (MQ=255)
                   tGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGt                                                                                                                                 >  1:86167/1‑149 (MQ=255)
                          gAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAacacac                                                                                                                          >  4:79611/1‑149 (MQ=255)
                                 tAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTttat                                                                                                                   <  2:446906/149‑1 (MQ=255)
                                     cTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACg                                                                                                               >  3:64158/1‑149 (MQ=255)
                                                atgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTCCTGTCCTGTAac                                                                                                                              >  2:15079/1‑123 (MQ=38)
                                                atgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTCCTGTCCTGTAac                                                                                                                              <  1:15079/123‑1 (MQ=38)
                                                   aaCGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGTCAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACt                                                                                                 <  4:197969/149‑1 (MQ=255)
                                                   aaCGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACt                                                                                                 >  1:428226/1‑149 (MQ=255)
                                                   aaCGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACt                                                                                                 <  1:494811/149‑1 (MQ=255)
                                                    aCGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTg                                                                                                <  2:329117/149‑1 (MQ=255)
                                                     cGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGc                                                                                               <  2:431787/149‑1 (MQ=255)
                                                      gTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCg                                                                                              >  1:315405/1‑149 (MQ=255)
                                                       tttGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGt                                                                                             <  2:24854/149‑1 (MQ=255)
                                                           aaGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCAc                                                                                         <  4:47915/149‑1 (MQ=255)
                                                             gTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTttatt                                                                                                                  >  3:204524/1‑122 (MQ=255)
                                                             gTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTttatt                                                                                                                  <  4:204524/122‑1 (MQ=255)
                                                             ctttCTACTGAAAGATAAACCCATTTAAATTCATAAAATTACAATCACAAACCGGTCATTCTCTTTGCGCATCTTTAATATTACATTACCGGAAATTTCTGCCCTGTCACACATATTTTATTACGCTCCTTCCGCTACTGTGTCCACtt                                                                                       >  4:513571/2‑149 (MQ=255)
                                                                     tGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGt                                                                               <  4:73035/149‑1 (MQ=255)
                                                                     tGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGt                                                                               >  4:55926/1‑149 (MQ=255)
                                                                     tGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGt                                                                               >  4:109353/1‑149 (MQ=255)
                                                                      gAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGtt                                                                              <  3:505925/149‑1 (MQ=255)
                                                                         aGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCa                                                                           <  4:338001/149‑1 (MQ=255)
                                                                              aaaCCATTTAAATTCATAAAATTACATTTACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAg                                                                      >  4:203656/1‑149 (MQ=255)
                                                                               aaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAgg                                                                     >  3:260109/1‑149 (MQ=255)
                                                                                aCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACATGACAGGAAATTTCTTCCCTGTAACACACATTTTATTATGCTGCATTCGCTACTGCGTCCCCTTCCGGACGTTTCAGGAAGGc                                                                    >  4:45231/1‑149 (MQ=255)
                                                                                aCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGc                                                                    <  1:546082/149‑1 (MQ=255)
                                                                                   aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTcc                                                                                           <  4:248711/123‑1 (MQ=255)
                                                                                   aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTcc                                                                                           >  3:248711/1‑123 (MQ=255)
                                                                                   aTTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATa                                                                 >  2:430444/1‑149 (MQ=255)
                                                                                    tttAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAg                                                                <  4:6466/149‑1 (MQ=255)
                                                                                     ttAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAgg                                                               <  3:357603/149‑1 (MQ=255)
                                                                                         aTTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaa                                                           >  3:148388/1‑149 (MQ=255)
                                                                                          ttCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaaa                                                          <  1:219262/149‑1 (MQ=255)
                                                                                            cATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAcc                                                        >  3:62119/1‑149 (MQ=255)
                                                                                             aTAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAccc                                                       >  2:438811/1‑149 (MQ=255)
                                                                                                  aTTACATTCACAAACCGGTCTTTCTCTTTGCGCTCATCAAACGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGccac                                                  >  3:370236/1‑149 (MQ=255)
                                                                                                       aTTCACAACCCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCg                                             <  2:280244/149‑1 (MQ=255)
                                                                                                       aTTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCg                                             >  1:12859/1‑149 (MQ=255)
                                                                                                           acaAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAg                                                                      >  1:571690/1‑120 (MQ=255)
                                                                                                           acaAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAg                                                                      <  2:571690/120‑1 (MQ=255)
                                                                                                            caAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCg                                        >  1:81967/1‑149 (MQ=255)
                                                                                                              aaCCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGc                                      >  2:556891/1‑149 (MQ=255)
                                                                                                                   gTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCaataa                                 >  1:239597/1‑149 (MQ=255)
                                                                                                                     cATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAAtt                               <  3:63396/149‑1 (MQ=255)
                                                                                                                      aTTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTg                              >  3:186960/1‑149 (MQ=255)
                                                                                                                            tttGCGCTCATCAAATGTTACAGGGCAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCa                        <  2:153387/149‑1 (MQ=255)
                                                                                                                            tttGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCa                        <  2:275898/149‑1 (MQ=255)
                                                                                                                            tttGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCa                        >  1:575239/1‑149 (MQ=255)
                                                                                                                               gcgcTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGt                     <  4:348523/149‑1 (MQ=255)
                                                                                                                                 gcTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTaa                   >  2:278380/1‑149 (MQ=255)
                                                                                                                                    catcaAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAAccc                <  1:4321/149‑1 (MQ=255)
                                                                                                                                           tGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCg         <  2:224020/149‑1 (MQ=255)
                                                                                                                                             ttACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGc       <  1:360185/149‑1 (MQ=255)
                                                                                                                                             ttACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGc       >  4:504131/1‑149 (MQ=255)
                                                                                                                                             ttACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGc       >  2:360185/1‑149 (MQ=255)
                                                                                                                                             ttACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGc       <  3:504131/149‑1 (MQ=255)
                                                                                                                                                  ggacaggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCGTAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAAt  >  4:320742/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAAT  >  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: