Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 847,118 G→T 15.8% I312I (ATC→ATA ybiP ← putative EptAB family phosphoethanolamine transferase, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273847,1180GT15.8% 924.0 / 103.0 310I312I (ATC→ATAybiPputative EptAB family phosphoethanolamine transferase, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (134/127);  new base T (23/26);  total (157/153)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 312 total aligned reads displayed.

GTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGT  >  CP009273/846970‑847267
                                                                                                                                                    |                                                                                                                                                     
  aGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATa                                                                                                                                                    >  2:1505891/1‑150 (MQ=255)
   gAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAg                                                                                                                                                   <  1:3392202/150‑1 (MQ=255)
    aaCAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGc                                                                                                                                                  >  2:490033/1‑150 (MQ=255)
     aCAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCg                                                                                                                                                 >  2:1086387/1‑150 (MQ=255)
     aCAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCg                                                                                                                                                 <  1:3793623/150‑1 (MQ=255)
       aaTCAGCTGCGGTTGCGCGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCTTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGt                                                                                                                                               >  1:3585809/1‑150 (MQ=255)
       aaTCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGt                                                                                                                                               <  2:187470/150‑1 (MQ=255)
       aaTCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGt                                                                                                                                               <  2:2873657/150‑1 (MQ=255)
        aTCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTa                                                                                                                                              <  2:3545087/150‑1 (MQ=255)
           aGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCg                                                                                                                                           <  1:1064740/150‑1 (MQ=255)
                cGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTc                                                                                                                                      >  2:704497/1‑150 (MQ=255)
                 ggTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCg                                                                                                                                     <  2:2885649/150‑1 (MQ=255)
                   ttGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGcc                                                                                                                                   <  2:3827510/150‑1 (MQ=255)
                   ttGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGcc                                                                                                                                   <  1:867102/150‑1 (MQ=255)
                      cGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGAt                                                                                                                                >  1:2868399/1‑150 (MQ=255)
                      cGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGAt                                                                                                                                >  2:4171915/1‑150 (MQ=255)
                       gAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGAtt                                                                                                                               <  2:4328632/150‑1 (MQ=255)
                        aGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGAttt                                                                                                                              >  2:4536917/1‑150 (MQ=255)
                        aGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGAttt                                                                                                                              >  2:4098776/1‑150 (MQ=255)
                         gTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTg                                                                                                                             >  2:577971/1‑150 (MQ=255)
                         gTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTg                                                                                                                             >  1:314960/1‑150 (MQ=255)
                            ctctTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGAcc                                                                                                                          <  1:2043542/150‑1 (MQ=255)
                             tctTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGAccc                                                                                                                         <  1:5313403/150‑1 (MQ=255)
                             tctTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGAccc                                                                                                                         >  2:3879236/1‑150 (MQ=255)
                               ttGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTg                                                                                                                       >  1:1947306/1‑150 (MQ=255)
                                   gcCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACACTGGTTg                                                                                                                   >  2:2767353/1‑150 (MQ=255)
                                    ccAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCGTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTgg                                                                                                                  <  2:889310/150‑1 (MQ=255)
                                     cAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGa                                                                                                                 >  2:5137252/1‑150 (MQ=255)
                                     cAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGa                                                                                                                 >  2:808628/1‑150 (MQ=255)
                                          cTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGAACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAacc                                                                                                            >  1:1334055/1‑150 (MQ=255)
                                             gAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAaccacc                                                                                                         >  1:4601954/1‑150 (MQ=255)
                                             gAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAaccacc                                                                                                         <  1:2448028/150‑1 (MQ=255)
                                             gAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATATGACCCTGGTTGGAAAaccacc                                                                                                         <  1:2159436/150‑1 (MQ=255)
                                                cGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAgg                                                                                                      <  2:605687/150‑1 (MQ=255)
                                                cGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAgg                                                                                                      <  1:2813734/150‑1 (MQ=255)
                                                 ggTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGt                                                                                                     >  2:3430170/1‑150 (MQ=255)
                                                  gTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTc                                                                                                    >  1:2236621/1‑150 (MQ=255)
                                                     atatCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTTCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTgg                                                                                                 <  1:4101628/150‑1 (MQ=255)
                                                     atatCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTgg                                                                                                 <  2:3052249/150‑1 (MQ=255)
                                                     atatCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTgg                                                                                                 >  1:842497/1‑150 (MQ=255)
                                                         ccAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGc                                                                                             >  1:4681538/1‑150 (MQ=255)
                                                         ccAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGc                                                                                             >  1:2660897/1‑150 (MQ=255)
                                                          cAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGcc                                                                                            <  1:1853966/150‑1 (MQ=255)
                                                             tAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCcgc                                                                                         <  1:3673355/150‑1 (MQ=255)
                                                             tAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCcgc                                                                                         >  1:2807288/1‑150 (MQ=255)
                                                                  cTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGAt                                                                                    <  1:1139840/150‑1 (MQ=255)
                                                                   ttCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGAtt                                                                                   <  2:4747999/150‑1 (MQ=255)
                                                                   ttCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGAtt                                                                                   >  2:2817833/1‑150 (MQ=255)
                                                                        cTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCa                                                                              <  1:2257607/150‑1 (MQ=255)
                                                                         tttGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAg                                                                             >  1:293352/1‑150 (MQ=255)
                                                                          ttGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAgg                                                                            <  2:5278333/150‑1 (MQ=255)
                                                                              tgtTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGa                                                                        >  2:4938764/1‑150 (MQ=255)
                                                                              tgtTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGa                                                                        <  2:2466484/150‑1 (MQ=255)
                                                                              tgtTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGa                                                                        >  1:670302/1‑150 (MQ=255)
                                                                               gtTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGAc                                                                       <  1:1505891/150‑1 (MQ=255)
                                                                                tttttATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACa                                                                      <  1:2491245/150‑1 (MQ=255)
                                                                                ttattATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACa                                                                      <  1:5136267/147‑1 (MQ=255)
                                                                                      tCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTa                                                                <  2:2979334/150‑1 (MQ=255)
                                                                                      tCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTa                                                                <  2:2629650/150‑1 (MQ=255)
                                                                                        tGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATc                                                              >  1:3941768/1‑150 (MQ=255)
                                                                                          cTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCt                                                            <  2:336046/150‑1 (MQ=255)
                                                                                           ttCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTg                                                           >  2:839087/1‑150 (MQ=255)
                                                                                           ttCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTg                                                           <  2:4646502/150‑1 (MQ=255)
                                                                                               aaaTTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAAc                                                       >  2:2009527/1‑150 (MQ=255)
                                                                                                aaTTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACt                                                      >  1:4354041/1‑150 (MQ=255)
                                                                                                 aTTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACtg                                                     >  1:583584/1‑150 (MQ=255)
                                                                                                 aTTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACtg                                                     >  1:3997210/1‑150 (MQ=255)
                                                                                                 aTTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACtg                                                     >  1:4438282/1‑150 (MQ=255)
                                                                                                   tACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACtgtg                                                   <  2:161014/150‑1 (MQ=255)
                                                                                                    aCCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTgg                                                  <  2:1914439/150‑1 (MQ=255)
                                                                                                     ccTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGt                                                 <  2:1830399/150‑1 (MQ=255)
                                                                                                        tcttTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTg                                              <  2:5180840/150‑1 (MQ=255)
                                                                                                          tttCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGcc                                            >  1:3837296/1‑150 (MQ=255)
                                                                                                             cAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATc                                         <  2:2867731/150‑1 (MQ=255)
                                                                                                              aGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCg                                        >  1:4932315/1‑150 (MQ=255)
                                                                                                               ggAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGa                                       >  2:1411459/1‑150 (MQ=255)
                                                                                                               ggAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGa                                       >  2:2665451/1‑150 (MQ=255)
                                                                                                               ggAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGa                                       >  2:1701692/1‑150 (MQ=255)
                                                                                                                gAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGAc                                      >  1:3530676/1‑150 (MQ=255)
                                                                                                                   gTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAAc                                   <  2:1358844/150‑1 (MQ=255)
                                                                                                                    tACACTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAAcg                                  <  2:2535469/150‑1 (MQ=255)
                                                                                                                      cacTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTAGCGTATTCGCCGAGTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCCGGGTGACACAGTTATCCTGAACCTGTGGTTTGCCATCGACAAcgcg                                >  2:3408317/1‑150 (MQ=255)
                                                                                                                       acTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGa                               >  2:1190799/1‑150 (MQ=255)
                                                                                                                          tCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTg                            <  1:299762/150‑1 (MQ=255)
                                                                                                                           cATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGa                           >  1:4792466/1‑150 (MQ=255)
                                                                                                                            aTCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAg                          >  2:4054631/1‑150 (MQ=255)
                                                                                                                            aTCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAg                          >  2:1826212/1‑150 (MQ=255)
                                                                                                                                gCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGtt                      <  1:263332/150‑1 (MQ=255)
                                                                                                                                 cTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTa                     >  2:4432986/1‑150 (MQ=255)
                                                                                                                                   cGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAg                   <  1:3960557/150‑1 (MQ=255)
                                                                                                                                   cGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAg                   <  1:1838767/150‑1 (MQ=255)
                                                                                                                                      ttGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCa                <  2:2280626/150‑1 (MQ=255)
                                                                                                                                       tGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCaa               <  1:1251392/150‑1 (MQ=255)
                                                                                                                                       tGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCaa               <  2:2164999/150‑1 (MQ=255)
                                                                                                                                            aTGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGAt          >  2:3304500/1‑150 (MQ=255)
                                                                                                                                               cTGGCTATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACACAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGCTTTc       >  2:3738697/1‑150 (MQ=255)
                                                                                                                                                tGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCt      <  1:4171915/150‑1 (MQ=255)
                                                                                                                                                 ggCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTg     <  1:4900165/150‑1 (MQ=255)
                                                                                                                                                   cGATAGCGGTATCGTTTTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCg   <  2:3112928/150‑1 (MQ=255)
                                                                                                                                                    gATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGt  >  1:2584578/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGTACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGT  >  CP009273/846970‑847267

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: