Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,257,764 C→A 13.1% D205Y (GAC→TAC)  tdcC ← L‑threonine/L‑serine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,257,7640CA13.1% 942.3 / 82.0 287D205Y (GAC→TAC) tdcCL‑threonine/L‑serine transporter
Reads supporting (aligned to +/- strand):  ref base C (113/136);  new base A (27/11);  total (140/147)
Fisher's exact test for biased strand distribution p-value = 4.80e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 287 total aligned reads displayed.

GGGAACATTTACGTTCGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCAGGTAGCTCATCACTTTAACCATCAGATCCTTACCA  >  CP009273/3257617‑3257912
                                                                                                                                                   |                                                                                                                                                    
  gAACATTTACGTTCGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGa                                                                                                                                                  <  2:889277/150‑1 (MQ=255)
  gAACATTTACGTTCGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTACACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGa                                                                                                                                                  <  2:1330690/150‑1 (MQ=255)
   aaCATTTACGTTCGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGAc                                                                                                                                                 <  2:3048128/150‑1 (MQ=255)
   aaCATTTACGTTCGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGAc                                                                                                                                                 >  1:402322/1‑150 (MQ=255)
   aaCATTTACGTTCGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGAc                                                                                                                                                 <  1:3354448/150‑1 (MQ=255)
    aCATTTACGTTCGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGAcc                                                                                                                                                <  1:2360783/150‑1 (MQ=255)
        ttACGTTCGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGtt                                                                                                                                            <  1:1536853/150‑1 (MQ=255)
        ttACGTTCGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGtt                                                                                                                                            <  2:3999192/150‑1 (MQ=255)
        ttACGTTCGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGtt                                                                                                                                            <  1:937188/150‑1 (MQ=255)
        ttACGTTCGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGtt                                                                                                                                            <  1:2195615/150‑1 (MQ=255)
        ttACGTTCGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGtt                                                                                                                                            <  1:2426582/150‑1 (MQ=255)
        ttACGTTCGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGtt                                                                                                                                            <  1:2896700/150‑1 (MQ=255)
            gTTCGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACg                                                                                                                                        >  2:1016333/1‑150 (MQ=255)
             ttCGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGa                                                                                                                                       >  1:2534932/1‑150 (MQ=255)
               cGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACa                                                                                                                                     <  1:3846371/150‑1 (MQ=255)
               cGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACa                                                                                                                                     <  2:4425014/150‑1 (MQ=255)
               cGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACa                                                                                                                                     <  2:2186267/150‑1 (MQ=255)
                ggTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAg                                                                                                                                    >  2:2491682/1‑150 (MQ=255)
                   gAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGa                                                                                                                                 >  1:4511193/1‑150 (MQ=255)
                      gTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAAcc                                                                                                                              <  2:1943125/150‑1 (MQ=255)
                         gcgACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAg                                                                                                                           <  2:4141664/150‑1 (MQ=255)
                          cgACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAgg                                                                                                                          >  2:514658/1‑150 (MQ=255)
                          cgACCGAAGTCTTTCTCATACTCTTCACGATTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAgg                                                                                                                          >  2:3805631/1‑150 (MQ=255)
                             ccGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCa                                                                                                                       <  2:2980114/150‑1 (MQ=255)
                             ccGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTAATGACCGGTTAACGACAGCGAACCGAGGTCa                                                                                                                       >  1:4506631/1‑150 (MQ=255)
                               gAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTAATGACCGGTTAACGACAGCGAACCGAGGTCAAc                                                                                                                     >  2:2256232/1‑150 (MQ=255)
                                aaGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAAcc                                                                                                                    <  2:2051670/150‑1 (MQ=255)
                                aaGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAAcc                                                                                                                    >  1:3285537/1‑150 (MQ=255)
                                 aGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCt                                                                                                                   <  2:3452852/150‑1 (MQ=255)
                                 aGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCt                                                                                                                   <  2:1182271/150‑1 (MQ=255)
                                  gTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTg                                                                                                                  >  2:3352816/1‑150 (MQ=255)
                                  gTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTg                                                                                                                  <  2:1090696/150‑1 (MQ=255)
                                  gTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTg                                                                                                                  <  1:1893195/150‑1 (MQ=255)
                                     tttCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTc                                                                                                               <  2:468444/150‑1 (MQ=255)
                                     tttCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTc                                                                                                               <  1:1985353/150‑1 (MQ=255)
                                     tttCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTc                                                                                                               <  2:4085245/150‑1 (MQ=255)
                                      ttCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCg                                                                                                              >  2:1364368/1‑150 (MQ=255)
                                      ttCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCg                                                                                                              >  2:1239276/1‑150 (MQ=255)
                                      ttCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCg                                                                                                              >  1:4103502/1‑150 (MQ=255)
                                       tctcATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGa                                                                                                             >  1:1088664/1‑150 (MQ=255)
                                          cATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATaa                                                                                                          >  2:4428393/1‑150 (MQ=255)
                                          cATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTAATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATaa                                                                                                          >  1:502495/1‑150 (MQ=255)
                                           aTACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAAc                                                                                                         >  1:1294520/1‑150 (MQ=255)
                                            tACTCTTCCCGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACt                                                                                                        <  1:1670682/150‑1 (MQ=255)
                                            tACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACt                                                                                                        >  1:2985766/1‑150 (MQ=255)
                                            tACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACt                                                                                                        <  1:1670507/150‑1 (MQ=255)
                                              ctctTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGc                                                                                                      <  2:1315349/150‑1 (MQ=255)
                                                   cACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAACGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGt                                                                                                 >  1:3695611/1‑150 (MQ=255)
                                                      gCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTcc                                                                                              >  1:281210/1‑150 (MQ=255)
                                                        ttAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTAATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAg                                                                                            >  1:4423272/1‑150 (MQ=255)
                                                         tAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGt                                                                                           >  1:3632/1‑150 (MQ=255)
                                                         tAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGt                                                                                           >  1:1042088/1‑150 (MQ=255)
                                                          aGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTa                                                                                          <  2:3598378/150‑1 (MQ=255)
                                                           gAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTaa                                                                                         >  2:141113/1‑150 (MQ=255)
                                                                    aaggaagAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGc                                                                                >  2:355671/1‑150 (MQ=255)
                                                                    aaggaagAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGc                                                                                >  1:1362795/1‑150 (MQ=255)
                                                                      ggaagAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGa                                                                              <  1:3315934/150‑1 (MQ=255)
                                                                       gaagAGACGATTGGCGAGAAGTTAAAGGAGAACACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGAGCCGCGAc                                                                             >  2:2917594/1‑150 (MQ=255)
                                                                        aagAGACGATTGGCTAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACa                                                                            <  2:3544930/150‑1 (MQ=255)
                                                                        aagAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCTGTTAACGACAGCGAACCGAGGTCAACCTGGACGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACa                                                                            >  1:2922841/1‑150 (MQ=255)
                                                                         agagACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTAATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACaa                                                                           <  1:3216281/150‑1 (MQ=255)
                                                                          gagaCGATTGGCGAGAAGTTAAAGGATAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACaaa                                                                          >  1:3955230/1‑150 (MQ=255)
                                                                            gaCGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATACCCAGCCAGACAGTGATCAGGATACCGTAATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGa                                                                        <  2:2524845/150‑1 (MQ=255)
                                                                              cGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTAATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGaaa                                                                      >  1:574009/1‑150 (MQ=255)
                                                                                aTTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATc                                                                    >  2:3312544/1‑150 (MQ=255)
                                                                                        gaAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAg                                                            <  2:2588396/150‑1 (MQ=255)
                                                                                        gaAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAg                                                            <  1:635684/150‑1 (MQ=255)
                                                                                        gaAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAg                                                            <  2:3897548/150‑1 (MQ=255)
                                                                                           gTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCt                                                         >  2:4273268/1‑150 (MQ=255)
                                                                                           gTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCt                                                         <  2:2988585/150‑1 (MQ=255)
                                                                                           gTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTAATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCt                                                         <  1:2072735/150‑1 (MQ=255)
                                                                                             tAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTgg                                                       <  1:3773715/150‑1 (MQ=255)
                                                                                              aaaGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGc                                                      >  1:330452/1‑150 (MQ=255)
                                                                                                    gaAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATaaa                                                <  1:3160555/150‑1 (MQ=255)
                                                                                                     aaaaCCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAAc                                               <  2:2399524/150‑1 (MQ=255)
                                                                                                     aaaaCCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTAATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAAc                                               >  1:3552309/1‑150 (MQ=255)
                                                                                                       aaCCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACgg                                             >  1:1726879/1‑150 (MQ=255)
                                                                                                        aCCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGc                                            <  1:4144207/150‑1 (MQ=255)
                                                                                                        aCCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGc                                            >  2:27865/1‑150 (MQ=255)
                                                                                                          cATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCa                                          <  2:1691239/150‑1 (MQ=255)
                                                                                                           atgatgGAAATCCCCAGCCAGACAGTGATCAGGATACCGTAATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCAt                                         >  2:2886927/1‑150 (MQ=255)
                                                                                                            tgatgGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATa                                        <  1:444319/150‑1 (MQ=255)
                                                                                                            tgatgGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATa                                        >  2:3950988/1‑150 (MQ=255)
                                                                                                               tgGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCTACCTGTTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCa                                     <  2:3013453/150‑1 (MQ=255)
                                                                                                               tgGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCa                                     <  2:4480176/150‑1 (MQ=255)
                                                                                                                  aaaTCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCAGGt                                  <  2:262114/150‑1 (MQ=255)
                                                                                                                    aTCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCAGGTAg                                <  2:799155/150‑1 (MQ=255)
                                                                                                                             cagACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCAGGTAGCTCATCACt                       <  2:4470522/150‑1 (MQ=255)
                                                                                                                             cagACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCAGGTAGCTCATCACt                       >  1:3486207/1‑150 (MQ=255)
                                                                                                                               gACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCAGGTAGCTCATCACttt                     <  1:3033759/150‑1 (MQ=255)
                                                                                                                                   gTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGGCAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCAGGTAGCTCATCACTTTAAcc                 <  2:2275311/150‑1 (MQ=255)
                                                                                                                                   gTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGTAATCAGCACCAGGCTGGCGATAAACGGCCATACCAGGTAGCTCATCACTTTAAcc                 <  1:1364368/150‑1 (MQ=255)
                                                                                                                                   gTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCAGGTAGCTCATCACTTTAAcc                 <  2:4511193/150‑1 (MQ=255)
                                                                                                                                     gATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCAGGTAGCTCATCACTTTAACCAt               >  2:1985979/1‑150 (MQ=255)
                                                                                                                                     gATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCAGGTAGCTCATCACTTTAACCAt               >  2:1003884/1‑150 (MQ=255)
                                                                                                                                        cAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCAGGTAGCTCATCACTTTAACCATCAg            >  2:4208570/1‑150 (MQ=255)
                                                                                                                                           gATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGAGAACTGCAGAGTTCCAGTAAGGCATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCAGGTAGCTCATCTCTTTAACCAGCAGATc         >  2:4379058/1‑150 (MQ=255)
                                                                                                                                            aTACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCAGGTAGCTCATCACTTTAACCATCAGATcc        >  2:4049898/1‑150 (MQ=255)
                                                                                                                                             tACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCAGGTAGCTCATCACTTTAACCATCAGATCCt       >  2:1839541/1‑150 (MQ=255)
                                                                                                                                                 gTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCAGGTAGCTCATCACTTTAACCATCAGATCCTTAcc   >  1:2751614/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGGAACATTTACGTTCGGTGAAGTCGCGACCGAAGTCTTTCTCATACTCTTCACGCTTAGAAACCACGAAGGAAGAGACGATTGGCGAGAAGTTAAAGGAGAAAACCATGATGGAAATCCCCAGCCAGACAGTGATCAGGATACCGTCATGACCGGTTAACGACAGCGAACCGAGGTCAACCTGGTCGATAACTGCAGAGTTCCAGTAAGGGATCAGCGACAAAGAAATCAGCACCAGGCTGGCGATAAACGGCCATACCAGGTAGCTCATCACTTTAACCATCAGATCCTTACCA  >  CP009273/3257617‑3257912

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: