Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,441,854 C→A 9.6% R27S (CGT→AGT)  feaB → phenylacetaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,441,8540CA 9.6% 984.5 / 51.9 281R27S (CGT→AGT) feaBphenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (140/114);  new base A (13/14);  total (153/128)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 292 total aligned reads displayed.

AGCGAAAAAAGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGG  >  CP009273/1441706‑1442001
                                                                                                                                                    |                                                                                                                                                   
 gCGAAAAAAGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGgtcg                                                                                                                                                    >  2:3567727/1‑149 (MQ=255)
    aaaaaaGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTcc                                                                                                                                                 <  2:1431964/149‑1 (MQ=255)
     aaaaaGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCt                                                                                                                                                <  1:1194653/149‑1 (MQ=255)
     aaaaaGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCt                                                                                                                                                >  1:1983354/1‑149 (MQ=255)
      aaaaGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTg                                                                                                                                               <  1:2093072/149‑1 (MQ=255)
        gaGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGc                                                                                                                                             >  2:40128/2‑149 (MQ=255)
        aaGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGc                                                                                                                                             >  2:1933667/1‑149 (MQ=255)
           tGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGcccc                                                                                                                                          >  2:41915/1‑149 (MQ=255)
           tGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGcccc                                                                                                                                          >  1:1859611/1‑149 (MQ=255)
            gACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCg                                                                                                                                         >  2:3169089/1‑149 (MQ=255)
              cTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGca                                                                                                                                       >  1:31650/1‑149 (MQ=255)
                   cTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAg                                                                                                                                  <  1:349353/149‑1 (MQ=255)
                    ttGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGt                                                                                                                                 <  2:1292205/149‑1 (MQ=255)
                    ttGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGt                                                                                                                                 <  1:2991284/149‑1 (MQ=255)
                     tGTCGCTGCGTACACTTAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTg                                                                                                                                <  2:1603238/149‑1 (MQ=255)
                         gCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGaaaa                                                                                                                            <  1:2944954/149‑1 (MQ=255)
                           tGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACCGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCCCAAAGTGAAAAAc                                                                                                                          >  2:2698537/1‑149 (MQ=255)
                           tGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTCTATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAAc                                                                                                                          >  1:1052699/1‑149 (MQ=255)
                           tGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAAc                                                                                                                          >  2:305878/1‑149 (MQ=255)
                            gCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACg                                                                                                                         <  1:89614/149‑1 (MQ=255)
                                    tGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGa                                                                                                                 >  2:2603399/1‑149 (MQ=255)
                                     gAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGAt                                                                                                                >  2:2040124/1‑149 (MQ=255)
                                      aaaTCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATc                                                                                                               <  2:1663555/149‑1 (MQ=255)
                                       aaTCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCt                                                                                                              >  1:1126398/1‑149 (MQ=255)
                                        aTCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCtt                                                                                                             <  2:1785674/149‑1 (MQ=255)
                                          cacaCTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTg                                                                                                           <  2:2548915/149‑1 (MQ=255)
                                           acacTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGa                                                                                                          <  2:516174/149‑1 (MQ=255)
                                              cTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATcc                                                                                                       <  2:1206201/149‑1 (MQ=255)
                                               tGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCg                                                                                                      >  2:1233275/1‑149 (MQ=255)
                                                gggTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCgg                                                                                                     <  1:2172144/149‑1 (MQ=255)
                                                 ggTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGc                                                                                                    <  1:2812737/149‑1 (MQ=255)
                                                   tAAATCATCAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCa                                                                                                  <  1:4010465/149‑1 (MQ=255)
                                                      ataataAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCg                                                                                               >  1:1115328/1‑149 (MQ=255)
                                                       taataaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCgg                                                                                              <  1:4454540/149‑1 (MQ=255)
                                                        aataaGGAAAAGTGATGGCAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCggg                                                                                             >  2:74732/1‑149 (MQ=255)
                                                         ataaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGc                                                                                            >  1:2489365/1‑149 (MQ=255)
                                                         ataaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGc                                                                                            <  1:4688581/149‑1 (MQ=255)
                                                         ataaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGc                                                                                            >  2:1473693/1‑149 (MQ=255)
                                                         ataaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGc                                                                                            <  2:3372179/149‑1 (MQ=255)
                                                          taaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCa                                                                                           >  2:1029468/1‑149 (MQ=255)
                                                          taaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCa                                                                                           >  1:291199/1‑149 (MQ=255)
                                                           aaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCaa                                                                                          <  2:272017/149‑1 (MQ=255)
                                                           aaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCaa                                                                                          >  2:61036/1‑149 (MQ=255)
                                                            aGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAg                                                                                         <  1:4108932/149‑1 (MQ=255)
                                                             ggAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAATGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGa                                                                                        >  1:2379618/1‑149 (MQ=255)
                                                               aaaaGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGaaa                                                                                      >  1:3122578/1‑149 (MQ=255)
                                                                aaaGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACATTTTCTCGATCGTCAACACGGTCTTTATATTGATGTTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAAt                                                                                     >  1:1220542/1‑149 (MQ=255)
                                                                 aaGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAAtt                                                                                    <  1:3815812/149‑1 (MQ=255)
                                                                    tgatgaCAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCg                                                                                 >  1:3918719/1‑149 (MQ=255)
                                                                    tgatgaCAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCg                                                                                 >  1:668001/1‑149 (MQ=255)
                                                                     gatgaCAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGt                                                                                >  2:900115/1‑149 (MQ=255)
                                                                      atgaCAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTc                                                                               >  1:565182/1‑149 (MQ=255)
                                                                       tgaCAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCt                                                                              <  2:4364278/149‑1 (MQ=255)
                                                                         aCAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTAc                                                                            >  1:1686441/1‑149 (MQ=255)
                                                                         aCAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTAc                                                                            >  2:3124253/1‑149 (MQ=255)
                                                                          cAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTAct                                                                           >  2:151978/1‑149 (MQ=255)
                                                                           agagCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTActg                                                                          >  1:4353270/1‑149 (MQ=255)
                                                                             agCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTActgct                                                                        >  1:3963562/1‑149 (MQ=255)
                                                                             agCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTActgct                                                                        <  1:1889425/149‑1 (MQ=255)
                                                                                cGCATGTAGCAGTATTAAGCCAGGTCCAAAAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACTGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGAt                                                                     <  1:2135046/149‑1 (MQ=255)
                                                                                 gCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATg                                                                    <  2:3021609/149‑1 (MQ=255)
                                                                                     gTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCaa                                                                <  2:2373304/149‑1 (MQ=255)
                                                                                      tAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAAc                                                               <  2:3660345/149‑1 (MQ=255)
                                                                                       aGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGCTCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACg                                                              >  1:1607415/1‑149 (MQ=255)
                                                                                         cAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGaa                                                            <  2:1356500/149‑1 (MQ=255)
                                                                                         cAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGaa                                                            <  2:2649084/149‑1 (MQ=255)
                                                                                           gTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGc                                                          <  1:2023754/149‑1 (MQ=255)
                                                                                              ttAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGa                                                       <  1:1802669/149‑1 (MQ=255)
                                                                                                aaGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATg                                                     <  2:4180630/149‑1 (MQ=255)
                                                                                                 aGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGt                                                    >  2:3810134/1‑149 (MQ=255)
                                                                                                 aGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGt                                                    >  2:812479/1‑149 (MQ=255)
                                                                                                        tCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGACATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTCGATAAc                                             >  1:4209662/1‑149 (MQ=255)
                                                                                                            aCAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAg                                         >  2:1739296/1‑149 (MQ=255)
                                                                                                            aCAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAg                                         >  2:3261407/1‑149 (MQ=255)
                                                                                                            aCAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAg                                         >  1:4513513/1‑149 (MQ=255)
                                                                                                             cAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGt                                        >  1:350940/1‑149 (MQ=255)
                                                                                                               gTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCa                                      <  2:4307928/149‑1 (MQ=255)
                                                                                                               gTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCa                                      <  2:27973/149‑1 (MQ=255)
                                                                                                               gTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCa                                      <  2:4072426/149‑1 (MQ=255)
                                                                                                                tttCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCAt                                     <  2:3235716/149‑1 (MQ=255)
                                                                                                                       aTCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGcc                              >  2:4615077/1‑149 (MQ=255)
                                                                                                                        tCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCt                             <  1:2040124/149‑1 (MQ=255)
                                                                                                                         cGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTg                            <  2:4746640/149‑1 (MQ=255)
                                                                                                                           tCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGc                          >  2:2008703/1‑149 (MQ=255)
                                                                                                                            cAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCg                         >  1:177365/1‑149 (MQ=255)
                                                                                                                            cAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCg                         >  1:272371/1‑149 (MQ=255)
                                                                                                                               cacGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGc                      >  2:4502937/1‑149 (MQ=255)
                                                                                                                                acGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGcc                     >  1:4420255/1‑149 (MQ=255)
                                                                                                                                   gTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCttt                  <  2:292852/149‑1 (MQ=255)
                                                                                                                                      tttATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGtc               >  1:3135246/1‑149 (MQ=255)
                                                                                                                                       ttATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGtct              <  2:108908/149‑1 (MQ=255)
                                                                                                                                          tatTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTcgc           >  1:3585304/1‑149 (MQ=255)
                                                                                                                                          tatTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAATCGGTTGGCGATCTGTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTcgc           <  2:2470339/149‑1 (MQ=255)
                                                                                                                                           atTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTcgcg          >  2:3925000/1‑149 (MQ=255)
                                                                                                                                            ttGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGt         <  2:190214/149‑1 (MQ=255)
                                                                                                                                            ttGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGt         >  1:2444930/1‑149 (MQ=255)
                                                                                                                                              gATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCg       >  1:413903/1‑149 (MQ=255)
                                                                                                                                                 ggTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTg    <  1:459496/149‑1 (MQ=255)
                                                                                                                                                  gtcgtcCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTgg   <  2:3736788/149‑1 (MQ=255)
                                                                                                                                                  gtcgtcCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTgg   >  1:1088759/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGCGAAAAAAGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGG  >  CP009273/1441706‑1442001

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: