Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,150,054 G→T 100% S222Y (TCT→TAT)  sthA ← pyridine nucleotide transhydrogenase, soluble

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,150,0540GT98.7% 995.8 / 2.1 320S222Y (TCT→TAT) sthApyridine nucleotide transhydrogenase, soluble
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base T (178/138);  total (178/142)
Fisher's exact test for biased strand distribution p-value = 3.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 326 total aligned reads displayed.

GCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGG  >  CP009273/4149906‑4150200
                                                                                                                                                    |                                                                                                                                                  
gCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGgagagaga                                                                                                                                                   >  1:1461318/1‑150 (MQ=255)
  gACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAAt                                                                                                                                                 <  1:4754787/150‑1 (MQ=255)
  gACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCAACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAAt                                                                                                                                                 <  2:1579576/150‑1 (MQ=255)
   aCCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATc                                                                                                                                                <  2:2400003/150‑1 (MQ=255)
    ccGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCt                                                                                                                                               <  2:89370/150‑1 (MQ=255)
    ccGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCt                                                                                                                                               <  1:5105108/150‑1 (MQ=255)
     cGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTg                                                                                                                                              <  1:4504678/150‑1 (MQ=255)
      gTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGa                                                                                                                                             >  1:555305/1‑150 (MQ=255)
            aTAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACAtctc                                                                                                                                       <  1:4296162/150‑1 (MQ=255)
             tAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATctct                                                                                                                                      <  1:5262924/150‑1 (MQ=255)
              agagCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCtt                                                                                                                                     <  1:3320240/150‑1 (MQ=255)
              agagCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCtt                                                                                                                                     >  2:242341/1‑150 (MQ=255)
              agagCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCtt                                                                                                                                     >  2:1597765/1‑150 (MQ=255)
                   aggcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCg                                                                                                                                <  1:1495992/150‑1 (MQ=255)
                        gTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGaaa                                                                                                                           >  2:270008/1‑150 (MQ=255)
                        gTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGaaa                                                                                                                           >  2:2150408/1‑150 (MQ=255)
                        gTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGaaa                                                                                                                           >  2:1717933/1‑150 (MQ=255)
                          cAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATg                                                                                                                         >  1:2218631/1‑150 (MQ=255)
                           aGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGc                                                                                                                        >  1:2524521/1‑150 (MQ=255)
                           aGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGc                                                                                                                        >  1:3154438/1‑150 (MQ=255)
                             cTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCca                                                                                                                      >  1:5184427/1‑150 (MQ=255)
                             cTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCca                                                                                                                      >  2:2025963/1‑150 (MQ=255)
                             cTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCca                                                                                                                      <  2:563392/150‑1 (MQ=255)
                                  aGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAgg                                                                                                                 >  1:358372/1‑150 (MQ=255)
                                  aGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAgg                                                                                                                 >  1:3823742/1‑150 (MQ=255)
                                  aGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGAGAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCTgg                                                                                                                 >  2:2749027/1‑150 (MQ=255)
                                    ttttttACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCg                                                                                                               <  2:2652988/150‑1 (MQ=255)
                                       tttACCCGACTTCAGATGCATGATCACACCTTCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGTGATATGCCAGCAGGCGATc                                                                                                            <  2:106015/150‑1 (MQ=255)
                                         tACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATcgc                                                                                                          >  1:409220/1‑150 (MQ=255)
                                          aCCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATcgcg                                                                                                         <  2:395836/150‑1 (MQ=255)
                                                  tCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGa                                                                                                 <  1:235016/150‑1 (MQ=255)
                                                   cAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGAt                                                                                                <  2:506108/150‑1 (MQ=255)
                                                    aGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATc                                                                                               <  1:1052723/150‑1 (MQ=255)
                                                    aGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATc                                                                                               >  2:2933802/1‑150 (MQ=255)
                                                    aGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATc                                                                                               <  1:3034337/150‑1 (MQ=255)
                                                    aGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATc                                                                                               >  2:4627131/1‑150 (MQ=255)
                                                      aTGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAg                                                                                             >  2:249900/1‑150 (MQ=255)
                                                       tGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATCGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGa                                                                                            >  2:3050361/1‑150 (MQ=255)
                                                       tGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGa                                                                                            <  1:4354199/150‑1 (MQ=255)
                                                       tGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGa                                                                                            >  2:4774150/1‑150 (MQ=255)
                                                         cATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATc                                                                                          >  1:957582/1‑150 (MQ=255)
                                                           tGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCa                                                                                        >  1:4154050/1‑150 (MQ=255)
                                                            gATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCAc                                                                                       <  2:1928805/150‑1 (MQ=255)
                                                             aTCACACCATCGTCACGGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACt                                                                                      <  1:4607241/150‑1 (MQ=255)
                                                             aTCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACt                                                                                      >  2:4415424/1‑150 (MQ=255)
                                                             aTCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACt                                                                                      <  2:1295473/150‑1 (MQ=255)
                                                             aTCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACt                                                                                      <  1:4434364/150‑1 (MQ=255)
                                                              tCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCACAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACtt                                                                                     >  2:2524642/1‑150 (MQ=255)
                                                               cacaCCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACttt                                                                                    <  1:642060/150‑1 (MQ=255)
                                                               cacaCCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACttt                                                                                    <  2:1258444/150‑1 (MQ=255)
                                                               cacaCCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACttt                                                                                    <  2:1025459/150‑1 (MQ=255)
                                                                acacCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACtttt                                                                                   <  2:2300726/150‑1 (MQ=255)
                                                                 cacCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTa                                                                                  <  2:3426721/150‑1 (MQ=255)
                                                                 cacCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTa                                                                                  >  1:1179629/1‑150 (MQ=255)
                                                                  acCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTAc                                                                                 >  1:1830162/1‑150 (MQ=255)
                                                                     aTCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATc                                                                              >  1:338051/1‑150 (MQ=255)
                                                                       cGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCa                                                                            >  2:2762038/1‑150 (MQ=255)
                                                                         tCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATCGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATa                                                                          >  2:314338/1‑150 (MQ=255)
                                                                         tCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCCGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATa                                                                          >  2:3630728/1‑150 (MQ=255)
                                                                         tCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATa                                                                          <  2:498809/150‑1 (MQ=255)
                                                                         tCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATa                                                                          <  2:1510571/150‑1 (MQ=255)
                                                                            caGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACcg                                                                       >  2:4282977/1‑150 (MQ=255)
                                                                               ccTTCGATCTTCTCGAACTCTTCGTTGTGACGAATCACTACGCCACTGATCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCTATCGCGGGTGTTGATCAGATCCACTTTTACATCCAAACCGCgg                                                                    >  2:3113925/1‑150 (MQ=255)
                                                                                cTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGa                                                                   >  1:1914718/1‑150 (MQ=255)
                                                                                 ttCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGaa                                                                  >  1:2169643/1‑150 (MQ=255)
                                                                                    gATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGAt                                                               >  2:4165274/1‑150 (MQ=255)
                                                                                    gATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGAt                                                               >  1:189216/1‑150 (MQ=255)
                                                                                    gATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGAt                                                               >  2:4576949/1‑150 (MQ=255)
                                                                                     aTCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATc                                                              <  1:3754202/150‑1 (MQ=255)
                                                                                       cttctCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGa                                                            >  2:133214/1‑150 (MQ=255)
                                                                                       cttctCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGa                                                            >  1:4438713/1‑150 (MQ=255)
                                                                                        ttctCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGAc                                                           <  1:3237535/150‑1 (MQ=255)
                                                                                          ctcGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGc                                                         >  2:3322337/1‑150 (MQ=255)
                                                                                          ctcGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGAGAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGc                                                         >  2:774698/1‑150 (MQ=255)
                                                                                            cGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCat                                                       >  2:4867178/1‑150 (MQ=255)
                                                                                              tACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCatat                                                     <  1:1247199/150‑1 (MQ=255)
                                                                                                 tctTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTca                                                  <  1:4052074/150‑1 (MQ=255)
                                                                                                    tCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAg                                               <  1:1597765/150‑1 (MQ=255)
                                                                                                     cGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGc                                              >  1:1341893/1‑150 (MQ=255)
                                                                                                         gtgACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGAt                                          <  1:3931798/150‑1 (MQ=255)
                                                                                                          tgACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATc                                         <  2:1742421/150‑1 (MQ=255)
                                                                                                           gACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCa                                        <  1:1731792/150‑1 (MQ=255)
                                                                                                              gAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTc                                     <  2:3823742/150‑1 (MQ=255)
                                                                                                               aaTCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTcc                                    >  2:541501/1‑150 (MQ=255)
                                                                                                                aTCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCa                                   >  1:752106/1‑150 (MQ=255)
                                                                                                                aTCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCa                                   <  1:3510832/150‑1 (MQ=255)
                                                                                                                     tACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCAcc                              >  1:1693093/1‑150 (MQ=255)
                                                                                                                       cGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCAt                            >  1:4990005/1‑150 (MQ=255)
                                                                                                                       cGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGAAGATCAGATCCACTTTAACATCCAGACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCaa                            >  2:4452160/1‑149 (MQ=255)
                                                                                                                        gCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATa                           <  1:4921113/150‑1 (MQ=255)
                                                                                                                           aCTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGAt                        <  2:2840722/150‑1 (MQ=255)
                                                                                                                            cTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATa                       <  2:1863236/150‑1 (MQ=255)
                                                                                                                               ttCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGt                    >  1:2784470/1‑150 (MQ=255)
                                                                                                                               ttCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGt                    <  2:668715/150‑1 (MQ=255)
                                                                                                                               ttCAAGAAGTGATAGGAGAGATAATCTGACATCAAATGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATACACAGCCGATCACTCCAGCACCATAGATAAGt                    >  2:4925936/1‑150 (MQ=255)
                                                                                                                                   aGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACAt                <  1:2901361/150‑1 (MQ=255)
                                                                                                                                     aaGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATgg              >  1:2545705/1‑150 (MQ=255)
                                                                                                                                        tGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGgcgc           <  2:5222525/150‑1 (MQ=255)
                                                                                                                                          aTAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTGTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCgg         >  2:4350030/1‑150 (MQ=255)
                                                                                                                                                gagataATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCgtg   >  1:2032145/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGG  >  CP009273/4149906‑4150200

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: