Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 950,648 C→T 5.8% V145I (GTT→ATT) ‡ pflA ← pyruvate formate‑lyase activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913950,6480CT 5.8% 442.1 / 10.2 138V145I (GTT→ATT) ‡pflApyruvate formate‑lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base C (58/72);  new base T (2/6);  total (60/78)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

Only 100 of 143 total aligned reads displayed.

AACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCTTTTTTGC  >  NC_000913/950499‑950797
                                                                                                                                                     |                                                                                                                                                     
aaCGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAAc                                                                                                                                                       <  1:349931/150‑1 (MQ=255)
       cccATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAAcc                                                                                                                                                <  1:373005/150‑1 (MQ=255)
           tatCACGGGTAAATTCACCGAGGCGATGCTCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGAGTGCGGTGGGTGGAGTTTCCAACCaccg                                                                                                                                            <  2:1800446/150‑4 (MQ=255)
           tagcACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTTCTGTGGTTGGAGTTTCCAACCaccg                                                                                                                                            <  2:1791564/147‑4 (MQ=255)
            atCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGTACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAGTTTCCAACCaccgc                                                                                                                                           <  2:1897883/150‑5 (MQ=255)
              cACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGAtttt                                                                                                                                         >  1:317780/1‑150 (MQ=255)
              cACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGAtttt                                                                                                                                         >  2:1013774/1‑150 (MQ=255)
                ggggTAAATACACCGAGGCGATGCTCTGAATCGTCATCGTCAGAACAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTtgg                                                                                                                                       <  2:1594773/149‑1 (MQ=255)
                cGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTtgg                                                                                                                                       >  2:44926/1‑150 (MQ=255)
                 gggTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTtggt                                                                                                                                      >  2:187478/1‑150 (MQ=255)
                  ggTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTtggtg                                                                                                                                     >  1:898252/1‑150 (MQ=255)
                   gTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTtggtgg                                                                                                                                    <  1:875369/150‑1 (MQ=255)
                     aaaTTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGAt                                                                                                                                  >  1:1064693/1‑150 (MQ=255)
                      aaTTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGAtc                                                                                                                                 >  2:2140132/1‑150 (MQ=255)
                       aTTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGAtct                                                                                                                                <  1:1093239/150‑1 (MQ=255)
                          cACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCtcgt                                                                                                                             <  2:383736/150‑1 (MQ=255)
                               aGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTc                                                                                                                        <  1:259730/150‑1 (MQ=255)
                                   taTGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCt                                                                                                                    <  2:384888/149‑1 (MQ=255)
                                       cgcTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGttt                                                                                                                <  1:712602/150‑1 (MQ=255)
                                        gcTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTg                                                                                                               <  1:777394/150‑1 (MQ=255)
                                        gcTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCACCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTg                                                                                                               >  1:2081412/1‑150 (MQ=255)
                                         cTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTga                                                                                                              <  2:638871/150‑1 (MQ=255)
                                         cTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTga                                                                                                              <  2:1238450/150‑1 (MQ=255)
                                            aaTCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGAt                                                                                                           <  2:1114449/150‑1 (MQ=255)
                                              tcgtcatcgtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCg                                                                                                         <  1:1169097/150‑1 (MQ=255)
                                                gtcatcgtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCaa                                                                                                                                                  <  1:808947/107‑1 (MQ=255)
                                                gtcatcgtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAg                                                                                                       >  1:196550/1‑150 (MQ=255)
                                                     cgtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTa                                                                                                  >  1:1393468/1‑150 (MQ=255)
                                                        cagaccagacTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCa                                                                                               <  2:1813213/150‑1 (MQ=255)
                                                        caGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCa                                                                                               <  2:531282/150‑1 (MQ=255)
                                                            ccAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTc                                                                                           >  1:2235742/1‑150 (MQ=255)
                                                            ccAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGGGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTc                                                                                           <  1:2455980/150‑1 (MQ=255)
                                                               gCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGCAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGt                                                                                        >  2:1343062/1‑150 (MQ=255)
                                                               gCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGt                                                                                        <  1:1939571/150‑1 (MQ=255)
                                                                    ggACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACtt                                                                                   >  2:1550182/1‑150 (MQ=255)
                                                                    ggACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGGGGATCTCGTCGGTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACtt                                                                                   <  1:764956/150‑1 (MQ=255)
                                                                     gACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTc                                                                                  <  2:842687/150‑1 (MQ=255)
                                                                      acaacaACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTACAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGATATCGAGCATTACCAGGTCGGTTACTTcc                                                                                 <  2:1596581/150‑1 (MQ=255)
                                                                       caacaaCGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCca                                                                                >  2:1127863/1‑150 (MQ=255)
                                                                       caacaaCGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCca                                                                                >  1:1975345/1‑150 (MQ=255)
                                                                         acaacGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCcagc                                                                              >  2:1458020/1‑150 (MQ=255)
                                                                         acaacGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCcagc                                                                              <  1:583994/150‑1 (MQ=255)
                                                                         acaacGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCcagc                                                                              >  1:1255084/1‑150 (MQ=255)
                                                                            acGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGt                                                                           <  2:1220315/150‑1 (MQ=255)
                                                                            acGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGt                                                                           <  2:898252/150‑1 (MQ=255)
                                                                             cGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGtt                                                                          <  2:115707/150‑1 (MQ=255)
                                                                               tAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTtca                                                                        <  1:44926/150‑1 (MQ=255)
                                                                                aGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTtcat                                                                       <  1:2021237/150‑1 (MQ=255)
                                                                                aGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTtcat                                                                       <  1:398114/150‑1 (MQ=255)
                                                                                 gCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTtcatc                                                                      <  1:1561997/150‑1 (MQ=255)
                                                                                  cGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGGTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTtcatca                                                                     >  1:719495/1‑150 (MQ=255)
                                                                                        cacaCCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCAcc                                                               >  1:2463396/1‑150 (MQ=255)
                                                                                            ccTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGAt                                                           >  2:1226496/1‑150 (MQ=255)
                                                                                             cTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATc                                                          >  2:750352/1‑150 (MQ=255)
                                                                                              ttCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCg                                                         >  2:1807466/1‑150 (MQ=255)
                                                                                              ttCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCg                                                         <  1:13642/150‑1 (MQ=255)
                                                                                              ttCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCg                                                         >  1:1638373/1‑150 (MQ=255)
                                                                                              ttCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCg                                                         >  1:162171/1‑150 (MQ=255)
                                                                                               tCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGt                                                        <  1:333621/150‑1 (MQ=255)
                                                                                               tCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGt                                                        >  1:110791/1‑150 (MQ=255)
                                                                                                cacaTTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTa                                                       <  2:2291239/150‑1 (MQ=255)
                                                                                                  caTTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGAGGTCGTTGCAAACTCCAACCAGATTTTGGTGGATCTCGACGTTCATCTGTTTGAGAAAGAGCATTCCCAGGTCGGTTACTTGAAGCACTTCATCAATCACCGGATCGTacc                                                     >  2:850871/1‑148 (MQ=255)
                                                                                                   aTTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAacg                                                    <  1:2433793/150‑1 (MQ=255)
                                                                                                    tttttGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAacga                                                   >  2:2268315/1‑150 (MQ=255)
                                                                                                       gagTTCGCTAAATATCTGGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAAcgacga                                                >  2:287948/3‑150 (MQ=255)
                                                                                                       gagTTCGCTAAATATCTGGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGATCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAAcgacga                                                >  2:1465092/3‑150 (MQ=255)
                                                                                                           tCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACaa                                            <  1:856469/150‑1 (MQ=255)
                                                                                                             gCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACaaaa                                          <  1:1541611/150‑1 (MQ=255)
                                                                                                              ccAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAAc                                         <  1:645691/150‑1 (MQ=255)
                                                                                                               cAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAAcc                                        <  2:1454580/150‑1 (MQ=255)
                                                                                                                  atatTTAGCGAACTCCAGCGTGCGGTGGCTGGAAACTCCAACCAGATTTTGCTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGtt                                     <  2:451477/150‑1 (MQ=255)
                                                                                                                   tatTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTg                                    <  1:999568/150‑1 (MQ=255)
                                                                                                                   tatTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTg                                    <  1:180596/150‑1 (MQ=255)
                                                                                                                   tatTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTCCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTg                                    >  2:295337/1‑150 (MQ=255)
                                                                                                                    atTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTgg                                   <  2:759702/150‑1 (MQ=255)
                                                                                                                    atTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTgg                                   <  2:2394797/150‑1 (MQ=255)
                                                                                                                      ttAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGgtg                                 >  2:1044028/1‑150 (MQ=255)
                                                                                                                      ttAGCGAACGCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGCCGAACAAAACCGTTGgtg                                 >  2:2436613/1‑150 (MQ=255)
                                                                                                                       tAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTgg                                   >  1:2337801/1‑147 (MQ=255)
                                                                                                                       tAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGgtgt                                >  2:1834370/1‑150 (MQ=255)
                                                                                                                        aGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTc                               <  1:2432415/150‑1 (MQ=255)
                                                                                                                        aGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTc                               <  2:2081412/150‑1 (MQ=255)
                                                                                                                         gCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTcc                              >  2:682550/1‑150 (MQ=255)
                                                                                                                         gCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTcc                              >  1:1263755/1‑150 (MQ=255)
                                                                                                                         gCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTcc                              <  2:891019/150‑1 (MQ=255)
                                                                                                                           gAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAg                            <  2:514540/150‑1 (MQ=255)
                                                                                                                           gAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAg                            >  1:796080/1‑150 (MQ=255)
                                                                                                                             aCTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGTTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGAc                          <  2:1413237/150‑1 (MQ=255)
                                                                                                                             aCTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGAc                          >  1:992219/1‑150 (MQ=255)
                                                                                                                                ccAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGGCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAgg                       <  1:1716725/150‑1 (MQ=255)
                                                                                                                                        cGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATg               <  1:1458020/150‑1 (MQ=255)
                                                                                                                                        cGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATg               <  1:1333083/150‑1 (MQ=255)
                                                                                                                                          gtggtTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGcc             <  2:2234113/150‑1 (MQ=255)
                                                                                                                                          ctggtTGGAGTTTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGcc             >  1:399566/2‑150 (MQ=255)
                                                                                                                                              ttGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCcttct         >  1:1927253/1‑150 (MQ=255)
                                                                                                                                                 gAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAgg                       >  1:1189804/1‑133 (MQ=255)
                                                                                                                                                    aCTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCTTTTTTg   <  2:317780/150‑1 (MQ=255)
                                                                                                                                                    aCTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCTTTTTTg   <  2:1549979/150‑1 (MQ=255)
                                                                                                                                                     cTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCTTTTTTGc  >  2:587443/1‑150 (MQ=255)
                                                                                                                                                     cTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCTTTTTTGc  >  1:20940/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCTTTTTTGC  >  NC_000913/950499‑950797

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: