Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,259,549 A→T 7.3% intergenic (‑253/+170) 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,5490AT 7.3% 291.8 / 10.1 82intergenic (‑253/+170)psuK/fruApseudouridine kinase/fused fructose‑specific PTS enzymes: IIBcomponent/IIC components
Reads supporting (aligned to +/- strand):  ref base A (36/40);  new base T (3/3);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAA  >  NC_000913/2259403‑2259675
                                                                                                                                                  |                                                                                                                              
tGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGaaaa                                                                                                                              >  2:1113974/1‑149 (MQ=255)
  cATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaata                                                                                                                            >  3:500152/1‑149 (MQ=255)
   aTGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaataa                                                                                                                           <  3:819377/149‑1 (MQ=255)
    tGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaataat                                                                                                                          >  2:1441575/1‑149 (MQ=255)
          aCGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGATTATTTTCTGccc                                                                                                                    <  2:1026154/149‑1 (MQ=255)
            gtaaagtaaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaataat                                                                                                                          <  3:120171/141‑1 (MQ=255)
            gtaaagtaaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaataat                                                                                                                          >  4:120171/1‑141 (MQ=255)
               aagtaaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAg                                                                                                               >  4:1172198/1‑149 (MQ=255)
                 gtaaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGa                                                                                                             >  2:1396773/1‑149 (MQ=255)
                  taaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAGTGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGaa                                                                                                            >  4:613513/1‑149 (MQ=255)
                    aaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAAt                                                                                                          <  1:1396773/149‑1 (MQ=255)
                      ttACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGc                                                                                                                      >  2:574863/1‑135 (MQ=255)
                      ttACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGc                                                                                                                      <  1:574863/135‑1 (MQ=255)
                       tACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAAc                                                                                                       <  4:500152/149‑1 (MQ=255)
                        aCCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAAcc                                                                                                      <  4:377466/149‑1 (MQ=255)
                        aCCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAAcc                                                                                                      >  4:197326/1‑149 (MQ=255)
                          cGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAAc                                                                                                       <  1:1089619/146‑1 (MQ=255)
                          cGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAAc                                                                                                       >  2:1089619/1‑146 (MQ=255)
                                 gagGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGttt                                                                                             <  1:1441575/149‑1 (MQ=255)
                                    gCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGATTATTTTCTGCCCt                                                                                                                   >  1:576460/1‑124 (MQ=255)
                                    gCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGATTATTTTCTGCCCt                                                                                                                   <  2:576460/124‑1 (MQ=255)
                                    gCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtgat                                                                                          >  4:308163/1‑149 (MQ=255)
                                    gCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtgat                                                                                          <  3:1172198/149‑1 (MQ=255)
                                    gCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtgat                                                                                          <  2:1303047/149‑1 (MQ=255)
                                    gCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtgat                                                                                          >  1:1303047/1‑149 (MQ=255)
                                            ggCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGttt                                                                                  <  3:308163/149‑1 (MQ=255)
                                            ggCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGttt                                                                                  <  2:355136/149‑1 (MQ=255)
                                             gCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtga                                                                                           <  1:811553/139‑1 (MQ=255)
                                             gCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtga                                                                                           >  2:811553/1‑139 (MQ=255)
                                                  tttGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAg                                                                                                               >  2:409527/1‑114 (MQ=255)
                                                  tttGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAg                                                                                                               <  1:409527/114‑1 (MQ=255)
                                                         aGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTg                                                                     >  1:194326/1‑149 (MQ=255)
                                                         aGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTg                                                                     <  3:197326/149‑1 (MQ=255)
                                                           gCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaataat                                                                                                                          <  3:949141/94‑1 (MQ=255)
                                                           gCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaataat                                                                                                                          >  4:342808/1‑94 (MQ=255)
                                                           gCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaataat                                                                                                                          <  3:342808/94‑1 (MQ=255)
                                                           gCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaataat                                                                                                                          >  4:949141/1‑94 (MQ=255)
                                                           gCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGaa                                                                   >  1:156863/1‑149 (MQ=255)
                                                            cTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGaaa                                                                  <  2:194326/149‑1 (MQ=255)
                                                             tCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGATTATTTTCTGccc                                                                                                                    >  1:46731/1‑98 (MQ=255)
                                                             tCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGATTATTTTCTGccc                                                                                                                    <  2:46731/98‑1 (MQ=255)
                                                             tCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGATTATTTTCTGCCCtttt                                                                                                                >  3:1392255/1‑100 (MQ=39)
                                                                   aaTGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATaaaa                                                           <  2:156863/149‑1 (MQ=255)
                                                                     tGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGaaa                                                                  >  1:431158/1‑140 (MQ=255)
                                                                     tGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGaaa                                                                  <  2:431158/140‑1 (MQ=255)
                                                                      gCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCa                                                        <  1:1113974/149‑1 (MQ=255)
                                                                        gATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTg                                                                                                                 <  1:1392569/90‑1 (MQ=255)
                                                                        gATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTg                                                                                                                 >  2:1392569/1‑90 (MQ=255)
                                                                        gATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtgat                                                                                          <  1:1281818/113‑1 (MQ=255)
                                                                        gATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtgat                                                                                          >  2:1281818/1‑113 (MQ=255)
                                                                        gATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCAtt                                                      >  4:620074/1‑149 (MQ=255)
                                                                         aTTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGaaa                                                                  >  2:775677/1‑136 (MQ=255)
                                                                         aTTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGaaa                                                                  <  1:775677/136‑1 (MQ=255)
                                                                         aTTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCAt                                                       <  4:683747/147‑1 (MQ=255)
                                                                         aTTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCAt                                                       >  3:683747/1‑147 (MQ=255)
                                                                                gCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCa                                              >  2:877426/1‑149 (MQ=255)
                                                                                  cTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATaa                                                                                                        >  4:1163402/1‑89 (MQ=255)
                                                                                  cTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATaa                                                                                                        <  3:1163402/89‑1 (MQ=255)
                                                                                   tGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATaa                                           >  4:827558/1‑149 (MQ=255)
                                                                                    gAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGaaa                                                                  >  2:552745/1‑125 (MQ=255)
                                                                                    gAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGaaa                                                                  <  1:552745/125‑1 (MQ=255)
                                                                                    gAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATaaa                                          >  2:1429783/1‑149 (MQ=255)
                                                                                        tACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTa                                      <  3:620074/149‑1 (MQ=255)
                                                                                          cAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaataat                                                                                                                          >  3:126915/1‑63 (MQ=255)
                                                                                          cAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaataat                                                                                                                          <  4:126915/63‑1 (MQ=255)
                                                                                          cAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAgt                                                                                                 <  2:435128/88‑1 (MQ=255)
                                                                                          cAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAgt                                                                                                 >  1:435128/1‑88 (MQ=255)
                                                                                          cAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACa                                    <  3:827558/149‑1 (MQ=255)
                                                                                           aaaaGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGtt                                                                                              <  3:479110/90‑1 (MQ=255)
                                                                                           aaaaGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGtt                                                                                              >  4:479110/1‑90 (MQ=255)
                                                                                             aaGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATc                                                                                                                         <  4:766193/61‑1 (MQ=255)
                                                                                             aaGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATc                                                                                                                         >  3:766193/1‑61 (MQ=255)
                                                                                                  aGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACaaa                            <  1:877426/149‑1 (MQ=255)
                                                                                                       tttGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGt                       <  1:1429783/149‑1 (MQ=255)
                                                                                                               gggATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAGTCTGCCCTTGAGGAAATAACCTGAGTGTTtgat                                                                                          <  3:573811/74‑1 (MQ=255)
                                                                                                               gggATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAGTCTGCCCTTGAGGAAATAACCTGAGTGTTtgat                                                                                          >  4:573811/1‑74 (MQ=255)
                                                                                                                      acagaacagaACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTc                                                                          <  3:1414727/83‑1 (MQ=255)
                                                                                                                      acagaacagaACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTc                                                                          >  4:1414727/1‑83 (MQ=255)
                                                                                                                      acagaacagaACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCt        >  4:300560/1‑149 (MQ=255)
                                                                                                                       cagaacagaaCAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATaaaa                                         <  1:1296561/115‑1 (MQ=255)
                                                                                                                       cagaacagaaCAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATaaaa                                         >  2:1296561/1‑115 (MQ=255)
                                                                                                                        agaacagaacAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCtca      <  3:300560/149‑1 (MQ=255)
                                                                                                                            cagaacaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCaa  <  3:613513/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                              
TGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAA  >  NC_000913/2259403‑2259675

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