Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,681,969 C→G 6.3% Q259H (CAG→CAC) ‡ 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,9690CG 6.3% 351.5 / 12.9 111Q259H (CAG→CAC) ‡yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base C (52/52);  new base G (4/3);  total (56/55)
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.89e-01

Only 100 of 114 total aligned reads displayed.

CCGTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGT  >  NC_000913/2681820‑2682118
                                                                                                                                                     |                                                                                                                                                      
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:936870/150‑1 (MQ=255)
        gTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGAttt                                                                                                                                                >  1:578074/1‑150 (MQ=255)
           ggTCAGACTCTTATCCCACGCCTGGCCAAATTCAGAGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAAAATGGTGGAGGGCAACGTGCAGCATACCGCCATCTGCATCGTGCCACAACGCGCCGACAAAATCGTACTGCGGTTTTTc                                                                                                                                            >  1:277024/1‑150 (MQ=255)
           ggTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTc                                                                                                                                            <  2:671440/150‑1 (MQ=255)
            gTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGCAGTACGAttttgtc                                                                                                                                            <  2:1254153/150‑4 (MQ=255)
             tCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCgg                                                                                                                                          <  1:1350881/150‑1 (MQ=255)
             tCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCgg                                                                                                                                          <  2:1306540/150‑1 (MQ=255)
              cAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGCAGTACGAttttgtcgg                                                                                                                                          <  1:1454517/150‑6 (MQ=255)
               aGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGAAAAATCGCAGTACGATTTTGTCGGcg                                                                                                                                        <  1:34184/150‑2 (MQ=255)
                    ctTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAGCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCatata                                                                                                                                   >  2:1418867/1‑150 (MQ=255)
                     ttATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCatatat                                                                                                                                  <  1:1657588/150‑1 (MQ=255)
                        tCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCACACCGCCATCTTCATCGTGACCCCACTCGCCGACAAAAGCGAACTGCGATTTTTCGGCCAGATATGAg                                                                                                                               >  1:299484/1‑150 (MQ=255)
                         cccACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAgg                                                                                                                              >  2:722292/1‑150 (MQ=255)
                         cccACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAgg                                                                                                                              >  2:1153478/1‑150 (MQ=255)
                          ccACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGt                                                                                                                             <  1:1559698/150‑1 (MQ=255)
                            aCGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAATATCGTACTGCGATTTTTCGGCCATATATGAGGTTg                                                                                                                           <  1:1462517/150‑1 (MQ=255)
                            aCGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTg                                                                                                                           <  2:879160/150‑1 (MQ=255)
                            aCGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTg                                                                                                                           <  1:424206/150‑1 (MQ=255)
                             cGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTgg                                                                                                                          <  2:551620/150‑1 (MQ=255)
                                cTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAAc                                                                                                                       <  1:880368/150‑1 (MQ=255)
                                   gCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAgg                                                                                                                    >  2:529778/1‑150 (MQ=255)
                                   gCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAgg                                                                                                                    >  2:1294268/1‑150 (MQ=255)
                                    ccAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGc                                                                                                                   <  2:1537618/150‑1 (MQ=255)
                                          tCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACAttt                                                                                                             >  1:545051/1‑150 (MQ=255)
                                           cACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACAtttt                                                                                                            <  2:1218901/150‑1 (MQ=255)
                                           cACTGTGTCACCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACAtttt                                                                                                            >  1:1153011/1‑150 (MQ=255)
                                                tgtCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAg                                                                                                       <  1:1193786/150‑1 (MQ=255)
                                                tgtCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGAGGTTTGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAg                                                                                                       <  2:1079760/150‑1 (MQ=255)
                                                 gtCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGc                                                                                                      <  1:395954/150‑1 (MQ=255)
                                                   ccccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCga                                                                                                    >  1:1459333/1‑150 (MQ=255)
                                                    cccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACCCCACGCGCCGACAAAATCGTACTGCGATTTTGCGGCCATATATGAGGTTGGAACAGGCACATTTTTAGAGCgag                                                                                                   >  2:1101041/1‑150 (MQ=255)
                                                    cccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCgag                                                                                                   <  1:1590695/150‑1 (MQ=255)
                                                    cccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCgag                                                                                                   <  2:1487953/150‑1 (MQ=255)
                                                    cccAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCgag                                                                                                   >  1:624282/1‑150 (MQ=255)
                                                    cccAACTCCACTGTGTTTTAGCTGGGGCAATATGGTGTTTGGCAACGTGCAGCAAACCGACATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCgag                                                                                                   <  2:258796/150‑1 (MQ=255)
                                                        aCTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAAt                                                                                               <  1:168904/150‑1 (MQ=255)
                                                         cTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATg                                                                                              <  2:687240/150‑1 (MQ=255)
                                                         cTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATg                                                                                              <  1:594275/150‑1 (MQ=255)
                                                          tCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGt                                                                                             >  2:284518/1‑150 (MQ=255)
                                                           ccACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCATCAAACCGCAATCTTCATCTTGAAACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAATTGTc                                                                                            <  1:240458/150‑1 (MQ=255)
                                                            cACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTcc                                                                                           >  1:234033/1‑150 (MQ=255)
                                                              cTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCAc                                                                                         <  2:1725752/150‑1 (MQ=255)
                                                               tGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACt                                                                                        >  2:510607/1‑150 (MQ=255)
                                                                    tttACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGc                                                                                   >  1:423831/1‑150 (MQ=255)
                                                                        ccTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGag                                                                               >  1:374399/1‑150 (MQ=255)
                                                                        acTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGag                                                                               <  1:850847/149‑1 (MQ=255)
                                                                            gcgcAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCAAGCGGagtagt                                                                           >  2:1387608/1‑150 (MQ=255)
                                                                              gcAATCCGGTGGCGGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGAACTGCGAGTTATCGGCCATATATGAGGATGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTcc                                                                         <  2:88944/150‑1 (MQ=255)
                                                                                aaTATGGTTGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCAc                                                                       <  2:1180878/150‑1 (MQ=255)
                                                                                aaTATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCAc                                                                       <  2:1586406/150‑1 (MQ=255)
                                                                                  tatGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACtt                                                                     <  1:666540/150‑1 (MQ=255)
                                                                                    tggtggTTGGCAACGTGCAGCATACCGGCATCTTCTTCGTGACACGACGCGGCGAGAAAATCGTCCTGCGATTTTTCGGCCTTATATGAGGTTGGATCTGGCTCATTTTTTTAGCGTGATATGTGCACTCCATCGGAGTAGTCCACTTTg                                                                   <  2:1730296/150‑1 (MQ=255)
                                                                                    tggtggTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTg                                                                   >  1:965469/1‑150 (MQ=255)
                                                                                     ggtggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTGGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgt                                                                  <  1:1600424/150‑1 (MQ=255)
                                                                                      gtggtTGGCAACGTGCAGCAAACCGCAATCTTCATCTTGACACCACGCGCCGACAAAATCGTCCTGAGATTTTTAGGCCATATATCAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCCCTCCAGCTGAGTAGTCCACTTAgta                                                                 >  2:963582/1‑150 (MQ=255)
                                                                                       tggtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgtag                                                                <  1:1060263/150‑1 (MQ=255)
                                                                                         gtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTTACACCACGCGCCGTCAAAATCGTACTGCTATTTTTCGGCCATATATGAGGTTGTAACAGGCACATTTTTATAGCGAGAATTGTCCACTCCAGCGGAGTAGTCCACTTTgtagta                                                              <  2:1199705/150‑1 (MQ=255)
                                                                                         gtTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTgtagta                                                              <  2:1121899/150‑1 (MQ=255)
                                                                                              cAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGc                                                         >  1:1557559/1‑150 (MQ=255)
                                                                                                     cagcaAACCGCCATCTTCATCGTGACACCACGCGCCGACTAAATCGTACTGCGCTTTTTCGGCCATTTATGCGGTTGGAACATGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGc                                                  <  1:1597768/150‑1 (MQ=255)
                                                                                                     cagcaAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGc                                                  <  2:1771440/150‑1 (MQ=255)
                                                                                                          aaCCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTTGCGATgg                                             <  2:1153011/150‑1 (MQ=255)
                                                                                                            ccGCCATCTTCATCGTGACGCCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgggg                                           >  1:411881/1‑150 (MQ=255)
                                                                                                            ccGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATgggg                                           >  1:1324866/1‑150 (MQ=255)
                                                                                                             cGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTATGTAGTAAGTGCCGGTAGCGATGGGGa                                          <  1:528243/150‑1 (MQ=255)
                                                                                                              gCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGaa                                         >  1:1684708/1‑150 (MQ=255)
                                                                                                                cATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAg                                       >  1:1734022/1‑150 (MQ=255)
                                                                                                                   cTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCgg                                    >  2:921617/1‑150 (MQ=255)
                                                                                                                    ttCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGga                                   <  1:1711511/150‑1 (MQ=255)
                                                                                                                     tCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATCTATGAGGTTGGAACAGGCAGATTTTTATAGCGAGTAATTTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGgag                                  <  1:343379/150‑1 (MQ=255)
                                                                                                                     tCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGGTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGGGGAGTAGTCCACTTTGTAGTAAGTGCCGGTATCGATGGGGAAAGCGgag                                  <  1:1331346/150‑1 (MQ=255)
                                                                                                                     tCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTGCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCTCTCCAGCGGAGTAGTCCACTTTGTAGTAATTGCCGGTAGCGATGGGGAAAGCGgag                                  <  1:348537/150‑1 (MQ=255)
                                                                                                                      cATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGTCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGTAAAGCGgaga                                 >  2:1697426/1‑150 (MQ=255)
                                                                                                                      cATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGgaga                                 >  2:725836/1‑150 (MQ=255)
                                                                                                                      cATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGgaga                                 >  1:1169518/1‑150 (MQ=255)
                                                                                                                        tCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTTCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACg                               >  1:181371/1‑150 (MQ=255)
                                                                                                                        tCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTTGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACg                               >  1:75061/1‑150 (MQ=255)
                                                                                                                         cGTGACACCACGCGCCGACACAATCGTACTGCTAGTTTTCGGCCATATCTGAGGTTGGCACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACgg                              <  1:1615174/150‑1 (MQ=255)
                                                                                                                          gTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGc                             >  1:173200/1‑150 (MQ=255)
                                                                                                                           tGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGcc                            >  1:415139/1‑150 (MQ=255)
                                                                                                                              caccacGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGt                         >  2:1066313/1‑150 (MQ=255)
                                                                                                                                  acGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGc                     <  1:1216938/150‑1 (MQ=255)
                                                                                                                                   cgcgCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACATGCACATTTTTATAGCGAGATATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTATGTTCCGGTAGCGATGGGGAATGCTGAGACGGCGCGTTTGcc                    <  2:952069/150‑1 (MQ=255)
                                                                                                                                       gccgaaaaaTCGCAGTACGATTTTGTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTg                 >  1:719740/6‑150 (MQ=255)
                                                                                                                                       gccgaaaaaTCGCAGTACGATTTTGTCGGCCAGATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCAACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGACGGTAGCGATGGGGAAAGAGGAGAAGGCCCGTTTGCCGTg                 >  2:1024385/6‑150 (MQ=255)
                                                                                                                                        ccgaaaaaTCGCAGTACGATTTTGTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGa                >  2:1504916/5‑150 (MQ=255)
                                                                                                                                        cGACAAAATCGTACTGAGATTTATCGGCCATATATGATGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCATACGAGTAGTAAACTTTGAAGTAAGTGACGGTATCCATGGGGAAAACGGATACTGCCCGTTTGcctttaa               >  2:328623/1‑145 (MQ=255)
                                                                                                                                           cAAAATCGTACTGCGATTTTTCGGCCATATATGAGTTTGGCACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTGTGTTGTAAGTGCCGGTAGCGATGGGGAAAGCGGCGACGGCCCGTTTGACGTGATCaa            >  2:1137474/1‑150 (MQ=255)
                                                                                                                                           cAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATTTCCAGTCCAGCGGAGTAGGCAACTTTGTTGTAAGTGCCGGTAGGGAGGGGGAAGGCGGACCCGGCCCGTTTCCCGTGATaaa            >  2:794615/1‑150 (MQ=255)
                                                                                                                                           cAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCaa            >  2:1120665/1‑150 (MQ=255)
                                                                                                                                           cAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCaa            >  1:916788/1‑150 (MQ=255)
                                                                                                                                             aactCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAac          <  1:1393755/147‑1 (MQ=255)
                                                                                                                                             aaaTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAac          >  2:19453/1‑150 (MQ=255)
                                                                                                                                               aTCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGCCGGCCCGTTTGCCGTGATCAAacac        >  1:1254612/1‑150 (MQ=255)
                                                                                                                                                     cTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGt  >  1:1024754/1‑150 (MQ=255)
                                                                                                                                                     cTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGt  <  1:1729410/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCGTTATTGTCGGTCAGACTCTTATCCCACGCCTGGCCAAATTCACTGTGTCCCCAACTCCACTGTTTTTTACCTGGCGCAATATGGTGGTTGGCAACGTGCAGCAAACCGCCATCTTCATCGTGACACCACGCGCCGACAAAATCGTACTGCGATTTTTCGGCCATATATGAGGTTGGAACAGGCACATTTTTATAGCGAGAAATGTCCACTCCAGCGGAGTAGTCCACTTTGTAGTAAGTGCCGGTAGCGATGGGGAAAGCGGAGACGGCCCGTTTGCCGTGATCAAACACCGCCGT  >  NC_000913/2681820‑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: