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

GCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTC  >  NC_000913/2259404‑2259697
                                                                                                                                                 |                                                                                                                                                    
gCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaat                                                                                                                                                   >  1:589223/1‑149 (MQ=255)
 cATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaata                                                                                                                                                  <  3:135442/149‑1 (MQ=255)
 cATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaata                                                                                                                                                  >  4:135442/1‑149 (MQ=255)
    gTTGTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAAAGCGATGAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATc                                                                                                                                               <  3:591373/149‑1 (MQ=255)
     ttttACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCt                                                                                                                                              >  1:295258/1‑149 (MQ=255)
         aCGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGccc                                                                                                                                          <  2:172182/149‑1 (MQ=255)
           gtaaagtaaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATGAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCtt                                                                                                                                        >  4:546529/1‑149 (MQ=255)
            taaagtaaaTTACCGTCAGAGAGGTAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGATTATTTTCTGCCCttt                                                                                                                                       <  1:497628/149‑2 (MQ=255)
            taaagtaaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTg                                                                                                                                       >  4:130655/1‑149 (MQ=255)
             aaagtaaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGa                                                                                                                                      <  1:165223/149‑1 (MQ=255)
              aagtaaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAg                                                                                                                                     <  1:689393/149‑1 (MQ=255)
                gtaaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCTTTGAGGa                                                                                                                                   >  3:8143/1‑149 (MQ=255)
                gtaaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGa                                                                                                                                      <  4:699580/146‑1 (MQ=255)
                gtaaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGa                                                                                                                                      >  3:699580/1‑146 (MQ=255)
                 taaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAGATAATCTGCCCtt                                                                                                                                        <  2:588686/143‑1 (MQ=255)
                 taaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCtt                                                                                                                                        >  1:588686/1‑143 (MQ=255)
                   aaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAAt                                                                                                                                >  2:731481/1‑149 (MQ=255)
                     ttACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCTGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGACCAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATaa                                                                                                                              <  3:546529/149‑1 (MQ=255)
                     ttACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATaa                                                                                                                              <  3:365186/149‑1 (MQ=255)
                     ttACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATaa                                                                                                                              >  4:540591/1‑149 (MQ=255)
                          gCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGa                                                                                                                         <  3:346698/149‑1 (MQ=255)
                              gagagGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAg                                                                                                                                     <  1:222957/133‑1 (MQ=255)
                              gagagGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAg                                                                                                                                     >  2:222957/1‑133 (MQ=255)
                               agagGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGtt                                                                                                                    <  2:589223/149‑1 (MQ=255)
                                gagGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGttt                                                                                                                   <  1:35498/149‑1 (MQ=255)
                                  ggCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtga                                                                                                                 <  1:274688/149‑1 (MQ=255)
                                   tctaaaTGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGttt                                                                                                                   <  4:691309/143‑1 (MQ=255)
                                   gCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACGATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGttt                                                                                                                   >  3:691309/1‑146 (MQ=255)
                                   gCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCt                                                                                                                           <  1:64328/138‑1 (MQ=255)
                                   gCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCt                                                                                                                           >  2:64328/1‑138 (MQ=255)
                                    ctaaaTGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtgatg                                                                                                               >  4:693443/3‑149 (MQ=255)
                                       aaTGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCTTTGAGGAAATAACCTGAGTGTTTGATGAAc                                                                                                            <  4:8143/149‑1 (MQ=255)
                                         tGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGa                                                                                                                         >  1:2469/1‑134 (MQ=255)
                                         tGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGa                                                                                                                         <  2:2469/134‑1 (MQ=255)
                                            gCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTg                                                                                                       >  2:562036/1‑149 (MQ=255)
                                            gCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTg                                                                                                       >  2:779261/1‑149 (MQ=255)
                                               gATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAg                                                                                                    <  2:295258/149‑1 (MQ=255)
                                                aTTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGt                                                                                                   >  2:731582/1‑149 (MQ=255)
                                                 tttGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGtt                                                                                                  >  3:418121/1‑149 (MQ=255)
                                                      gCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTAc                                                                                             >  2:736537/1‑149 (MQ=255)
                                                              acaAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGAt                                                                                     >  3:351761/1‑149 (MQ=255)
                                                                aaaaTGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATaa                                                                                   >  1:539199/1‑149 (MQ=255)
                                                                  aaTGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATaaaa                                                                                 >  4:328189/1‑149 (MQ=255)
                                                                  aaTGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATaaaa                                                                                 <  1:770968/149‑1 (MQ=255)
                                                                  aaTGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATaaaa                                                                                 <  4:309658/149‑1 (MQ=255)
                                                                      cGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCAt                                                                             <  1:779261/149‑1 (MQ=255)
                                                                       gATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCAtt                                                                            <  2:539199/149‑1 (MQ=255)
                                                                       gATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCAtt                                                                            <  1:731481/149‑1 (MQ=255)
                                                                        aTTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTa                                                                                              >  3:616491/1‑130 (MQ=255)
                                                                        aTTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTa                                                                                              <  4:616491/130‑1 (MQ=255)
                                                                         ttAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTa                                                                          >  1:325161/1‑149 (MQ=255)
                                                                            aaTGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAAt                                                                       <  3:707407/149‑1 (MQ=255)
                                                                                ccTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCAt                                                                   >  3:441351/1‑149 (MQ=255)
                                                                                ccTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCAt                                                                   <  4:134739/149‑1 (MQ=255)
                                                                                 cTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATa                                                                  >  1:548725/1‑149 (MQ=255)
                                                                                  tGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATaa                                                                 >  4:468612/1‑149 (MQ=255)
                                                                                         cAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACa                                                          <  3:693443/149‑1 (MQ=255)
                                                                                           aaaGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACAtt                                                        <  4:418121/149‑1 (MQ=255)
                                                                                                gAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACaa                                                   >  2:197401/1‑149 (MQ=255)
                                                                                                        tgtgACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCa                                           <  3:130655/149‑1 (MQ=255)
                                                                                                            aCGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATtct                                       <  3:540591/149‑1 (MQ=255)
                                                                                                              gggATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTctct                                     <  4:351761/149‑1 (MQ=255)
                                                                                                                   gCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCAtt                                         >  4:745765/1‑140 (MQ=255)
                                                                                                                   gCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCAtt                                         <  3:745765/140‑1 (MQ=255)
                                                                                                                     acagaacagaACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTc                                                       <  4:372722/124‑1 (MQ=255)
                                                                                                                     acagaacagaACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTc                                                       >  3:372722/1‑124 (MQ=255)
                                                                                                                     acagaacagaACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCt                              >  4:584063/1‑149 (MQ=255)
                                                                                                                     acagaacagaACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCt                              >  2:629936/1‑149 (MQ=255)
                                                                                                                     acagaacagaACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCt                              >  2:749854/1‑149 (MQ=255)
                                                                                                                        gaacagaacaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCtcat                           >  2:46898/1‑149 (MQ=255)
                                                                                                                         aacagaacaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCtcatc                          >  1:424515/1‑149 (MQ=255)
                                                                                                                             gaacaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAAt                      <  2:548725/149‑1 (MQ=255)
                                                                                                                              aacaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTgcgc                               <  4:12020/139‑1 (MQ=255)
                                                                                                                              aacaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTgcgc                               >  3:12020/1‑139 (MQ=255)
                                                                                                                                 aaTAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGa                                                                                                                         <  3:326899/46‑1 (MQ=255)
                                                                                                                                 aaTAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGa                                                                                                                         >  4:326899/1‑46 (MQ=255)
                                                                                                                                 aaTAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTa                  <  3:328189/149‑1 (MQ=255)
                                                                                                                                 aaTAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTa                  >  4:159832/1‑149 (MQ=255)
                                                                                                                                 aaTAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTa                  >  2:156830/1‑149 (MQ=255)
                                                                                                                                  aTAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTac                 <  1:197401/149‑1 (MQ=255)
                                                                                                                                     ctcaaGGGCAGATTATTTTCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGATCGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTc                                                       <  4:673427/105‑1 (MQ=38)
                                                                                                                                     ctcaaGGGCAGATTATTTTCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTc                                                       >  3:673427/4‑108 (MQ=39)
                                                                                                                                     ctcaaGGGCAGATTATTTTCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTaca                <  4:315323/144‑1 (MQ=255)
                                                                                                                                     ctcaaGGGCAGATTATTTTCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTaca                >  3:315323/4‑147 (MQ=255)
                                                                                                                                      aaaaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCaaa                       >  1:476610/1‑139 (MQ=255)
                                                                                                                                      aaaaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCaaa                       <  2:476610/139‑1 (MQ=255)
                                                                                                                                      aaaaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacagga             >  4:281288/1‑149 (MQ=255)
                                                                                                                                      aaaaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacagga             <  3:584063/149‑1 (MQ=255)
                                                                                                                                      aaaaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacagga             <  3:281288/149‑1 (MQ=255)
                                                                                                                                        aaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacaggaca           <  1:562036/149‑1 (MQ=255)
                                                                                                                                        aaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacaggaca           <  1:46898/149‑1 (MQ=255)
                                                                                                                                           ggCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTAcaggacagga        >  2:429939/1‑149 (MQ=255)
                                                                                                                                                 aaaTAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTc  <  1:731582/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTC  >  NC_000913/2259404‑2259697

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