Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,681,971 G→A 6.2% Q259* (CAG→TAG) ‡ yphG ← DUF4380 domain‑containing TPR repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,681,9710GA 6.2% 347.3 / 8.1 113Q259* (CAG→TAG) ‡yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (52/53);  new base A (4/3);  total (57/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.33e-01

Only 100 of 113 total aligned reads displayed.

ATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGT  >  NC_000913/2681825‑2682118
                                                                                                                                                  |                                                                                                                                                    
aTTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGa                                                                                                                                                   >  1:1586406/1‑150 (MQ=255)
  tGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGAtt                                                                                                                                                 >  2:1216938/1‑150 (MQ=255)
   gTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCCCTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAATTGGTACTGCGAttt                                                                                                                                                >  2:828277/1‑150 (MQ=255)
   gTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGAttt                                                                                                                                                >  1:578074/1‑150 (MQ=255)
      ggTCAGACTCTTATCCCACGCCTGGCCAAATTCAGAGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAAAATGGTGGAGGGCAACGTGCAGCATACCGCCATCTGCATCGTGCCACAACGCGCCGACAAAATCGTACTGCGGTTTT‑Tc                                                                                                                                            >  1:277024/1‑150 (MQ=255)
      ggTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑Tc                                                                                                                                            <  2:671440/150‑1 (MQ=255)
      ggTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCCTACTGCGATTTT‑Tc                                                                                                                                            >  2:1331346/1‑150 (MQ=255)
       gTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGCAGTACGAttttgtc                                                                                                                                            <  2:1254153/150‑4 (MQ=255)
        tCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCgg                                                                                                                                          <  1:1350881/150‑1 (MQ=255)
        tCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCgg                                                                                                                                          <  2:1306540/150‑1 (MQ=255)
         cAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGCAGTACGAttttgtcgg                                                                                                                                          <  1:1454517/150‑6 (MQ=255)
          aGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGA‑AAAATCGCAGTACGATTTTGTCGGcg                                                                                                                                        <  1:34184/150‑2 (MQ=255)
           taCTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCa                                                                                                                                       <  1:155249/149‑1 (MQ=255)
               ctTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAGCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCatata                                                                                                                                   >  2:1418867/1‑150 (MQ=255)
                ttATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCatatat                                                                                                                                  <  1:1657588/150‑1 (MQ=255)
                   tCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCACACCGCCATCTTCATCGTGACCCCACTCGCCGACAAAAGCGAACTGCGATTTT‑TCGGCCAGATATGAg                                                                                                                               >  1:299484/1‑150 (MQ=255)
                    cccACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAgg                                                                                                                              >  2:1153478/1‑150 (MQ=255)
                    cccACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAgg                                                                                                                              >  2:1711511/1‑150 (MQ=255)
                    cccACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAgg                                                                                                                              >  2:722292/1‑150 (MQ=255)
                       aCGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAATATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTg                                                                                                                           <  1:1462517/150‑1 (MQ=255)
                       aCGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTg                                                                                                                           <  1:424206/150‑1 (MQ=255)
                       aCGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTg                                                                                                                           <  2:879160/150‑1 (MQ=255)
                        cGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTgg                                                                                                                          <  2:551620/150‑1 (MQ=255)
                           cTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAAc                                                                                                                       <  1:880368/150‑1 (MQ=255)
                              gCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCCCGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACggg                                                                                                                    >  1:1634054/1‑150 (MQ=255)
                              gCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAgg                                                                                                                    >  2:529778/1‑150 (MQ=255)
                               ccAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGc                                                                                                                   <  2:1537618/150‑1 (MQ=255)
                                     tCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACAttt                                                                                                             >  1:545051/1‑150 (MQ=255)
                                      cACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACAtttt                                                                                                            <  2:1218901/150‑1 (MQ=255)
                                      cACTGTGTCACCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACAtttt                                                                                                            >  1:1153011/1‑150 (MQ=255)
                                           tgtCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAg                                                                                                       <  1:1193786/150‑1 (MQ=255)
                                           tgtCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGAGGTTTGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAg                                                                                                       <  2:1079760/150‑1 (MQ=255)
                                            gtCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGc                                                                                                      <  1:395954/150‑1 (MQ=255)
                                              ccccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCga                                                                                                    >  1:1459333/1‑150 (MQ=255)
                                               cccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACCCCACGCGCCGACAAAATCGTACTGCGATTTT‑GCGGCCATATATGAGGTTGGAACAGGCACATTTTTAGAGCgag                                                                                                   >  2:1101041/1‑150 (MQ=255)
                                               cccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCgag                                                                                                   <  1:1590695/150‑1 (MQ=255)
                                               cccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCgag                                                                                                   <  2:1487953/150‑1 (MQ=255)
                                               cccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCgag                                                                                                   >  1:1273903/1‑150 (MQ=255)
                                               cccAACTCCACTGTGTTTTAGCTGGGGCAATATGGTGTTTGGCAACGTGCAGCAAACCGACATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCgag                                                                                                   <  2:258796/150‑1 (MQ=255)
                                                    cTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATg                                                                                              <  2:687240/150‑1 (MQ=255)
                                                    cTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATg                                                                                              <  1:594275/150‑1 (MQ=255)
                                                     tCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGt                                                                                             >  2:284518/1‑150 (MQ=255)
                                                       cACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTcc                                                                                           >  1:234033/1‑150 (MQ=255)
                                                        aCTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCa                                                                                          >  1:1744299/1‑150 (MQ=255)
                                                         cTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCAc                                                                                         <  2:1725752/150‑1 (MQ=255)
                                                          ttttttttACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCTCCATCTTCATCGTGACACCACGCGCAGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACATGCACATTTTTATAGCGAGAAATGTCCACt                                                                                        <  1:348026/148‑1 (MQ=255)
                                                          tGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACt                                                                                        >  2:510607/1‑150 (MQ=255)
                                                               tttACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGc                                                                                   >  1:423831/1‑150 (MQ=255)
                                                                 tACCTCGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCgg                                                                                 <  2:1174186/150‑1 (MQ=255)
                                                                   ccTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGag                                                                               >  1:374399/1‑150 (MQ=255)
                                                                   acTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGag                                                                               <  1:850847/149‑1 (MQ=255)
                                                                       gcgcAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCAAGCGGagtagt                                                                           >  2:1387608/1‑150 (MQ=255)
                                                                         gcAATCCGGTGGCGGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGAACTGCGAGTTA‑TCGGCCATATATGAGGATGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTcc                                                                         <  2:88944/150‑1 (MQ=255)
                                                                           aaTATGGTTGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCAc                                                                       <  2:1180878/150‑1 (MQ=255)
                                                                           aaTATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCCAGAAATGTCCACTCCAGCGGAGTAGTCCAc                                                                       <  1:722292/150‑1 (MQ=255)
                                                                             tatGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACtt                                                                     <  1:666540/150‑1 (MQ=255)
                                                                               tggtggTTGGCAACGTGCAGCATACCGGCATCTTCTTCGTGACACGACGCGGCGAGAAAATCGTCCTGCGATTTT‑TCGGCCTTATATGAGGTTGGATCTGGCTCATTTTTTTAGCGTGATATGTGCACTCCATCGGAGTAGTCCACTTTg                                                                   <  2:1730296/150‑1 (MQ=255)
                                                                               tggtggTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTg                                                                   >  1:965469/1‑150 (MQ=255)
                                                                                ggtggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTTCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATTTCCACTCCAGCGGAGTAGTCCACTTTgt                                                                  >  2:435602/1‑150 (MQ=255)
                                                                                ggtggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTGGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgt                                                                  <  1:1600424/150‑1 (MQ=255)
                                                                                 gtggtTGGCAACGTGCAGCAAACCGCAATCTTCATCTTGACACCACGCGCCGACAAAATCGTCCTGAGATTTT‑TAGGCCATATATCAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCCCTCCAGCTGAGTAGTCCACTTAgta                                                                 >  2:963582/1‑150 (MQ=255)
                                                                                  tggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgtag                                                                <  1:1060263/150‑1 (MQ=255)
                                                                                    gtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTTACACCACGCGCCGTCAAAATCGTACTGCTATTTT‑TCGGCCATATATGAGGTTGTAACAGGCACATTTTTATAGCGAGAATTGTCCACTCCAGCGGAGTAGTCCACTTTgtagta                                                              <  2:1199705/150‑1 (MQ=255)
                                                                                    gtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgtagta                                                              <  2:1121899/150‑1 (MQ=255)
                                                                                         cAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGc                                                         >  1:1557559/1‑150 (MQ=255)
                                                                                                cagcaAACCGCCATCTTCATCGTGACACCACGCGCCGACTAAATCGTACTGCGCTTTT‑TCGGCCATTTATGCGGTTGGAACATGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGc                                                  <  1:1597768/150‑1 (MQ=255)
                                                                                                cagcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGc                                                  <  2:1771440/150‑1 (MQ=255)
                                                                                                     aaCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTTGCGATgg                                             <  2:1153011/150‑1 (MQ=255)
                                                                                                        cGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTATGTAGTAAGTGCCGGTAGCGATGGGGa                                          <  1:528243/150‑1 (MQ=255)
                                                                                                         gCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGaa                                         >  1:1684708/1‑150 (MQ=255)
                                                                                                           cATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAg                                       >  1:1734022/1‑150 (MQ=255)
                                                                                                              cTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCgg                                    >  2:921617/1‑150 (MQ=255)
                                                                                                               ttCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGga                                   <  1:1711511/150‑1 (MQ=255)
                                                                                                                tCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATCTATGAGGTTGGAACAGGCAGATTTTTATAGCGAGTAATTTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGgag                                  <  1:343379/150‑1 (MQ=255)
                                                                                                                tCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGGTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGGGGAGTAGTCCACTTTGTAGTAAGTGCCGGTATCGATGGGGAAAGCGgag                                  <  1:1331346/150‑1 (MQ=255)
                                                                                                                tCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑GCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCTCTCCAGCGGAGTAGTCCACTTTGTAGTAATTGCCGGTAGCGATGGGGAAAGCGgag                                  <  1:348537/150‑1 (MQ=255)
                                                                                                                 cATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGTCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGTAAAGCGgaga                                 >  2:1697426/1‑150 (MQ=255)
                                                                                                                 cATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGgaga                                 >  1:1169518/1‑150 (MQ=255)
                                                                                                                 cATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGgaga                                 >  2:725836/1‑150 (MQ=255)
                                                                                                                   tCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTTCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACg                               >  1:181371/1‑150 (MQ=255)
                                                                                                                   tCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTTGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACg                               >  1:75061/1‑150 (MQ=255)
                                                                                                                    cGTGACACCACGCGCCGACACAATCGTACTGCTAGTTT‑TCGGCCATATCTGAGGTTGGCACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACgg                              <  1:1615174/150‑1 (MQ=255)
                                                                                                                     gTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGc                             >  1:173200/1‑150 (MQ=255)
                                                                                                                      tGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGcc                            >  1:415139/1‑150 (MQ=255)
                                                                                                                         caccacGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGt                         >  2:1066313/1‑150 (MQ=255)
                                                                                                                             acGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGc                     <  1:1216938/150‑1 (MQ=255)
                                                                                                                              cgcgCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACATGCACATTTTTATAGCGAGATATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTATGTTCCGGTAGCGATGGGGAATGCTGAGACGGCGCGTTTGcc                    <  2:952069/150‑1 (MQ=255)
                                                                                                                                  gccgaaaaaTCGCAGTACGATTTTGTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTg                 >  1:503653/6‑150 (MQ=255)
                                                                                                                                  gccgaaaaaTCGCAGTACGATTTTGTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTg                 >  1:719740/6‑150 (MQ=255)
                                                                                                                                  gccgaaaaaTCGCAGTACGATTTTGTCGGCCAGATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCAACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGACGGTAGCGATGGGGAAAGAGGAGAAGGCCCGTTTGCCGTg                 >  2:1024385/6‑150 (MQ=255)
                                                                                                                                   ccgaaaaaTCGCAGTACGATTTTGTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGa                >  2:1504916/5‑150 (MQ=255)
                                                                                                                                   cGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGAt               <  2:1069322/150‑1 (MQ=255)
                                                                                                                                      cAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATTTCCAGTCCAGCGGAGTAGGCAACTTTGTTGTAAGTGCCGGTAGGGAGGGGGAAGGCGGACCCGGCCCGTTTCCCGTGATaaa            >  2:794615/1‑150 (MQ=255)
                                                                                                                                      cAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCaa            >  2:1120665/1‑150 (MQ=255)
                                                                                                                                      cAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCaa            >  1:916788/1‑150 (MQ=255)
                                                                                                                                        aactCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAac          <  1:1393755/147‑1 (MQ=255)
                                                                                                                                        aaaTCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAac          >  2:19453/1‑150 (MQ=255)
                                                                                                                                          aTCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGCCGGCCCGTTTGCCGTGATCAAacac        >  1:1254612/1‑150 (MQ=255)
                                                                                                                                                cTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGt  <  1:1729410/150‑1 (MQ=255)
                                                                                                                                                cTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGt  >  1:1024754/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGT  >  NC_000913/2681825‑2682118

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