Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,146,589 A→G 15.8% D111G (GAC→GGC)  pflC → putative [formate‑C‑acetyltransferase 2]‑activating enzyme; pyruvate formate‑lyase 1‑activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,146,5890AG15.8% 1799.4 / 273.3 711D111G (GAC→GGC) pflCputative [formate‑C‑acetyltransferase 2]‑activating enzyme; pyruvate formate‑lyase 1‑activating enzyme
Reads supporting (aligned to +/- strand):  ref base A (310/287);  new base G (57/56);  total (367/344)
Fisher's exact test for biased strand distribution p-value = 8.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-01

Only 100 of 720 total aligned reads displayed.

GGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACTGCCGGAGACGCACCAGCCAGC  >  NC_000913/4146441‑4146737
                                                                                                                                                    |                                                                                                                                                    
 gCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGAc                                                                                                                                                     <  3:5134393/149‑1 (MQ=255)
 gCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGAc                                                                                                                                                     <  1:4234339/149‑1 (MQ=255)
      aTCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACAtttt                                                                                                                                                >  2:1483895/1‑149 (MQ=255)
       tCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGGCATtttc                                                                                                                                               >  4:3881156/1‑149 (MQ=255)
       tCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACAtttt                                                                                                                                                >  2:2763387/1‑148 (MQ=255)
       tCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATtttc                                                                                                                                               >  2:559150/1‑149 (MQ=255)
       tCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATtttc                                                                                                                                               <  3:389569/149‑1 (MQ=255)
       tCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATtttc                                                                                                                                               <  3:65469/149‑1 (MQ=255)
          aGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATtttcttt                                                                                                                                            <  4:5657868/149‑1 (MQ=255)
                 aCGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGt                                                                                                                                     <  1:1483895/149‑1 (MQ=255)
                 aCGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGt                                                                                                                                     >  4:333746/1‑149 (MQ=255)
                   gCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTcc                                                                                                                                   <  3:5408254/149‑1 (MQ=255)
                         agGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCcggcggc                                                                                                                             <  2:5393850/149‑1 (MQ=255)
                           gCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCcggcggcgg                                                                                                                           >  1:5509017/1‑149 (MQ=255)
                           gCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCcggcggcgg                                                                                                                           >  2:4944488/1‑149 (MQ=255)
                           gCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCcggcggcgg                                                                                                                           >  1:2544972/1‑149 (MQ=255)
                            cGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCggcggcggc                                                                                                                          <  4:3579879/149‑1 (MQ=255)
                             gAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGGCATtttcttt                                                                                                                                            >  4:1177363/1‑130 (MQ=255)
                                 tGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCAc                                                                                                                     >  3:5735241/1‑149 (MQ=255)
                                    cTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCt                                                                                                                  <  2:3428065/149‑1 (MQ=255)
                                     tGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAATGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCtt                                                                                                                 >  4:4817651/1‑149 (MQ=255)
                                     tGCACAGTGCGAACTGTTTGCGAGATGCGGATAAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGAGGGCATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCtt                                                                                                                 >  2:4264922/1‑149 (MQ=255)
                                      gCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCttt                                                                                                                <  4:6832105/149‑1 (MQ=255)
                                       cACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTc                                                                                                               <  4:2791037/149‑1 (MQ=255)
                                        aCTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCt                                                                                                              >  1:3537737/1‑149 (MQ=255)
                                           gtgCGAAATGTATGCGTGATGCGGATGAATGCCCCACCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGTTCTGGAACGGAAAGTCATCAAATATGACATTTTCTTTCGCACGTCCAGCGGCGGAGTCACGCTTTATggc                                                                                                           <  2:4518930/149‑1 (MQ=255)
                                                  aTGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAg                                                                                                    <  2:2179558/149‑1 (MQ=255)
                                                   tGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAgt                                                                                                   <  3:1549318/149‑1 (MQ=255)
                                                   tGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAgt                                                                                                   <  3:3042300/149‑1 (MQ=255)
                                                     tttGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAgtgt                                                                                                 >  2:6753581/1‑149 (MQ=255)
                                                     tttGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAgtgt                                                                                                 >  3:6917484/1‑149 (MQ=255)
                                                      ttGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGtt                                                                                                <  2:2459745/149‑1 (MQ=255)
                                                             aTGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAg                                                                                         >  3:2112137/1‑149 (MQ=255)
                                                              tGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGGCATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAgg                                                                                        <  3:184802/149‑1 (MQ=255)
                                                              tGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAgg                                                                                        >  3:2541158/1‑149 (MQ=255)
                                                               gCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGc                                                                                       >  3:1299566/1‑149 (MQ=255)
                                                                    tGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGt                                                                                  <  2:5696026/149‑1 (MQ=255)
                                                                     gAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGtt                                                                                 >  4:3155066/1‑149 (MQ=255)
                                                                        tgacccTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGc                                                                              <  4:2599361/146‑1 (MQ=255)
                                                                        tGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGc                                                                              >  3:2244304/1‑149 (MQ=255)
                                                                         gCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGGCATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCt                                                                             >  3:6802111/1‑149 (MQ=255)
                                                                           cccTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGGCATTTTCTTTCGc                                                                                                                                         <  1:3368307/87‑1 (MQ=255)
                                                                           cccTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTAc                                                                           >  3:453747/1‑149 (MQ=255)
                                                                            ccTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGGCATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTAcc                                                                          >  3:3244179/1‑149 (MQ=255)
                                                                              tCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCg                                                                        >  2:3196981/1‑149 (MQ=255)
                                                                                 ggggCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGttt                                                                     <  2:7025747/149‑1 (MQ=255)
                                                                                 ggggCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGttt                                                                     <  4:2937477/149‑1 (MQ=255)
                                                                                  gggCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGtttt                                                                    <  1:1671863/149‑1 (MQ=255)
                                                                                    gCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGAAGGCGGAGTTTGCTACCCGtttttt                                                                  >  3:4669851/1‑149 (MQ=255)
                                                                                     cGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTa                                                                 <  4:4369256/149‑1 (MQ=255)
                                                                                      gTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGGCATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTAc                                                                >  4:41770/1‑149 (MQ=255)
                                                                                      gTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTAc                                                                <  4:743720/149‑1 (MQ=255)
                                                                                       tttGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACa                                                               >  2:5995374/1‑149 (MQ=255)
                                                                                       tttGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACa                                                               >  2:3591941/1‑149 (MQ=255)
                                                                                        ttGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAg                                                              <  3:3588317/149‑1 (MQ=255)
                                                                                         tGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGc                                                             <  1:4803032/149‑1 (MQ=255)
                                                                                         tGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGc                                                             <  1:4421076/149‑1 (MQ=255)
                                                                                          gAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCg                                                            >  2:1297119/1‑149 (MQ=255)
                                                                                          gAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTCTACAGCg                                                            >  2:3991011/1‑149 (MQ=255)
                                                                                            aCGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGAc                                                          <  3:988322/149‑1 (MQ=255)
                                                                                            aCGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGAc                                                          <  3:723645/149‑1 (MQ=255)
                                                                                               gATTGGCCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTcc                                                                                                                                   >  2:5153841/1‑73 (MQ=255)
                                                                                                aTTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATtttctt                                                                                                                                             <  4:2275407/62‑1 (MQ=255)
                                                                                                 ttGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCgg                                                     <  1:1809474/149‑1 (MQ=255)
                                                                                                     tCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCtgt                                                 <  1:7203032/149‑1 (MQ=255)
                                                                                                      agcgATATCAGCCATGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCtgtg                                                >  4:4614447/2‑149 (MQ=255)
                                                                                                       gcgATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTgg                                               <  4:1730843/149‑1 (MQ=255)
                                                                                                       gcgATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTgg                                               >  3:1119211/1‑149 (MQ=255)
                                                                                                        cgATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTggg                                              <  1:3918670/149‑1 (MQ=255)
                                                                                                          atatCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGgt                                            <  1:827211/149‑1 (MQ=255)
                                                                                                           tatCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGttt                                                                     >  4:6302473/1‑123 (MQ=255)
                                                                                                             tCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGgtgtg                                         <  3:4841061/149‑1 (MQ=255)
                                                                                                              cAGCCTTGACGCTCTGGAAAGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGgtgtgt                                        >  3:4183135/1‑149 (MQ=255)
                                                                                                               aGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTc                                       >  1:5698587/1‑149 (MQ=255)
                                                                                                               aGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTc                                       <  1:210962/149‑1 (MQ=255)
                                                                                                               aGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTc                                       >  1:6095593/1‑149 (MQ=255)
                                                                                                                gCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCa                                      >  2:950710/1‑149 (MQ=255)
                                                                                                                 ccTTGACGCTCTGAAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGACGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCAt                                     <  1:6695608/149‑1 (MQ=255)
                                                                                                                   ttGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATgc                                   <  2:476946/149‑1 (MQ=255)
                                                                                                                    tGACGCTCTGGAACGGGAAGTGATGAAAGATGGCATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATgcg                                  <  3:6062250/149‑1 (MQ=255)
                                                                                                                     gACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATgcgc                                 >  3:5892801/1‑149 (MQ=255)
                                                                                                                       cGCTCTGGAACGGGAAGTGATGAAAGATGGCATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCa                               <  2:7288773/149‑1 (MQ=255)
                                                                                                                             ggAACGGGAAGTGATGAAAGATGGCATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGaaa                         <  2:2998304/149‑1 (MQ=255)
                                                                                                                             ggAACGGGAAGTGATGAAAGATGGCATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGaaa                         <  2:6470830/149‑1 (MQ=255)
                                                                                                                             ggAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGaaa                         <  4:1932948/149‑1 (MQ=255)
                                                                                                                              gAACGGGAAGTGATGAAAGATGGCATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAAc                        >  3:2808312/1‑149 (MQ=255)
                                                                                                                                aCGGGAAGTGATGAAAGATGGCATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACTg                      >  1:6420050/1‑149 (MQ=255)
                                                                                                                                aCGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCggcggcggc                                                                                                                          >  2:6260244/1‑49 (MQ=255)
                                                                                                                                  ggaaaGTGATGAAAGATGGCATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTACTGCTGTGGGGTGTGTCATGCGCCATTGAAACTGcc                    >  4:2070035/4‑149 (MQ=255)
                                                                                                                                     aaGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACTGCCgta                 <  2:1717652/149‑3 (MQ=255)
                                                                                                                                          atgaAAGATGGCATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGaa                                                                                                     <  4:531763/60‑1 (MQ=255)
                                                                                                                                          atgaAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACTGCCGGAGACGc            >  1:1044819/1‑149 (MQ=255)
                                                                                                                                           tgaAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACTGCCGGAGACGCa           >  3:6270137/1‑149 (MQ=255)
                                                                                                                                           tgaAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACTGCCGGAGACGCa           >  4:5960841/1‑149 (MQ=255)
                                                                                                                                            gaAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGCGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACTGCCGGAGACGCAc          >  4:1457696/1‑149 (MQ=255)
                                                                                                                                            gaAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGgt                                            <  3:3385698/115‑1 (MQ=255)
                                                                                                                                               aGATGGCATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACTGCCGGAGACGCAccag       >  1:517599/1‑149 (MQ=255)
                                                                                                                                                  tGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACTGCCGGAGACGCAccagcca    >  4:661059/1‑149 (MQ=255)
                                                                                                                                                    aCATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGcc                                >  4:4975389/1‑119 (MQ=255)
                                                                                                                                                    aCATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACTGCCGGAGACGCACcagccagc  <  1:2881253/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGGCGGCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACTGCCGGAGACGCACCAGCCAGC  >  NC_000913/4146441‑4146737

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: