Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 950,646 A→G 7.0% V145V (GTT→GTC) ‡ pflA ← pyruvate formate‑lyase activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913950,6460AG 7.0% 251.2 / 7.5 85V145V (GTT→GTC) ‡pflApyruvate formate‑lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base A (32/47);  new base G (3/3);  total (35/50)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCTTTT  >  NC_000913/950505‑950793
                                                                                                                                             |                                                                                                                                                   
gCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAAc                                                                                                                                             <  1:1045843/150‑1 (MQ=255)
gCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAAc                                                                                                                                             <  1:1051189/150‑1 (MQ=255)
gCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAAc                                                                                                                                             >  1:1091467/1‑150 (MQ=255)
    atatCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAGTTTCCAACCacc                                                                                                                                         <  1:1115865/150‑3 (MQ=255)
    atatCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAGTTTCCAACCacc                                                                                                                                         <  1:322777/150‑3 (MQ=255)
     tatCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAGTTTCCAACCaccg                                                                                                                                        <  1:608236/150‑4 (MQ=255)
      atCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGAtt                                                                                                                                       >  1:540222/1‑150 (MQ=255)
        cACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATAGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGTGGTGGTTGGAAACTCCAACCAGAtttt                                                                                                                                     >  2:1145945/1‑150 (MQ=255)
         aCGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTtg                                                                                                                                    >  2:1452157/1‑150 (MQ=255)
          cGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTtgg                                                                                                                                   <  2:1500780/150‑1 (MQ=255)
          cGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTtgg                                                                                                                                   <  1:594607/150‑1 (MQ=255)
             gTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTtggtgg                                                                                                                                <  2:26950/150‑1 (MQ=255)
               aaaTTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGAt                                                                                                                              <  1:1114416/150‑1 (MQ=255)
                aaTTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGAtc                                                                                                                             >  2:1122062/1‑150 (MQ=255)
                 aTTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGAtct                                                                                                                            <  2:104815/150‑1 (MQ=255)
                 aTTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGAtct                                                                                                                            <  1:1444411/150‑1 (MQ=255)
                    cACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCtcgt                                                                                                                         >  1:1399894/1‑150 (MQ=255)
                    cACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCtcgt                                                                                                                         >  1:362849/1‑150 (MQ=255)
                     aCCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCtcgtc                                                                                                                        <  2:920719/150‑1 (MQ=255)
                       cGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTcgtcgt                                                                                                                      <  1:612090/150‑1 (MQ=255)
                          ggCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCa                                                                                                                   >  1:817655/1‑150 (MQ=255)
                            cGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATc                                                                                                                 <  1:1070067/150‑1 (MQ=255)
                            cGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATc                                                                                                                 <  1:541656/150‑1 (MQ=255)
                                gcgcTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGtt                                                                                                             >  2:714423/1‑150 (MQ=255)
                                 cgcTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGttt                                                                                                            <  2:953647/150‑1 (MQ=255)
                                  gcTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTg                                                                                                           <  2:45226/150‑1 (MQ=255)
                                  gcTGAATCGTCATCGTCAGACCAGCCGGGGACAACAACCTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTg                                                                                                           <  2:801180/150‑1 (MQ=255)
                                    tGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTgag                                                                                                         <  1:887033/150‑1 (MQ=255)
                                     gAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTtgg                                                                                                                                   <  1:59622/123‑1 (MQ=255)
                                     gAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTtgg                                                                                                                                   >  2:59622/1‑123 (MQ=255)
                                       atcgtcatcgtcAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATc                                                                                                      >  2:354271/1‑150 (MQ=255)
                                        tcgtcatcgtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCg                                                                                                     >  1:335270/1‑150 (MQ=255)
                                            catcgtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCa                                                                                                 >  1:591594/1‑150 (MQ=255)
                                            catcgtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCa                                                                                                 <  2:540222/150‑1 (MQ=255)
                                               cgtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTa                                                                                              <  1:1242336/150‑1 (MQ=255)
                                               cgtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTa                                                                                              >  1:1170410/1‑150 (MQ=255)
                                                  caGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCa                                                                                           <  2:1170410/150‑1 (MQ=255)
                                                   agcccAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCCGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCAATACCAg                                                                                          >  1:1126920/4‑150 (MQ=255)
                                                    gACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAgg                                                                                         >  2:425121/1‑150 (MQ=255)
                                                    gACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAgg                                                                                         >  2:387840/1‑150 (MQ=255)
                                                         gCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGt                                                                                    <  2:898942/150‑1 (MQ=255)
                                                            tGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTAcc                                                                                            >  2:117052/1‑139 (MQ=255)
                                                            tGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTAcc                                                                                            <  1:117052/139‑1 (MQ=255)
                                                                acaacaACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTcc                                                                             <  1:692931/150‑1 (MQ=255)
                                                                   acaacGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCcagc                                                                          <  2:111818/150‑1 (MQ=255)
                                                                    caacGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCcagca                                                                         <  1:1082187/150‑1 (MQ=255)
                                                                    caacGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCcagca                                                                         >  1:1300585/1‑150 (MQ=255)
                                                                        gTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCa                                                                                                 >  2:528129/1‑122 (MQ=255)
                                                                        gTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCa                                                                                                 <  1:528129/122‑1 (MQ=255)
                                                                        gTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTtc                                                                     >  1:970827/1‑150 (MQ=255)
                                                                         tAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTtca                                                                    >  2:628158/1‑150 (MQ=255)
                                                                         tAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTtca                                                                    >  2:1385062/1‑150 (MQ=255)
                                                                           gCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTtcatc                                                                  <  2:1399894/150‑1 (MQ=255)
                                                                             ggATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCatcaa                                                                >  1:202683/1‑150 (MQ=255)
                                                                                tCCCCACCTTCACCTTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCatcaatca                                                             >  2:857953/1‑150 (MQ=255)
                                                                                       cTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATc                                                      <  1:1085082/150‑1 (MQ=255)
                                                                                          cacaTTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTa                                                   <  2:961081/150‑1 (MQ=255)
                                                                                           acaTTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTaa                                                  <  1:692660/150‑1 (MQ=255)
                                                                                            caTTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAac                                                 <  1:109095/150‑1 (MQ=255)
                                                                                            caTTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAac                                                 >  2:512558/1‑150 (MQ=255)
                                                                                                tttGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAacgacg                                             >  1:378044/1‑150 (MQ=255)
                                                                                                  tGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGaa                                           <  1:354271/150‑1 (MQ=255)
                                                                                                     tCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGGCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACaa                                        <  2:890996/150‑1 (MQ=255)
                                                                                                             tatTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTg                                <  2:817655/150‑1 (MQ=255)
                                                                                                              atTTAGCGAACTCCAGCGTGTGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACTAACAAAACCGTTgg                               <  1:1145945/150‑1 (MQ=255)
                                                                                                              atTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTgg                               <  1:425121/150‑1 (MQ=255)
                                                                                                                  aGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTc                           >  2:252644/1‑150 (MQ=255)
                                                                                                                    cGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCa                         <  2:202683/150‑1 (MQ=255)
                                                                                                                    cGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGGGGATCTCGTCGGTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCa                         <  2:1468929/150‑1 (MQ=255)
                                                                                                                     gAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAg                        <  1:387840/150‑1 (MQ=255)
                                                                                                                        cTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACa                     <  1:1047203/150‑1 (MQ=255)
                                                                                                                        cTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACa                     <  2:307143/150‑1 (MQ=255)
                                                                                                                          ccAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAgg                   <  1:714423/150‑1 (MQ=255)
                                                                                                                              ggTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGATATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATg               <  2:807705/149‑1 (MQ=255)
                                                                                                                              caaaatctggtTGGAGTTTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATg               >  2:757067/8‑150 (MQ=255)
                                                                                                                              cGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATg               <  2:1300585/150‑1 (MQ=255)
                                                                                                                                 gCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAAt            >  2:631393/1‑150 (MQ=255)
                                                                                                                                  cGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATg           <  1:252644/150‑1 (MQ=255)
                                                                                                                                   tctggtTGGAGTTTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGc          >  2:1524520/3‑150 (MQ=255)
                                                                                                                                   tctggtTGGAGTTTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTAATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCAAGACAGGTATGAATGc          >  2:1511334/3‑150 (MQ=255)
                                                                                                                                   ggtggtTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGc          >  2:1420969/1‑150 (MQ=255)
                                                                                                                                      ggtTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCcag                                                                           >  1:415853/1‑82 (MQ=255)
                                                                                                                                      ggtTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCcag                                                                           <  2:415853/82‑1 (MQ=255)
                                                                                                                                         tGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCcttctt    >  1:913339/1‑150 (MQ=255)
                                                                                                                                           gAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCtttt  <  2:335270/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCTTTT  >  NC_000913/950505‑950793

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