Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,043,595 C→A 5.0% M144I (ATG→ATT 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,5950CA 5.0% 1568.7 / 34.1 393M144I (ATG→ATTgcvTaminomethyltransferase, tetrahydrofolate‑dependent, subunit (T protein) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base C (178/194);  new base A (15/6);  total (193/200)
Fisher's exact test for biased strand distribution p-value = 4.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 398 total aligned reads displayed.

CCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCAC  >  CP009273/3043446‑3043744
                                                                                                                                                     |                                                                                                                                                     
   gCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATg                                                                                                                                                    >  2:418059/1‑150 (MQ=255)
    cTTCACCGGTAGAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATgg                                                                                                                                                   <  2:4659125/150‑1 (MQ=255)
     ttCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGa                                                                                                                                                  <  1:71818/150‑1 (MQ=255)
     ttCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGa                                                                                                                                                  >  1:4345489/1‑150 (MQ=255)
          cGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGt                                                                                                                                             <  1:4770734/150‑1 (MQ=255)
           ggTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGtc                                                                                                                                            >  1:3756556/1‑150 (MQ=255)
           ggTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATAATGGAAAGGtc                                                                                                                                            <  2:77120/150‑1 (MQ=255)
            gTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGtca                                                                                                                                           <  1:1638046/150‑1 (MQ=255)
             tataACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGtcat                                                                                                                                          <  2:1378277/150‑1 (MQ=255)
             tataACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGtcat                                                                                                                                          <  1:335534/150‑1 (MQ=255)
              ataACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGtcatc                                                                                                                                         >  2:1079238/1‑150 (MQ=255)
              ataACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGtcatc                                                                                                                                         >  1:49316/1‑150 (MQ=255)
              ataACCAGTGGTGGCAATAAACAGATAGCCCGCCAGCACGCCAAAGAACGGAATCATCCCTTCCACCGCCTGACGCTGGGCGTCAATAAACAGTGTGGCAGCTTTTGCCAGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGtcatc                                                                                                                                         >  2:3651636/1‑150 (MQ=255)
               taACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGtcatca                                                                                                                                        <  1:4211451/150‑1 (MQ=255)
                aaCCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCAc                                                                                                                                       >  2:459941/1‑150 (MQ=255)
                aaCCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCAc                                                                                                                                       <  2:3795183/150‑1 (MQ=255)
                aaCCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCAc                                                                                                                                       >  2:2402128/1‑150 (MQ=255)
                  ccAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGa                                                                                                                                     >  1:3124919/1‑150 (MQ=255)
                    attggtggCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAAc                                                                                                                                   >  2:3564867/3‑150 (MQ=255)
                     gtggtgGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACg                                                                                                                                  >  2:1083430/1‑150 (MQ=255)
                     gtggtgGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACg                                                                                                                                  >  2:1148275/1‑150 (MQ=255)
                     gtggtgGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACg                                                                                                                                  <  2:3310723/150‑1 (MQ=255)
                        gtggCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTAGCACGGCACTAATGGAAAGGTCATCACGAACGGTa                                                                                                                               >  2:5109675/1‑150 (MQ=255)
                         tggCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTaa                                                                                                                              <  2:4300492/150‑1 (MQ=255)
                           gCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAAtt                                                                                                                            <  2:1401247/150‑1 (MQ=255)
                            cAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAAttt                                                                                                                           <  1:2801774/150‑1 (MQ=255)
                              aTAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCg                                                                                                                         >  1:3160208/1‑150 (MQ=255)
                               tAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCCCGAACGGTAATTTCGa                                                                                                                        <  2:2611836/150‑1 (MQ=255)
                                  aCAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGc                                                                                                                     <  2:151957/150‑1 (MQ=255)
                                   cAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTa                                                                                                                               >  1:878795/1‑139 (MQ=255)
                                   cAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGcc                                                                                                                    <  1:4054896/150‑1 (MQ=255)
                                   cAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGcc                                                                                                                    >  1:3794509/1‑150 (MQ=255)
                                   cAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATAATGGAAAGGTCATCACGAACGGTAATTTCGATGcc                                                                                                                    <  1:120195/150‑1 (MQ=255)
                                     gATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGa                                                                                                                  >  2:2029350/1‑150 (MQ=255)
                                          cccgccTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGt                                                                                                             >  2:169543/1‑150 (MQ=255)
                                          cccgccTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGt                                                                                                             <  1:3486217/150‑1 (MQ=255)
                                           ccgccTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGtt                                                                                                            <  1:4609376/150‑1 (MQ=255)
                                               cTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGc                                                                                                        <  1:231071/150‑1 (MQ=255)
                                                tGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCg                                                                                                       >  2:492588/1‑150 (MQ=255)
                                                 gCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCgt                                                                                                      <  1:3856470/150‑1 (MQ=255)
                                                   aCGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCgtgt                                                                                                    >  1:3316854/1‑150 (MQ=255)
                                                    cGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGtt                                                                                                   >  1:566603/1‑150 (MQ=255)
                                                      ccAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTgg                                                                                                 >  1:1927347/1‑150 (MQ=255)
                                                      ccAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTgg                                                                                                 >  2:2091142/1‑150 (MQ=255)
                                                       cAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTggg                                                                                                >  1:2723686/1‑150 (MQ=255)
                                                         aaGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTa                                                                                              >  2:4713208/1‑150 (MQ=255)
                                                               ggTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAg                                                                                        >  2:1249070/1‑150 (MQ=255)
                                                                gTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAgg                                                                                       <  2:3451942/150‑1 (MQ=255)
                                                                 tttCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAACGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTGGGGTAATCCAGga                                                                                      >  2:5007062/1‑150 (MQ=255)
                                                                 tttCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGga                                                                                      >  2:4849122/1‑150 (MQ=255)
                                                                  ttCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGgag                                                                                     <  1:3599884/150‑1 (MQ=255)
                                                                   tCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGgaga                                                                                    >  2:2969793/1‑150 (MQ=255)
                                                                     aTCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAgg                                                                                  >  1:3429881/1‑150 (MQ=255)
                                                                     aTCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAgg                                                                                  >  2:2803353/1‑150 (MQ=255)
                                                                        ccTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCt                                                                               <  1:3237171/150‑1 (MQ=255)
                                                                         cTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATAATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCtt                                                                              >  1:940212/1‑150 (MQ=255)
                                                                          ttCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCttt                                                                             <  2:559487/150‑1 (MQ=255)
                                                                               ccGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTcgc                                                                        >  2:2880872/1‑150 (MQ=255)
                                                                               ccGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTcgc                                                                        >  2:405587/1‑150 (MQ=255)
                                                                               ccGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTcgc                                                                        >  2:52035/1‑150 (MQ=255)
                                                                                  ccTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCggg                                                                     <  1:739041/150‑1 (MQ=255)
                                                                                   cTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGt                                                                    >  2:4141953/1‑150 (MQ=255)
                                                                                   cTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGt                                                                    >  1:2521416/1‑150 (MQ=255)
                                                                                   cTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGt                                                                    >  1:2917429/1‑150 (MQ=255)
                                                                                    tGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTg                                                                   >  2:2582571/1‑150 (MQ=255)
                                                                                     gACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTgg                                                                  >  1:677496/1‑150 (MQ=255)
                                                                                      aCGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGc                                                                 <  2:2850510/150‑1 (MQ=255)
                                                                                         cTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGa                                                              >  1:2682259/1‑150 (MQ=255)
                                                                                         cTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGa                                                              >  1:1242587/1‑150 (MQ=255)
                                                                                           gggCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGt                                                            <  1:169543/150‑1 (MQ=255)
                                                                                              cGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTGATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTaa                                                         <  1:1379039/150‑1 (MQ=255)
                                                                                              cGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTaa                                                         <  1:3653531/150‑1 (MQ=255)
                                                                                              cGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTa                                                          >  2:4941029/1‑149 (MQ=255)
                                                                                                tCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTaaca                                                       <  2:424579/150‑1 (MQ=255)
                                                                                                 cATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTaacaa                                                      >  2:1470559/1‑150 (MQ=255)
                                                                                                   ttAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACg                                                    <  2:1918369/150‑1 (MQ=255)
                                                                                                      aaCAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAgg                                                 >  1:1161813/1‑150 (MQ=255)
                                                                                                        cAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCg                                               >  2:686709/1‑150 (MQ=255)
                                                                                                         aGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCgg                                              >  2:2459265/1‑150 (MQ=255)
                                                                                                            gtgGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGgaag                                           <  1:2404061/150‑1 (MQ=255)
                                                                                                             tgGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATAATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGgaaga                                          <  2:3097998/150‑1 (MQ=255)
                                                                                                             tgGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATAATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGgaaga                                          <  2:1203721/150‑1 (MQ=255)
                                                                                                                 aGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATc                                      <  1:2117885/150‑1 (MQ=255)
                                                                                                                   cTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCtt                                    <  2:1073286/150‑1 (MQ=255)
                                                                                                                       tGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGt                                <  1:2204266/150‑1 (MQ=255)
                                                                                                                         ccTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTaa                              >  1:3052740/1‑150 (MQ=255)
                                                                                                                           tGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAag                            <  2:5055661/150‑1 (MQ=255)
                                                                                                                              gcATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAagtag                         <  1:2943097/150‑1 (MQ=255)
                                                                                                                                aTTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAgtagta                       >  1:2047889/1‑150 (MQ=255)
                                                                                                                                  tCGGCCCTTGCACGGCAATCATGGCAAGGTCATCACGAACGGGAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACCACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGa                     <  1:2659735/150‑1 (MQ=255)
                                                                                                                                      cccTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGAt                 >  2:3948401/1‑150 (MQ=255)
                                                                                                                                       ccTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATg                <  2:4690348/150‑1 (MQ=255)
                                                                                                                                        cTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGa               >  1:2444056/1‑150 (MQ=255)
                                                                                                                                        cTTGCACGGCAATAATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGa               <  1:151729/150‑1 (MQ=255)
                                                                                                                                         ttGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAg              <  1:418059/150‑1 (MQ=255)
                                                                                                                                            cACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTc           >  1:5030642/1‑150 (MQ=255)
                                                                                                                                                gCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGAAAAGTAGTAGACGATGAGGTCATCt       >  2:4528795/1‑150 (MQ=255)
                                                                                                                                                  aaTCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTAt     >  2:4398673/1‑150 (MQ=255)
                                                                                                                                                   aTCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAGAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATc    >  1:1425712/1‑150 (MQ=255)
                                                                                                                                                   aTCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATc    <  2:656138/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCGCTTCACCGGTATAACCAGTGGTGGCAATAAACAGATCGCCCGCCTGCACGCCAAAGAACGGTTTCATCCCTTCCACCGCCTGACGCTGGGCGTCATTAAACAGTGTGGCAGCTTTTGCCTGCGCATTCGGCCCTTGCACGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCAC  >  CP009273/3043446‑3043744

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: