Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,710,703 G→A 6.4% intergenic (‑228/+96) eptB ← / ← yhjX Kdo2‑lipid A phosphoethanolamine 7''‑transferase/putative pyruvate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,710,7030GA 6.4% 248.5 / 4.5 78intergenic (‑228/+96)eptB/yhjXKdo2‑lipid A phosphoethanolamine 7''‑transferase/putative pyruvate transporter
Reads supporting (aligned to +/- strand):  ref base G (42/31);  new base A (3/2);  total (45/33)
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.67e-01

GATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAATCGT  >  NC_000913/3710560‑3710852
                                                                                                                                               |                                                                                                                                                     
gATAAGTAAAATCTACACGCTCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAAtt                                                                                                                                                 <  1:381854/150‑1 (MQ=255)
gATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAAtt                                                                                                                                                 <  2:439134/150‑1 (MQ=255)
  tAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTgg                                                                                                                                               >  2:972694/1‑150 (MQ=255)
  tAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTgg                                                                                                                                               <  1:760989/150‑1 (MQ=255)
    aGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGcc                                                                                                                                             <  1:31324/150‑1 (MQ=255)
         aaTCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGAAAAAAGGGCCAATTAACTGGccctttt                                                                                                                                        <  2:805081/150‑5 (MQ=255)
            cTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAAtttt                                                                                                                                     <  1:972694/150‑1 (MQ=255)
             aacacGCGCGAATAAGTAGAACACAGCAAACGGAAAAAGTGTGGGGTAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTa                                                                                                                                    <  2:410153/149‑1 (MQ=255)
                 cgcgcgAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAg                                                                                                                                >  2:381088/1‑150 (MQ=255)
                 cgcgcgAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAg                                                                                                                                <  1:342366/150‑1 (MQ=255)
                         tAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGtt                                                                                                                        <  1:547701/150‑1 (MQ=255)
                          aaGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGttt                                                                                                                       >  2:44872/1‑150 (MQ=255)
                           aGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAATGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTta                                                                                                                      <  1:259266/150‑1 (MQ=255)
                           aGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTta                                                                                                                      <  1:1247912/150‑1 (MQ=255)
                           aGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTta                                                                                                                      >  2:129835/1‑150 (MQ=255)
                            gTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTtat                                                                                                                     <  1:540665/150‑1 (MQ=255)
                            gTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTCCATTAATCAAAAAATTACAAAAATCCAATTTTTAAAAAAAAATTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTtat                                                                                                                     >  2:994246/1‑150 (MQ=255)
                             tAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTtata                                                                                                                    <  1:962725/150‑1 (MQ=255)
                              aGACCACAGCAAAGTGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTGATTGTTtattc                                                                                                                   <  1:179756/150‑3 (MQ=255)
                              aGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATAc                                                                                                                   <  1:348012/150‑1 (MQ=255)
                               gACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTgg                                                                                                                                               >  2:250486/1‑121 (MQ=255)
                               gACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTgg                                                                                                                                               <  1:250486/121‑1 (MQ=255)
                               caCCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACa                                                                                                                  <  1:728674/149‑1 (MQ=255)
                                       cAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGGAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTg                                                                                                          <  2:660394/150‑1 (MQ=255)
                                        aaaGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAATCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTAGa                                                                                                         >  2:837129/1‑150 (MQ=255)
                                             gAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGt                                                                                                    >  2:1182647/1‑150 (MQ=255)
                                              aaaaaGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGtt                                                                                                   <  2:813393/150‑1 (MQ=255)
                                                  aGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGAt                                                                                               >  2:164946/1‑150 (MQ=255)
                                                      tgGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGcac                                                                                           <  1:1182647/150‑1 (MQ=255)
                                                          gtaaatttgtgCTTCATGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACgg                                                                                       >  1:1352970/8‑150 (MQ=255)
                                                          gAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACgg                                                                                       >  2:347664/1‑150 (MQ=255)
                                                                 tgtgCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGa                                                                                >  2:93075/1‑150 (MQ=255)
                                                                   tgCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAg                                                                              <  1:1243933/150‑1 (MQ=255)
                                                                   tgCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAg                                                                              <  1:93075/150‑1 (MQ=255)
                                                                   tgCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAg                                                                              >  2:1276775/1‑150 (MQ=255)
                                                                    gCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGtt                                                                                                           >  2:1350143/1‑120 (MQ=255)
                                                                    gCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGtt                                                                                                           <  1:1350143/120‑1 (MQ=255)
                                                                    gCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGt                                                                                                            <  2:1390806/119‑1 (MQ=255)
                                                                    gCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGt                                                                                                            >  1:1390806/1‑119 (MQ=255)
                                                                    gCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGc                                                                             >  2:74569/1‑150 (MQ=255)
                                                                        gAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGAACACGGAAGATGAAGCCCAg                                                                         >  1:1241474/1‑150 (MQ=255)
                                                                         aaGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGt                                                                        <  1:381088/150‑1 (MQ=255)
                                                                              caAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTc                                                                   >  1:1056831/1‑150 (MQ=255)
                                                                              caAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTc                                                                   >  2:858177/1‑150 (MQ=255)
                                                                                      cAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGa                                                           >  2:434213/1‑150 (MQ=255)
                                                                                      cAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGa                                                           >  2:929988/1‑150 (MQ=255)
                                                                                      cAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGa                                                           >  2:1371838/1‑150 (MQ=255)
                                                                                       aaaaaTGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAg                                                          <  1:978115/150‑1 (MQ=255)
                                                                                                aGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAgg                                                 >  2:387867/1‑150 (MQ=255)
                                                                                                  tttGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGa                                               <  1:129835/150‑1 (MQ=255)
                                                                                                  tttGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGa                                               >  2:717131/1‑150 (MQ=255)
                                                                                                     gAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGcc                                            >  2:374887/1‑150 (MQ=255)
                                                                                                       aaaGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCAt                                          >  2:817951/1‑150 (MQ=255)
                                                                                                        aaGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATg                                         >  2:402255/1‑150 (MQ=255)
                                                                                                         agagACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATgc                                        >  2:1294606/1‑150 (MQ=255)
                                                                                                          gagaCTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATgcg                                       >  1:226239/1‑150 (MQ=255)
                                                                                                           agaCTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATgcgc                                      <  2:1056831/150‑1 (MQ=255)
                                                                                                            gaCTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGcc                                     >  2:1086128/1‑150 (MQ=255)
                                                                                                                  aGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGc                               >  2:103052/1‑150 (MQ=255)
                                                                                                                        aCAAAAATCAGTAAAAAGGGCCAGTTAATAGCCGCAAAATTTACCAGACATTGTTTATACATTTCGTAGATGCCTTTGATGCACACGGAACCTAAAGCCCAGTAGCAATCGGCTGAGCATTAAAGGGAGAAATGCGCCTCAAGCAACAtt                         >  2:344875/1‑150 (MQ=255)
                                                                                                                           aaaaGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTc                      >  1:1293744/1‑150 (MQ=255)
                                                                                                                            aaaGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCt                     <  1:164946/150‑1 (MQ=255)
                                                                                                                             aaGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTg                    >  1:34678/1‑150 (MQ=255)
                                                                                                                             aaGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTg                    <  1:44872/150‑1 (MQ=255)
                                                                                                                              aGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGc                   <  1:1371838/150‑1 (MQ=255)
                                                                                                                                cAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGctc                 >  2:489182/1‑150 (MQ=255)
                                                                                                                                  gTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTg               >  2:342465/1‑150 (MQ=255)
                                                                                                                                  gTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTg               >  1:970749/1‑150 (MQ=255)
                                                                                                                                    aatttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGc             >  2:893458/6‑150 (MQ=255)
                                                                                                                                    aaaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGc             >  2:506313/1‑150 (MQ=255)
                                                                                                                                    aaaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGc             <  1:402255/150‑1 (MQ=255)
                                                                                                                                    aaaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGc             >  2:1239556/1‑150 (MQ=255)
                                                                                                                                    aaaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGc             >  2:856385/1‑150 (MQ=255)
                                                                                                                                     atttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCt            >  1:1004065/5‑150 (MQ=255)
                                                                                                                                     aaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCt            >  2:1028578/1‑150 (MQ=255)
                                                                                                                                      aaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTg           >  1:1043813/1‑150 (MQ=255)
                                                                                                                                        tgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACAATTCGTTGATGCGATTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTAACCCAACATTTACTGCTCTGGCTGAc         >  2:266462/2‑150 (MQ=255)
                                                                                                                                              aGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAAtcg   <  1:506313/150‑1 (MQ=255)
                                                                                                                                               gTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAAtcgt  <  1:1086128/150‑1 (MQ=255)
                                                                                                                                               gTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAAtcgt  <  1:103052/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAATCGT  >  NC_000913/3710560‑3710852

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