Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,043,612 G→T 6.1% R139S (CGT→AGT)  gcvT ← aminomethyltransferase, tetrahydrofolate‑dependent, subunit (T protein) of glycine cleavage complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,043,6120GT 6.1% 762.9 / 15.8 203R139S (CGT→AGT) gcvTaminomethyltransferase, tetrahydrofolate‑dependent, subunit (T protein) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base G (84/106);  new base T (7/6);  total (91/112)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

AGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGC  >  CP009273/3043466‑3043756
                                                                                                                                                  |                                                                                                                                                
  tggtggCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAAc                                                                                                                                               <  2:1973774/148‑1 (MQ=255)
      ggCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTa                                                                                                                                           <  1:2503182/148‑1 (MQ=255)
      ggCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTa                                                                                                                                           >  2:2181049/1‑148 (MQ=255)
      ggCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTa                                                                                                                                           >  2:3440742/1‑148 (MQ=255)
      ggCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCAGTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAACGGTCATCACGAACGGTa                                                                                                                                           >  2:3441229/1‑148 (MQ=255)
       gCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTaa                                                                                                                                          <  1:78919/148‑1 (MQ=255)
       gCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTaa                                                                                                                                          <  1:3546023/148‑1 (MQ=255)
         aaTAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAAtt                                                                                                                                        <  2:624881/148‑1 (MQ=255)
          aTAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAAttt                                                                                                                                       <  1:1297336/148‑1 (MQ=255)
            aaaCAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCg                                                                                                                                     >  2:2039738/1‑148 (MQ=255)
            aaaCAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCg                                                                                                                                     >  2:3473211/1‑148 (MQ=255)
            aaaCAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTCGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCg                                                                                                                                     >  2:2529726/1‑148 (MQ=255)
             aaCAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGa                                                                                                                                    <  2:1805214/148‑1 (MQ=255)
             aaCAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGa                                                                                                                                    >  1:2965466/1‑148 (MQ=255)
              aCAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGAt                                                                                                                                   >  2:3392118/1‑148 (MQ=255)
               cAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATg                                                                                                                                  <  1:1468922/148‑1 (MQ=255)
               cAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATg                                                                                                                                  >  1:2344935/1‑148 (MQ=255)
                    cgcccgccTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGaa                                                                                                                             >  2:2158041/1‑148 (MQ=255)
                     gcccgccTGCTCGCCAAAGATCGGTTTCATCCCTTCCACCGCCTGACGCTGGGAGTCATTAAACAGTGTGGCAGCTTTTTCCTGCGCATTCTGCCCTTGCACGGCAATCATGGAAAGGTCATCACGATCGGTAATTTCGATGCCGaaa                                                                                                                            <  2:2236226/148‑1 (MQ=255)
                        cgccTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGt                                                                                                                         >  1:2197148/1‑148 (MQ=255)
                          ccTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTc                                                                                                                       <  1:2334707/148‑1 (MQ=255)
                            tGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAg                                                                                                                     <  2:279086/148‑1 (MQ=255)
                            tGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAg                                                                                                                     <  1:1761884/148‑1 (MQ=255)
                                 gCCTAAGAAAGTTTTCATCCCTTCCACCGCCTGACGCTTGTGGTCATTAAACAGTGTGTCATCTTTTTCCTGCGCATTAGTCAGTTGCACGGCAATCATGGTAAGGTCATAACTAACGGTATTTTCGATGCCTAAAGGTTCAGCgttt                                                                                                                <  2:2270677/148‑3 (MQ=255)
                                 gCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCgtgt                                                                                                                <  1:3038916/148‑1 (MQ=255)
                                   cAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGTCGTAATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTg                                                                                                              <  2:2959846/148‑1 (MQ=255)
                                   cAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTg                                                                                                              >  2:1447289/1‑148 (MQ=255)
                                   cAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTg                                                                                                              <  2:951784/148‑1 (MQ=255)
                                     aaGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTggg                                                                                                            <  1:2354499/148‑1 (MQ=255)
                                         aCGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACTAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAAt                                                                                                        >  1:2157002/1‑148 (MQ=255)
                                          cGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATc                                                                                                       <  1:881064/148‑1 (MQ=255)
                                           ggTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATcc                                                                                                      >  2:355443/1‑148 (MQ=255)
                                             tttCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAg                                                                                                    >  2:3155507/1‑148 (MQ=255)
                                                cATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGgag                                                                                                 <  1:2223394/148‑1 (MQ=255)
                                                 aTCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGgaga                                                                                                >  1:2545782/1‑148 (MQ=255)
                                                 aTCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGgaga                                                                                                <  1:1411552/148‑1 (MQ=255)
                                                 aTCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGgaga                                                                                                <  1:2421620/148‑1 (MQ=255)
                                                  tCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGagag                                                                                               <  1:2868735/148‑1 (MQ=255)
                                                    ccTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGTAATCCAGGAGAGGTc                                                                                            <  1:659239/148‑1 (MQ=255)
                                                      ttCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCt                                                                                           <  1:3265256/148‑1 (MQ=255)
                                                         cACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCtttt                                                                                        <  2:2190578/148‑1 (MQ=255)
                                                           ccGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTc                                                                                      >  1:237344/1‑148 (MQ=255)
                                                                tGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCggg                                                                                 <  1:2158041/148‑1 (MQ=255)
                                                                 gACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGt                                                                                <  1:3210062/148‑1 (MQ=255)
                                                                 gACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGt                                                                                <  1:1367833/148‑1 (MQ=255)
                                                                       gggCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGa                                                                          >  1:914485/1‑148 (MQ=255)
                                                                         gCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGt                                                                        <  1:3486347/148‑1 (MQ=255)
                                                                          cGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGtt                                                                       >  1:769051/1‑148 (MQ=255)
                                                                          cGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGtt                                                                       >  2:1320624/1‑148 (MQ=255)
                                                                           gTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTa                                                                      >  2:2341236/1‑148 (MQ=255)
                                                                              aTTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTaaca                                                                   <  2:961454/148‑1 (MQ=255)
                                                                               ttAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTaacaa                                                                  <  1:513601/148‑1 (MQ=255)
                                                                               ttAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTaacaa                                                                  <  1:306081/148‑1 (MQ=255)
                                                                                tAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTaacaac                                                                 >  2:1276599/1‑148 (MQ=255)
                                                                                  aaCAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGa                                                               >  1:531301/1‑148 (MQ=255)
                                                                                    cAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAgg                                                             <  2:3250092/148‑1 (MQ=255)
                                                                                       tgtgGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCgg                                                          <  1:346680/148‑1 (MQ=255)
                                                                                           gCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGgaaga                                                      >  1:483134/1‑148 (MQ=255)
                                                                                           gCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGgaaga                                                      >  2:2462793/1‑148 (MQ=255)
                                                                                           gCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGgaaga                                                      <  1:2018506/148‑1 (MQ=255)
                                                                                            cAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACTAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGgaagaa                                                     >  2:435717/1‑148 (MQ=255)
                                                                                              gCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGTGGGTGGCGGAGTTAACAACGAGGCGGAAGAAAt                                                   <  2:526373/148‑1 (MQ=255)
                                                                                              gCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACg                                                                >  1:1088353/1‑135 (MQ=255)
                                                                                              gCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACg                                                                <  2:1088353/135‑1 (MQ=255)
                                                                                              gCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAAt                                                   >  2:1220671/1‑148 (MQ=255)
                                                                                                  ttGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTc                                               >  1:1875423/1‑148 (MQ=255)
                                                                                                  ttGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTc                                               <  2:141597/148‑1 (MQ=255)
                                                                                                    gCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAg                                             <  1:170859/148‑1 (MQ=255)
                                                                                                      cTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTa                                           >  1:950692/1‑148 (MQ=255)
                                                                                                       tGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTaa                                          >  2:3081874/1‑148 (MQ=255)
                                                                                                          gcATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTCGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAagt                                       <  1:3392118/148‑1 (MQ=255)
                                                                                                              tCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACTAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAgtagta                                   >  1:1027259/1‑148 (MQ=255)
                                                                                                               cGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGtagtag                                  <  1:2772653/148‑1 (MQ=255)
                                                                                                                 gCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGAc                                >  1:2752154/1‑148 (MQ=255)
                                                                                                                  cccTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACg                               <  2:994569/148‑1 (MQ=255)
                                                                                                                  cccTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACg                               <  2:1606116/148‑1 (MQ=255)
                                                                                                                  cccTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACg                               <  2:1515022/148‑1 (MQ=255)
                                                                                                                   ccTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTCATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGACGAAATCTTCAGTAAAGTAGTAGACGa                              >  2:1844495/1‑148 (MQ=255)
                                                                                                                   ccTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAGGAAATCTTCAGTAAAGTAGTAGACGa                              >  1:1757924/1‑148 (MQ=255)
                                                                                                                   ccTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGa                              <  2:1828887/148‑1 (MQ=255)
                                                                                                                   ccTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGa                              >  2:2121240/1‑148 (MQ=255)
                                                                                                                   ccTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGa                              >  2:329544/1‑148 (MQ=255)
                                                                                                                    cTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGAt                             >  2:251599/1‑148 (MQ=255)
                                                                                                                     ttGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATg                            <  1:1507416/148‑1 (MQ=255)
                                                                                                                      tGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGa                           >  1:1472382/1‑148 (MQ=255)
                                                                                                                      tGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGa                           >  1:2642555/1‑148 (MQ=255)
                                                                                                                       gCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAg                          <  1:972687/148‑1 (MQ=255)
                                                                                                                            gCAATCATGGAAAGGTCATCACTAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAt                     <  2:2157002/148‑1 (MQ=255)
                                                                                                                            gCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAt                     <  1:3440742/148‑1 (MQ=255)
                                                                                                                            gCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAt                     <  1:2181049/148‑1 (MQ=255)
                                                                                                                            gCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAt                     <  1:3441229/148‑1 (MQ=255)
                                                                                                                             cAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATc                    >  2:70734/1‑148 (MQ=255)
                                                                                                                               aTCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCTGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTa                  <  2:2724720/148‑1 (MQ=255)
                                                                                                                                tCATGGAAAGGTCATCACTAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTAt                 >  1:1219745/1‑148 (MQ=255)
                                                                                                                                    ggAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATcaca             <  1:1425427/148‑1 (MQ=255)
                                                                                                                                    ggAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATcaca             <  2:914529/148‑1 (MQ=255)
                                                                                                                                     gAAAGGTCATCACGAACTGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCacac            >  1:1850217/1‑148 (MQ=255)
                                                                                                                                              tcaCGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAgg   <  2:1498851/148‑1 (MQ=255)
                                                                                                                                              tcaCGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAgg   >  1:2568615/1‑148 (MQ=255)
                                                                                                                                              tcaCGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAgg   <  2:1498715/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
AGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGC  >  CP009273/3043466‑3043756

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