Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 2,931,717 C→A N307K (AAC→AAA fucK → L‑fuculokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,931,7170CA99.6% 772.8 / NA 232N307K (AAC→AAAfucKL‑fuculokinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (124/107);  minor base T (0/1);  total (124/108)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 232 total aligned reads displayed.

GCGCTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATT  >  CP009273/2931572‑2931863
                                                                                                                                                 |                                                                                                                                                  
gcgcTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAg                                                                                                                                                <  2:1067567/150‑1 (MQ=255)
gcgcTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAg                                                                                                                                                <  2:1707556/150‑1 (MQ=255)
 cgcTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAgg                                                                                                                                               >  1:1026568/1‑150 (MQ=255)
     ggTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATg                                                                                                                                           <  2:3723940/150‑1 (MQ=255)
       tGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCa                                                                                                                                         <  2:1768886/150‑1 (MQ=255)
       tGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCa                                                                                                                                         >  2:3114005/1‑150 (MQ=255)
             aCAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAAtggct                                                                                                                                   <  1:46296/150‑1 (MQ=255)
              cAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAAtggctg                                                                                                                                  <  1:916954/150‑1 (MQ=255)
               aaaaTGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAAtggctgg                                                                                                                                 <  1:602884/150‑1 (MQ=255)
               aaaaTGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAAtggctgg                                                                                                                                 <  2:239868/150‑1 (MQ=255)
                 aaTGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCa                                                                                                                               >  1:2307937/1‑150 (MQ=255)
                 aaTGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCACAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCa                                                                                                                               >  1:1861879/1‑150 (MQ=255)
                    gAACCCGTGCTCTCTTCCGGTACATGGTAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTACGCTGTTAAGTCAGTACGCCGGTTCCACCTCCGAACTGGATAGCCAGGCAGAATTGTATAAACCAAGAATGCAATGGCAGGCATca                                                                                                                            >  1:1243805/1‑149 (MQ=255)
                     aaCCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCACTGGCTGGCATCCg                                                                                                                           >  2:358535/1‑150 (MQ=255)
                     aaCCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCg                                                                                                                           >  2:1102730/1‑150 (MQ=255)
                     aaCCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCg                                                                                                                           >  2:1307206/1‑150 (MQ=255)
                      aCCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCgg                                                                                                                          >  2:2620169/1‑150 (MQ=255)
                         tgTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGt                                                                                                                       >  2:2979886/2‑150 (MQ=255)
                          gTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTg                                                                                                                      >  2:2314207/1‑150 (MQ=255)
                          gTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTg                                                                                                                      <  1:1182897/150‑1 (MQ=255)
                             ctctctTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTg                                                                                                                   <  2:463499/150‑1 (MQ=255)
                             ctctctTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTg                                                                                                                   >  1:3400108/1‑150 (MQ=255)
                              tctctTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTgg                                                                                                                  >  1:3359241/1‑150 (MQ=255)
                               ctctTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACTCCGGTTCCACCTGCGACCTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGa                                                                                                                 >  1:254073/1‑150 (MQ=255)
                               ctctTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGa                                                                                                                 >  1:1692038/1‑150 (MQ=255)
                                tctTCCGGTACATTGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTAAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATCGCTGGCATCCGGCGTGCTgtca                                                                                                                >  2:2251913/1‑147 (MQ=255)
                                tctTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGaa                                                                                                                <  1:2232580/150‑1 (MQ=255)
                                        tACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTgag                                                                                                        >  2:2171873/1‑150 (MQ=255)
                                        tACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTgag                                                                                                        >  1:2095541/1‑150 (MQ=255)
                                        tACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTgag                                                                                                        <  2:686895/150‑1 (MQ=255)
                                         aCATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTgaga                                                                                                       <  2:2823013/150‑1 (MQ=255)
                                          cATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGaa                                                                                                      >  2:2542132/1‑150 (MQ=255)
                                           aTGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGaaa                                                                                                     <  2:472321/150‑1 (MQ=255)
                                             gggAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAAc                                                                                                   <  1:1096675/150‑1 (MQ=255)
                                             gggAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAAc                                                                                                   <  2:1605336/150‑1 (MQ=255)
                                                   ttttAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCt                                                                                             >  2:817342/1‑150 (MQ=255)
                                                   ttttAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCCGGAATGGGTGAGAAAACTGTTCt                                                                                             <  2:2702502/150‑1 (MQ=255)
                                                      tAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGa                                                                                          <  1:1165405/150‑1 (MQ=255)
                                                      tAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGa                                                                                          >  2:803298/1‑150 (MQ=255)
                                                       aaTGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGAc                                                                                         <  1:443130/150‑1 (MQ=255)
                                                        aTGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACg                                                                                        >  2:2651082/1‑150 (MQ=255)
                                                           gTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCt                                                                                     <  2:207904/150‑1 (MQ=255)
                                                            ttCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTg                                                                                    <  1:1608746/150‑1 (MQ=255)
                                                              cGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGaa                                                                                  <  2:1968262/150‑1 (MQ=255)
                                                               gCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGaaa                                                                                 <  1:3560128/150‑1 (MQ=255)
                                                                cAGCGCCCAGGTTGATGCTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAac                                                                                <  2:559929/150‑1 (MQ=255)
                                                                cAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAac                                                                                <  2:821072/150‑1 (MQ=255)
                                                                       cAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTgg                                                                         >  1:2186605/1‑150 (MQ=255)
                                                                        aGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGc                                                                        >  2:1243610/1‑150 (MQ=255)
                                                                        aGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGc                                                                        >  2:1380924/1‑150 (MQ=255)
                                                                         ggTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCa                                                                       <  2:1175271/150‑1 (MQ=255)
                                                                          gTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCaa                                                                      >  1:408426/1‑150 (MQ=255)
                                                                             gATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATg                                                                   <  2:3359062/150‑1 (MQ=255)
                                                                              aTACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGt                                                                  >  1:3109619/1‑150 (MQ=255)
                                                                              aTACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGt                                                                  <  1:1306573/150‑1 (MQ=255)
                                                                                 cTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGttga                                                               >  1:1065352/1‑150 (MQ=255)
                                                                                   tCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGttgatt                                                             >  1:3416205/1‑150 (MQ=255)
                                                                                   tCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGttgatt                                                             >  2:1887776/1‑150 (MQ=255)
                                                                                    cGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGttgattg                                                            >  1:854813/1‑150 (MQ=255)
                                                                                     gCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGttgattga                                                           >  2:3694040/1‑150 (MQ=255)
                                                                                      cTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTgaa                                                          >  2:1354801/1‑150 (MQ=255)
                                                                                       tGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTgaag                                                         >  2:2664523/1‑150 (MQ=255)
                                                                                          tAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGaagaag                                                      <  1:584085/150‑1 (MQ=255)
                                                                                            aGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCACATGTTGATTGAAGAAGCt                                                    >  2:2122465/1‑150 (MQ=255)
                                                                                            aGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCt                                                    <  2:2373913/150‑1 (MQ=255)
                                                                                             gTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTc                                                   <  1:1030293/150‑1 (MQ=255)
                                                                                                aGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGTTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTc                                                >  1:820806/1‑150 (MQ=255)
                                                                                                 gTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCt                                               <  2:3492838/150‑1 (MQ=255)
                                                                                                 gTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCt                                               <  2:2204020/150‑1 (MQ=255)
                                                                                                    cGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGAt                                            >  1:942144/1‑150 (MQ=255)
                                                                                                     gCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATc                                           <  1:634524/150‑1 (MQ=255)
                                                                                                     gCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATc                                           >  2:3240755/1‑150 (MQ=255)
                                                                                                     gCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATc                                           <  1:2423614/150‑1 (MQ=255)
                                                                                                         gTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCgcgc                                       >  1:971953/1‑150 (MQ=255)
                                                                                                              aCCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGgc                                  <  1:2072649/150‑1 (MQ=255)
                                                                                                               ccTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGgcg                                 <  2:1542169/150‑1 (MQ=255)
                                                                                                                cTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGgcgc                                <  1:1813021/150‑1 (MQ=255)
                                                                                                                cTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGgcgc                                <  1:3707112/150‑1 (MQ=255)
                                                                                                                   cGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGa                             >  2:1042933/1‑150 (MQ=255)
                                                                                                                    gAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGAt                            <  1:144039/150‑1 (MQ=255)
                                                                                                                     aaCTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATg                           >  2:1831421/1‑150 (MQ=255)
                                                                                                                      aCTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATgg                          >  1:2834244/1‑150 (MQ=255)
                                                                                                                           aTAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTaa                     <  2:3078876/150‑1 (MQ=255)
                                                                                                                               ccAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAAt                 <  2:3400108/150‑1 (MQ=255)
                                                                                                                                caggcaggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATg                >  1:1325996/1‑150 (MQ=255)
                                                                                                                                  ggcaggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCa              >  1:3266758/1‑150 (MQ=255)
                                                                                                                                    caggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAgt            <  1:3114005/150‑1 (MQ=255)
                                                                                                                                    caggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAgt            >  2:2007331/1‑150 (MQ=255)
                                                                                                                                     aggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAgtg           >  1:3467443/1‑150 (MQ=255)
                                                                                                                                     aggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAgtg           <  1:2653846/150‑1 (MQ=255)
                                                                                                                                      gggTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAgtgt          >  2:3426967/1‑150 (MQ=255)
                                                                                                                                      gggTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAgtgt          >  1:2870534/1‑150 (MQ=255)
                                                                                                                                        gTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGa        >  2:1403875/1‑150 (MQ=255)
                                                                                                                                           gTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGAttt     <  2:567591/150‑1 (MQ=255)
                                                                                                                                            tataAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTa    >  1:2472410/1‑150 (MQ=255)
                                                                                                                                            tataAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTa    >  1:1484595/1‑150 (MQ=255)
                                                                                                                                             ataAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTAt   >  2:356319/1‑150 (MQ=255)
                                                                                                                                             ataAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTAt   <  1:1654443/150‑1 (MQ=255)
                                                                                                                                              taAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTAtt  <  1:2314227/150‑1 (MQ=255)
                                                                                                                                              taAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTAtt  >  1:2609969/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCGCTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATT  >  CP009273/2931572‑2931863

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: