Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,681,967 T→C 6.4% Y260C (TAC→TGC)  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,9670TC 6.4% 326.4 / 9.4 112Y260C (TAC→TGC) yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base T (51/52);  new base C (4/3);  total (56/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 = 7.21e-01

Only 100 of 114 total aligned reads displayed.

GGCCGTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACAC  >  NC_000913/2681818‑2682112
                                                                                                                                                     |                                                                                                                                                  
ggCCGTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCACCTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGt                                                                                                                                                    <  2:216380/150‑1 (MQ=255)
ggCCGTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGt                                                                                                                                                    <  1:1116677/150‑1 (MQ=255)
  ccGTTATTGTCGGTCAGACTCTTATCCCACGCCTGTCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGGGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTAc                                                                                                                                                  <  1:1109901/150‑1 (MQ=255)
       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)
          gTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGAttt                                                                                                                                          <  1:936870/150‑1 (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)
                      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)
                            ccACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGt                                                                                                                       <  1:1559698/150‑1 (MQ=255)
                              aCGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAATATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTg                                                                                                                     <  1:1462517/150‑1 (MQ=255)
                              aCGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTg                                                                                                                     <  2:879160/150‑1 (MQ=255)
                              aCGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTg                                                                                                                     <  1:424206/150‑1 (MQ=255)
                               cGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTgg                                                                                                                    <  2:551620/150‑1 (MQ=255)
                                  cTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAAc                                                                                                                 <  1:880368/150‑1 (MQ=255)
                                     gCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAgg                                                                                                              >  2:1294268/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)
                                                   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:624282/1‑150 (MQ=255)
                                                      cccAACTCCACTGTGTTTTAGCTGGGGCAATATGGTGTTTGGCAACGTGCAGCAAACCGACATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCgag                                                                                             <  2:258796/150‑1 (MQ=255)
                                                          aCTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAAt                                                                                         <  1:168904/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)
                                                             ccACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCATCAAACCGCAATCTTCATCTTGAAACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAATTGTc                                                                                      <  1:240458/150‑1 (MQ=255)
                                                              cACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTcc                                                                                     >  1:234033/1‑150 (MQ=255)
                                                               aCTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCa                                                                                    >  1:1744299/1‑150 (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)
                                                                                gcAATCCGGTGGCGGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGAACTGCGAGTTA‑TCGGCCATATATGAGGATGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTcc                                                                   <  2:88944/150‑1 (MQ=255)
                                                                                  aaTATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCAc                                                                 <  2:1586406/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)
                                                                                       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)
                                                                                                              ccGCCATCTTCATCGTGACGCCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgggg                                     >  1:411881/1‑150 (MQ=255)
                                                                                                              ccGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgggg                                     >  1:1324866/1‑150 (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)
                                                                                                                        cATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGTCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGTAAAGCGgaga                           >  2:1697426/1‑150 (MQ=255)
                                                                                                                        cATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGgaga                           >  1:1169518/1‑150 (MQ=255)
                                                                                                                        cATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGgaga                           >  2:725836/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)
                                                                                                                                          cGACAAAATCGTACTGAGATTTA‑TCGGCCATATATGATGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCATACGAGTAGTAAACTTTGAAGTAAGTGACGGTATCCATGGGGAAAACGGATACTGCCCGTTTGcctttaa         >  2:328623/1‑145 (MQ=255)
                                                                                                                                             cAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGTTTGGCACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTGTGTTGTAAGTGCCGGTAGCGATGGGGAAAGCGGCGACGGCCCGTTTGACGTGATCaa      >  2:1137474/1‑150 (MQ=255)
                                                                                                                                             cAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATTTCCAGTCCAGCGGAGTAGGCAACTTTGTTGTAAGTGCCGGTAGGGAGGGGGAAGGCGGACCCGGCCCGTTTCCCGTGATaaa      >  2:794615/1‑150 (MQ=255)
                                                                                                                                             cAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCaa      >  1:916788/1‑150 (MQ=255)
                                                                                                                                             cAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCaa      >  2:1120665/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)
                                                                                                                                                     |                                                                                                                                                  
GGCCGTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTT‑TCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACAC  >  NC_000913/2681818‑2682112

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: ‑