Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,405,946 T→G 94.3% Q43P (CAG→CCG)  yfbU ← UPF0304 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,405,9460TG94.3% 862.7 / 22.5 227Q43P (CAG→CCG) yfbUUPF0304 family protein
Reads supporting (aligned to +/- strand):  ref base T (9/5);  new base G (114/99);  total (123/104)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 227 total aligned reads displayed.

AAGGTGACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACT  >  CP009273/2405799‑2406092
                                                                                                                                                   |                                                                                                                                                  
aaGGTGACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGt                                                                                                                                                  >  1:3080764/1‑150 (MQ=255)
aaGGTGACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGt                                                                                                                                                  >  1:2504667/1‑150 (MQ=255)
 aGGTGACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTa                                                                                                                                                 <  1:2331539/150‑1 (MQ=255)
  ggTGACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTaa                                                                                                                                                <  1:3144479/150‑1 (MQ=255)
  ggTGACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTaa                                                                                                                                                <  1:3144942/150‑1 (MQ=255)
    tGACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATc                                                                                                                                              <  2:1230440/150‑1 (MQ=255)
         cgacgTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTaa                                                                                                                                         <  2:760768/150‑1 (MQ=255)
              ttCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACg                                                                                                                                    >  1:2173392/1‑150 (MQ=255)
               tCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGc                                                                                                                                   >  2:1006814/1‑150 (MQ=255)
                cATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCt                                                                                                                                  <  2:1552011/150‑1 (MQ=255)
                      tGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCaataa                                                                                                                            >  2:1267597/1‑150 (MQ=255)
                      tGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCaataa                                                                                                                            <  1:4308731/150‑1 (MQ=255)
                       ggATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCaataat                                                                                                                           <  2:875366/150‑1 (MQ=255)
                         aTTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTg                                                                                                                         <  2:4545921/150‑1 (MQ=255)
                          ttGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGt                                                                                                                        >  2:2222649/1‑150 (MQ=255)
                          ttGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGt                                                                                                                        >  1:948816/1‑150 (MQ=255)
                            gCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGttt                                                                                                                      <  2:3699741/150‑1 (MQ=255)
                             cTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTg                                                                                                                     <  1:3953891/150‑1 (MQ=255)
                              tGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGc                                                                                                                    >  1:3894436/1‑150 (MQ=255)
                              tGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGc                                                                                                                    >  2:3323533/1‑150 (MQ=255)
                              tGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGc                                                                                                                    >  1:3811390/1‑150 (MQ=255)
                              tGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATATCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGc                                                                                                                    >  1:4556623/1‑150 (MQ=255)
                                gTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAg                                                                                                                  <  2:1214243/150‑1 (MQ=255)
                                 tCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAgg                                                                                                                 <  1:4715907/150‑1 (MQ=255)
                                 tCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAgg                                                                                                                 <  2:196844/150‑1 (MQ=255)
                                   cTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCg                                                                                                               >  1:4155834/1‑150 (MQ=255)
                                   cTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCg                                                                                                               >  1:3527834/1‑150 (MQ=255)
                                    tGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGa                                                                                                              <  1:715262/150‑1 (MQ=255)
                                    tGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGa                                                                                                              >  1:2555578/1‑150 (MQ=255)
                                        aaTTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGt                                                                                                          >  1:4368970/1‑150 (MQ=255)
                                         aTTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTa                                                                                                         >  1:221992/1‑150 (MQ=255)
                                          ttAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTaa                                                                                                        <  1:1100900/150‑1 (MQ=255)
                                           tAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAAc                                                                                                       >  1:4527370/1‑150 (MQ=255)
                                               ccAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTc                                                                                                   <  1:1420530/150‑1 (MQ=255)
                                               ccAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTc                                                                                                   >  2:3540219/1‑150 (MQ=255)
                                               ccAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTc                                                                                                   >  1:281054/1‑150 (MQ=255)
                                                cAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCa                                                                                                  <  2:3801501/150‑1 (MQ=255)
                                                 aGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAg                                                                                                 >  1:3678327/1‑150 (MQ=255)
                                                 aGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAg                                                                                                 >  2:1584719/1‑150 (MQ=255)
                                                 aGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAg                                                                                                 >  2:2530155/1‑150 (MQ=255)
                                                  ggAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGc                                                                                                >  2:1423035/1‑150 (MQ=255)
                                                   gAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCa                                                                                               <  1:1206092/150‑1 (MQ=255)
                                                   gAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCa                                                                                               <  2:4054666/150‑1 (MQ=255)
                                                      aCATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCAttt                                                                                            <  1:671998/150‑1 (MQ=255)
                                                       cATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTg                                                                                           >  2:3307312/1‑150 (MQ=255)
                                                         tGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGcc                                                                                         <  2:319081/150‑1 (MQ=255)
                                                                cATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGa                                                                                  >  1:4079613/1‑150 (MQ=255)
                                                                  tGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGc                                                                                >  1:4453152/1‑150 (MQ=255)
                                                                  tGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGc                                                                                <  1:1481732/150‑1 (MQ=255)
                                                                    aTACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCAt                                                                              >  2:2368121/1‑150 (MQ=255)
                                                                    aTACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCAt                                                                              >  1:3261458/1‑150 (MQ=255)
                                                                    aTACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCAt                                                                              >  2:4747394/1‑150 (MQ=255)
                                                                     tACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATa                                                                             <  1:1574166/150‑1 (MQ=255)
                                                                      aCATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAg                                                                            >  2:1894589/1‑150 (MQ=255)
                                                                       cATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGGTCAGCATTTGCCGGATCGAGCATAGt                                                                           >  2:3086924/1‑150 (MQ=255)
                                                                        aTCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGtc                                                                          <  2:708513/150‑1 (MQ=255)
                                                                         tctcCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGtca                                                                         <  2:1648937/150‑1 (MQ=255)
                                                                         tctcCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGtca                                                                         <  2:1752317/150‑1 (MQ=255)
                                                                          ctcCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGtcat                                                                        <  2:54687/150‑1 (MQ=255)
                                                                          ctcCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGtcat                                                                        >  2:1737728/1‑150 (MQ=255)
                                                                           tcCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATACAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGtcatc                                                                       <  2:1513357/150‑1 (MQ=255)
                                                                              aTAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCatcatc                                                                    >  1:4422622/1‑150 (MQ=255)
                                                                               tAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCt                                                                   >  1:4510131/1‑150 (MQ=255)
                                                                               tAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCt                                                                   >  2:1990299/1‑150 (MQ=255)
                                                                                aaTGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCtt                                                                  >  2:569720/1‑150 (MQ=255)
                                                                                 aTGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTg                                                                 >  2:4363819/1‑150 (MQ=255)
                                                                                   gTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTa                                                               >  2:4629139/1‑150 (MQ=255)
                                                                                     cGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACt                                                             <  2:2653005/150‑1 (MQ=255)
                                                                                     cGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACt                                                             <  2:1324593/150‑1 (MQ=255)
                                                                                     cGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACt                                                             <  2:3997424/150‑1 (MQ=255)
                                                                                      gatgatAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATACTTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTTCTg                                                            >  1:1746501/1‑150 (MQ=255)
                                                                                      gatgatAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTg                                                            >  2:1194233/1‑150 (MQ=255)
                                                                                      gatgatAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTg                                                            >  1:835812/1‑150 (MQ=255)
                                                                                      gatgatAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTg                                                            >  2:46665/1‑150 (MQ=255)
                                                                                           tAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAg                                                       <  1:2421967/150‑1 (MQ=255)
                                                                                              tGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACa                                                    <  1:3357686/150‑1 (MQ=255)
                                                                                               gCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACaa                                                   <  2:4014793/150‑1 (MQ=255)
                                                                                               gCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACaa                                                   <  2:926117/150‑1 (MQ=255)
                                                                                                 ggCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACaaaa                                                 >  2:4310218/1‑150 (MQ=255)
                                                                                                  gCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAAt                                                >  1:1027761/1‑150 (MQ=255)
                                                                                                     ggTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAg                                             >  2:869828/1‑150 (MQ=255)
                                                                                                       tttCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGAc                                           <  1:2816899/150‑1 (MQ=255)
                                                                                                          cttcttTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGtt                                        >  1:1248758/1‑150 (MQ=255)
                                                                                                           ttcttTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTg                                       >  1:3580665/1‑150 (MQ=255)
                                                                                                                        ccAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCAtt                          <  2:1440993/150‑1 (MQ=255)
                                                                                                                         cAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCAttt                         >  2:2732576/1‑150 (MQ=255)
                                                                                                                         cAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCAttt                         >  1:932356/1‑150 (MQ=255)
                                                                                                                          aaaCTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTc                        >  1:391253/1‑150 (MQ=255)
                                                                                                                                  gATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAttgttg                >  1:3479210/1‑150 (MQ=255)
                                                                                                                                  gATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAttgttg                >  1:2682343/1‑150 (MQ=255)
                                                                                                                                  gATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAttgttg                >  1:261823/1‑150 (MQ=255)
                                                                                                                                   aTCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGa               >  1:3062005/1‑150 (MQ=255)
                                                                                                                                   aTCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGa               >  2:1149721/1‑150 (MQ=255)
                                                                                                                                   aTCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGa               <  2:3015496/150‑1 (MQ=255)
                                                                                                                                     ccAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACt             <  2:3762846/150‑1 (MQ=255)
                                                                                                                                        aTTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCt          <  2:4556623/150‑1 (MQ=255)
                                                                                                                                          tCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGt        >  1:3319741/1‑150 (MQ=255)
                                                                                                                                            aCGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTAt      <  2:2555578/150‑1 (MQ=255)
                                                                                                                                               cATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCAc   >  2:734022/1‑150 (MQ=255)
                                                                                                                                                aTCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCAct  <  1:2034445/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAGGTGACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACT  >  CP009273/2405799‑2406092

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: ‑