Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,991,234 T→A 7.5% L231* (TTG→TAG)  yigB → 5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase; pyrimidine phosphatase; riboflavin synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,991,2340TA 7.5% 813.3 / 42.8 240L231* (TTG→TAG) yigB5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase; pyrimidine phosphatase; riboflavin synthesis
Reads supporting (aligned to +/- strand):  ref base T (100/122);  new base A (11/7);  total (111/129)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 243 total aligned reads displayed.

GAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACA  >  CP009273/3991086‑3991383
                                                                                                                                                    |                                                                                                                                                     
gagaTCTTACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTg                                                                                                                                                      >  1:4040804/1‑150 (MQ=255)
gagaTCTTACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTg                                                                                                                                                      <  2:4630643/150‑1 (MQ=255)
  gaTCTTACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGTCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGc                                                                                                                                                    >  1:2421699/1‑150 (MQ=255)
   aTCTTACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCa                                                                                                                                                   <  1:3633956/150‑1 (MQ=255)
   aTCTTACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCa                                                                                                                                                   <  1:303468/150‑1 (MQ=255)
   aTCTGACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCa                                                                                                                                                   <  1:517381/150‑1 (MQ=255)
      ttACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCAtct                                                                                                                                                <  1:1865617/150‑1 (MQ=255)
        aCATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATctct                                                                                                                                              <  1:2652741/150‑1 (MQ=255)
        aCATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATctct                                                                                                                                              <  1:3419842/150‑1 (MQ=255)
              tGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTc                                                                                                                                        <  2:1116874/150‑1 (MQ=255)
              tGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTc                                                                                                                                        >  1:4992872/1‑150 (MQ=255)
              tGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTc                                                                                                                                        >  1:1339018/1‑150 (MQ=255)
               ggggACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCCTCTGGAAATTTCCCGGTTGGCATCTCTGACCACg                                                                                                                                       >  2:3299108/1‑150 (MQ=255)
               ggggACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCg                                                                                                                                       <  2:3210323/150‑1 (MQ=255)
                  gACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTg                                                                                                                                    <  2:454265/150‑1 (MQ=255)
                   aCGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGa                                                                                                                                   <  1:4682922/150‑1 (MQ=255)
                    cGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGat                                                                                                                                  <  2:795986/150‑1 (MQ=255)
                     gATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGata                                                                                                                                 >  1:485122/1‑150 (MQ=255)
                      aTCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGatat                                                                                                                                >  2:3829241/1‑150 (MQ=255)
                            ccacTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAg                                                                                                                          <  1:1939878/150‑1 (MQ=255)
                             cacTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGc                                                                                                                         >  1:454245/1‑150 (MQ=255)
                                tGACGTGGGTGGGGCAATTCGCAGCGGAAGGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCaaa                                                                                                                      <  2:1883400/150‑1 (MQ=255)
                                   cGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCt                                                                                                                   <  2:3184603/150‑1 (MQ=255)
                                    gtgggtggGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTg                                                                                                                  <  1:740352/150‑1 (MQ=255)
                                      gggtgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGta                                                                                                                >  1:5040677/1‑150 (MQ=255)
                                       ggtgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGtat                                                                                                               <  1:4140501/150‑1 (MQ=255)
                                         tgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGtatat                                                                                                             >  2:4307840/1‑150 (MQ=255)
                                               aagtCGCAGCGGAATGCATGCTTGTTGGATCAGACCTGAAAATGTCGATCTGATGCATACCTGGGACATCCGTTTACTGCCGAATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATAccc                                                                                                       <  2:3999024/147‑1 (MQ=255)
                                               aaTTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAACTCTGTATATATAccc                                                                                                       >  1:973430/1‑150 (MQ=255)
                                               aaTTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATAccc                                                                                                       >  1:3072919/1‑150 (MQ=255)
                                                aTTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCa                                                                                                      <  2:516256/150‑1 (MQ=255)
                                                aTTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCa                                                                                                      <  1:385296/150‑1 (MQ=255)
                                                 ttCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAg                                                                                                     <  2:905624/150‑1 (MQ=255)
                                                  tCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGc                                                                                                    >  2:2289481/1‑150 (MQ=255)
                                                   cGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCt                                                                                                   <  2:1334854/150‑1 (MQ=255)
                                                   cGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCt                                                                                                   >  1:3159255/1‑150 (MQ=255)
                                                    gCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCtt                                                                                                  >  1:3272089/1‑150 (MQ=255)
                                                       gCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCttttt                                                                                               <  2:3062862/150‑1 (MQ=255)
                                                        cGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTg                                                                                              >  1:90453/1‑150 (MQ=255)
                                                        cGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTg                                                                                              >  1:1737578/1‑150 (MQ=255)
                                                          gAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGc                                                                                            >  2:1457084/1‑150 (MQ=255)
                                                             tGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGa                                                                                         <  1:3985086/150‑1 (MQ=255)
                                                              gCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAg                                                                                        <  2:3888907/150‑1 (MQ=255)
                                                              gCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAg                                                                                        >  1:3786198/1‑150 (MQ=255)
                                                                aGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAggg                                                                                      >  1:2296324/1‑150 (MQ=255)
                                                                 ggCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGc                                                                                     >  2:1408624/1‑150 (MQ=255)
                                                                   cTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGt                                                                                   >  2:28817/1‑150 (MQ=255)
                                                                   cTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCAGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGt                                                                                   >  2:28361/1‑150 (MQ=255)
                                                                    ttgttgGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGtt                                                                                  >  2:1297469/1‑150 (MQ=255)
                                                                     tgttgGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTg                                                                                 <  1:1169249/150‑1 (MQ=255)
                                                                         ggATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCt                                                                             >  2:4226121/1‑150 (MQ=255)
                                                                         ggATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTAGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCt                                                                             >  1:2196515/1‑150 (MQ=255)
                                                                           aTCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTtc                                                                           <  1:489373/150‑1 (MQ=255)
                                                                              aGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTcc                                                                        <  2:4950454/150‑1 (MQ=255)
                                                                              aGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTcc                                                                        <  1:2521046/150‑1 (MQ=255)
                                                                                 ccGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGcc                                                                     >  1:2828400/1‑150 (MQ=255)
                                                                                  cGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGccc                                                                    <  2:3453144/150‑1 (MQ=255)
                                                                                   ggAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCCGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCa                                                                   >  1:1373433/1‑150 (MQ=255)
                                                                                         tGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTAt                                                             <  1:2465314/150‑1 (MQ=255)
                                                                                         tGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTAGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTAt                                                             >  1:4174062/1‑150 (MQ=255)
                                                                                            cGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTAtttt                                                          >  2:2238285/1‑150 (MQ=255)
                                                                                               tCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTAc                                                       >  1:5054360/1‑150 (MQ=255)
                                                                                               tCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTAc                                                       <  1:3393846/150‑1 (MQ=255)
                                                                                                cTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTAcg                                                      >  1:1619748/1‑150 (MQ=255)
                                                                                                  gATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTAGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTAcgcg                                                    >  2:2031848/1‑150 (MQ=255)
                                                                                                      cAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGcggcgg                                                <  2:2742346/150‑1 (MQ=255)
                                                                                                       aaaCCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGt                                               <  1:3234369/150‑1 (MQ=255)
                                                                                                         aCCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGc                                             <  2:3028126/150‑1 (MQ=255)
                                                                                                         aCCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTAGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGc                                             <  2:1797248/150‑1 (MQ=255)
                                                                                                          ccTGGGACAGCCGTTTAGTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGcc                                            <  1:3927701/150‑1 (MQ=255)
                                                                                                           cTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCa                                           <  1:628131/150‑1 (MQ=255)
                                                                                                           cTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCa                                           >  1:4340629/1‑150 (MQ=255)
                                                                                                            tGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCaa                                          >  2:4568969/1‑150 (MQ=255)
                                                                                                            tGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCaa                                          >  2:3850203/1‑150 (MQ=255)
                                                                                                             gggACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAAt                                         <  1:1271524/150‑1 (MQ=255)
                                                                                                             gggACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAAt                                         <  2:1189053/150‑1 (MQ=255)
                                                                                                              ggACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATg                                        <  1:274221/150‑1 (MQ=255)
                                                                                                              ggACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATg                                        <  2:352947/150‑1 (MQ=255)
                                                                                                                 cAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGAc                                     >  1:1471822/1‑150 (MQ=255)
                                                                                                                   gCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGt                                   >  2:1733430/1‑150 (MQ=255)
                                                                                                                      gTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTc                                <  2:2454435/150‑1 (MQ=255)
                                                                                                                      gTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTc                                >  1:2007291/1‑150 (MQ=255)
                                                                                                                        ttACTGCCGCATCTGGAAATTTCCCGGTAGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCtt                              <  2:4982340/150‑1 (MQ=255)
                                                                                                                        ttACTGCCGCATCTGGAAATGTCCCGGTAGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCtt                              <  2:503807/150‑1 (MQ=255)
                                                                                                                         tACTGCCGCATCTGGAAATTTCCCGGTAGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTa                             >  2:2256449/1‑150 (MQ=255)
                                                                                                                           cTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTAcc                           >  2:2701969/1‑150 (MQ=255)
                                                                                                                            tGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCt                          >  1:2106143/1‑150 (MQ=255)
                                                                                                                             gCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTg                         >  1:207262/1‑150 (MQ=255)
                                                                                                                                gCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTc                      <  1:346509/150‑1 (MQ=255)
                                                                                                                                gCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTc                      <  2:1950536/150‑1 (MQ=255)
                                                                                                                                  aTCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGa                    >  2:3625274/1‑150 (MQ=255)
                                                                                                                                    cTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACa                  >  2:4157671/1‑150 (MQ=255)
                                                                                                                                        aaaTTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCt              <  1:4667083/150‑1 (MQ=255)
                                                                                                                                           tttCCCGGTAGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTaa           >  2:361984/1‑150 (MQ=255)
                                                                                                                                             tCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATAGACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATg         <  1:494071/150‑1 (MQ=255)
                                                                                                                                              cccGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGa        <  1:3386720/150‑1 (MQ=255)
                                                                                                                                              cccGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGa        >  1:840718/1‑150 (MQ=255)
                                                                                                                                                cGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACa      <  1:2586243/150‑1 (MQ=255)
                                                                                                                                                 ggTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACaa     >  2:1695169/1‑150 (MQ=255)
                                                                                                                                                   ttGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAAc   >  2:278598/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACA  >  CP009273/3991086‑3991383

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: