Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,394,875 G→T 91.4% A124A (GCG→GCT purA → adenylosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,394,8750GT91.4% 936.4 / 79.2 325A124A (GCG→GCTpurAadenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base G (19/9);  new base T (135/162);  total (154/171)
Fisher's exact test for biased strand distribution p-value = 2.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 335 total aligned reads displayed.

GTCTCCGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTCCAGTTGG  >  CP009273/4394728‑4395023
                                                                                                                                                   |                                                                                                                                                    
  ctcCGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAg                                                                                                                                                  >  1:2214601/1‑150 (MQ=255)
   tcCGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGc                                                                                                                                                 <  1:70932/150‑1 (MQ=255)
   tcCGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGc                                                                                                                                                 >  1:3970228/1‑150 (MQ=255)
       gCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATc                                                                                                                                             <  2:3802225/150‑1 (MQ=255)
       gCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATc                                                                                                                                             <  1:1625926/150‑1 (MQ=255)
         cgcgCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCgg                                                                                                                                           >  2:2254822/1‑150 (MQ=255)
         cgcgCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCgg                                                                                                                                           >  2:4224541/1‑150 (MQ=255)
          gcgcTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGc                                                                                                                                          >  1:2802408/1‑150 (MQ=255)
           cgcTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGca                                                                                                                                         <  1:3260487/150‑1 (MQ=255)
             cTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGcacc                                                                                                                                       >  2:1732546/1‑150 (MQ=255)
              tgatgaAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGcacca                                                                                                                                      >  1:4195175/1‑150 (MQ=255)
                 tgaaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCg                                                                                                                                   <  1:4024217/150‑1 (MQ=255)
                 tgaaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCg                                                                                                                                   <  1:3502988/150‑1 (MQ=255)
                    aagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTc                                                                                                                                <  2:500270/150‑1 (MQ=255)
                    aagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTc                                                                                                                                <  2:520597/150‑1 (MQ=255)
                    aagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTc                                                                                                                                <  2:2277200/150‑1 (MQ=255)
                       ataagaaagaACTGGAAGACCGTGGCATCCACGTACGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCTTGGCGCTAAAGCGATCGGCACCACCGGTCGTg                                                                                                                             >  2:851057/5‑150 (MQ=255)
                       agatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTg                                                                                                                             >  2:1412741/1‑150 (MQ=255)
                       agatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTg                                                                                                                             <  2:392870/150‑1 (MQ=255)
                         atgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGt                                                                                                                           >  1:2522755/1‑150 (MQ=255)
                          tgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTa                                                                                                                          <  1:4121084/150‑1 (MQ=255)
                            aaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATc                                                                                                                        <  2:4593326/150‑1 (MQ=255)
                                aaCTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGc                                                                                                                    <  2:24352/150‑1 (MQ=255)
                                   tGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTg                                                                                                                 <  2:810617/150‑1 (MQ=255)
                                     gAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCt                                                                                                               <  2:1970981/150‑1 (MQ=255)
                                      aaGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCtt                                                                                                              <  1:766726/150‑1 (MQ=255)
                                      aaGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCtt                                                                                                              >  2:2847946/1‑150 (MQ=255)
                                       aGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTa                                                                                                             >  1:1953419/1‑150 (MQ=255)
                                       aGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTa                                                                                                             <  2:2813966/150‑1 (MQ=255)
                                         aCCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATg                                                                                                           <  1:1889497/150‑1 (MQ=255)
                                          ccGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGa                                                                                                          <  1:3699416/150‑1 (MQ=255)
                                           cGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGaa                                                                                                         <  2:2877068/150‑1 (MQ=255)
                                             tGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGa                                                                                                       <  2:2479239/150‑1 (MQ=255)
                                                cATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATaa                                                                                                    >  1:1849891/1‑150 (MQ=255)
                                                cATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATaa                                                                                                    >  1:926493/1‑150 (MQ=255)
                                                 aTCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATaaa                                                                                                   >  2:217648/1‑150 (MQ=255)
                                                  tCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAg                                                                                                  >  1:1817260/1‑150 (MQ=255)
                                                   ccccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGt                                                                                                 <  1:1728097/150‑1 (MQ=255)
                                                   ccccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGt                                                                                                 >  2:2375630/1‑150 (MQ=255)
                                                   ccccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGt                                                                                                 >  2:1474572/1‑150 (MQ=255)
                                                    cccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTa                                                                                                <  2:3814147/150‑1 (MQ=255)
                                                    cccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTa                                                                                                <  2:274420/150‑1 (MQ=255)
                                                         tCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACg                                                                                           <  1:4627911/150‑1 (MQ=255)
                                                           gTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTc                                                                                         <  2:1137474/150‑1 (MQ=255)
                                                            tGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTcg                                                                                        <  2:3299837/150‑1 (MQ=255)
                                                                 gTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTc                                                                                   <  1:1384422/150‑1 (MQ=255)
                                                                      ctgctgTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCgt                                                                              >  1:2336935/1‑150 (MQ=255)
                                                                         ctgTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGtt                                                                           >  1:3920331/1‑150 (MQ=255)
                                                                           gTCTGAAGCATGTCCGCAGATCCTTGATTCTCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAACGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTgg                                                                         <  1:3464043/150‑1 (MQ=255)
                                                                              tGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGa                                                                      >  2:2650561/1‑150 (MQ=255)
                                                                               gAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGAc                                                                     >  1:145140/1‑150 (MQ=255)
                                                                                aaGCATGTCCGCTGATCCTTTATTATCTCGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGAct                                                                    <  2:3553633/150‑2 (MQ=255)
                                                                                aaGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGAcc                                                                    >  1:831832/1‑150 (MQ=255)
                                                                                 aGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCt                                                                   >  1:3963237/1‑150 (MQ=255)
                                                                                  gCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCtt                                                                  <  2:2916119/150‑1 (MQ=255)
                                                                                    aTGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCtttt                                                                >  2:3500089/1‑150 (MQ=255)
                                                                                    aTGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCtttt                                                                <  1:4353964/150‑1 (MQ=255)
                                                                                      gTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCg                                                              <  2:3714620/150‑1 (MQ=255)
                                                                                       tCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGa                                                             >  2:33455/1‑150 (MQ=255)
                                                                                        ccGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGAc                                                            >  1:447976/1‑150 (MQ=255)
                                                                                         cGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACa                                                           <  1:535226/150‑1 (MQ=255)
                                                                                         cGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACa                                                           <  2:1770640/150‑1 (MQ=255)
                                                                                         cGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACa                                                           <  1:920489/150‑1 (MQ=255)
                                                                                          gCTGATCCTTGATTATCACGTTGCGCTGGATAACACGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAATATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACaa                                                          >  1:1939681/1‑150 (MQ=255)
                                                                                            tGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAg                                                        >  1:4399812/1‑150 (MQ=255)
                                                                                             gATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGa                                                       <  2:1019806/150‑1 (MQ=255)
                                                                                               tCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGaaa                                                     <  2:3850765/150‑1 (MQ=255)
                                                                                               tCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGaaa                                                     >  1:1112882/1‑150 (MQ=255)
                                                                                                     aTTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCg                                               >  2:3112088/1‑150 (MQ=255)
                                                                                                      ttATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGc                                              <  2:4487564/150‑1 (MQ=255)
                                                                                                         tctcGTTGAGCTGGATAACGCGCGTGAGATAGCGCGTGGCGCTAAATCGATCGGCACCACCGGTCGTGGTATCGGGCCTGATTATGATGATAAAGTATCACTTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGTAACCTTCGCTGa                                           <  2:2503612/147‑1 (MQ=255)
                                                                                                         tCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGa                                           <  1:1076662/150‑1 (MQ=255)
                                                                                                          cACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGaa                                          >  2:4054744/1‑150 (MQ=255)
                                                                                                            cGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGaaaa                                        >  2:4130736/1‑150 (MQ=255)
                                                                                                            cGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGaaaa                                        >  1:3084912/1‑150 (MQ=255)
                                                                                                               tGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACt                                     >  2:1907661/1‑150 (MQ=255)
                                                                                                               tGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACt                                     <  1:1131626/150‑1 (MQ=255)
                                                                                                                   cTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGaaa                                 <  2:1278365/150‑1 (MQ=255)
                                                                                                                   cTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGaaa                                 >  1:718197/1‑150 (MQ=255)
                                                                                                                    tGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaag                                >  1:670642/1‑150 (MQ=255)
                                                                                                                     ggATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaaga                               <  1:1197013/150‑1 (MQ=255)
                                                                                                                      gATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaagaa                              <  1:3934306/150‑1 (MQ=255)
                                                                                                                        tAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGt                            >  1:4532806/1‑150 (MQ=255)
                                                                                                                        tAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGt                            <  2:4393042/150‑1 (MQ=255)
                                                                                                                        tAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGt                            >  2:1813332/1‑150 (MQ=255)
                                                                                                                         aaCGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTg                           >  1:2918547/1‑150 (MQ=255)
                                                                                                                           cgcgcgTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGAt                         <  1:2103694/150‑1 (MQ=255)
                                                                                                                           cgcgcgTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGAt                         <  1:3582811/150‑1 (MQ=255)
                                                                                                                               cgTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGaa                     <  2:2945345/150‑1 (MQ=255)
                                                                                                                                gTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAat                    >  1:52943/1‑150 (MQ=255)
                                                                                                                                gTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAat                    <  2:480138/150‑1 (MQ=255)
                                                                                                                                 tGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAata                   <  2:2622570/150‑1 (MQ=255)
                                                                                                                                   agaAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATc                 <  2:4000747/150‑1 (MQ=255)
                                                                                                                                   agaAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATc                 <  2:1676310/150‑1 (MQ=255)
                                                                                                                                    gaAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATca                <  1:3333102/150‑1 (MQ=255)
                                                                                                                                    gaAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATca                >  1:3220605/1‑150 (MQ=255)
                                                                                                                                      aaGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATcaca              >  2:3051025/1‑150 (MQ=255)
                                                                                                                                         cgcgTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACt           <  2:3186354/150‑1 (MQ=255)
                                                                                                                                           cgTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTc         <  1:1032726/150‑1 (MQ=255)
                                                                                                                                                cgctAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTCCAGtt    >  2:4173569/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTCTCCGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTCCAGTTGG  >  CP009273/4394728‑4395023

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: