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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,441,8540CA52.1% 74.6 / inf 295R27S (CGT→AGT) feaBphenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (76/65);  new base A (97/57);  total (173/122)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

Only 100 of 299 total aligned reads displayed.

GAAAAAAGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCC  >  CP009273/1441709‑1442003
                                                                                                                                                 |                                                                                                                                                     
 aaaaaaGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCTGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGAGCGTCAACACGGTCTTTATATTGATGGTCGTCCt                                                                                                                                                  >  2:243602/1‑150 (MQ=255)
 aaaaaaGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCt                                                                                                                                                  <  1:1917862/150‑1 (MQ=255)
     aaGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGcc                                                                                                                                              >  1:2950089/1‑150 (MQ=255)
     aaGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATAGATGGTAGTCCTGGcc                                                                                                                                              >  2:164356/1‑150 (MQ=255)
     aaGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTGATATTGATGGTCGTCCTGGcc                                                                                                                                              >  1:2950094/1‑150 (MQ=255)
      aGTGACTTTTCTTGTCGCTGCTTACACTGAAATCACACTGGGTAAATAATAAGGAATAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGccc                                                                                                                                             <  1:4175576/150‑1 (MQ=255)
       gTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGcccc                                                                                                                                            <  1:1772229/150‑1 (MQ=255)
       gTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGcccc                                                                                                                                            <  2:395084/150‑1 (MQ=255)
        tGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCg                                                                                                                                           <  1:1147089/150‑1 (MQ=255)
          aCTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGca                                                                                                                                         >  2:2716489/1‑150 (MQ=255)
            ttttCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGcaca                                                                                                                                       >  2:4073720/1‑150 (MQ=255)
             tttCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACaa                                                                                                                                      <  1:469689/150‑1 (MQ=255)
             tttCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACaa                                                                                                                                      <  2:2819041/150‑1 (MQ=255)
              ttCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACaaa                                                                                                                                     >  1:3633277/1‑150 (MQ=255)
              ttCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACaaa                                                                                                                                     >  1:3252999/1‑150 (MQ=255)
                 ttGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTg                                                                                                                                  <  1:136013/150‑1 (MQ=255)
                 ttGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTg                                                                                                                                  >  1:80180/1‑150 (MQ=255)
                  tGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGa                                                                                                                                 <  2:3973104/150‑1 (MQ=255)
                      gCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGaaaaa                                                                                                                             >  1:1531354/1‑150 (MQ=255)
                      gCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGaaaaa                                                                                                                             >  1:1530906/1‑150 (MQ=255)
                      gCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGaaaaa                                                                                                                             <  2:55717/150‑1 (MQ=255)
                        tGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACg                                                                                                                           <  1:1344559/150‑1 (MQ=255)
                         gagTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACgg                                                                                                                          <  2:2652735/148‑1 (MQ=255)
                         gCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACgg                                                                                                                          >  1:4107719/1‑150 (MQ=255)
                            tACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTcc                                                                                                                                                   <  1:3362604/122‑1 (MQ=255)
                              cacTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGc                                                                                                                     >  2:526405/1‑150 (MQ=255)
                              cacTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGc                                                                                                                     >  2:2522945/1‑150 (MQ=255)
                              cacTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGc                                                                                                                     >  2:2955718/1‑150 (MQ=255)
                              cacTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGGTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGc                                                                                                                     >  2:4172734/1‑150 (MQ=255)
                               acTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCg                                                                                                                    <  1:327299/150‑1 (MQ=255)
                                  gAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGAATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATc                                                                                                                 <  2:2575319/150‑1 (MQ=255)
                                   aaaTCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCt                                                                                                                >  2:2570127/1‑150 (MQ=255)
                                    aaTCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCtt                                                                                                               <  1:2222783/150‑1 (MQ=255)
                                       cacaCTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGa                                                                                                            <  2:658273/150‑1 (MQ=255)
                                        acacTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGAt                                                                                                           >  2:1997541/1‑150 (MQ=255)
                                         cacTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATc                                                                                                          <  2:1934563/150‑1 (MQ=255)
                                          acTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATcc                                                                                                         <  2:1601114/150‑1 (MQ=255)
                                          acTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATcc                                                                                                         <  1:1718010/150‑1 (MQ=255)
                                           cTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCg                                                                                                        <  1:1993870/150‑1 (MQ=255)
                                           cTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCg                                                                                                        >  1:2584180/1‑150 (MQ=255)
                                            tGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCgg                                                                                                       >  2:3598526/1‑150 (MQ=255)
                                              ggaaaaTAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCATGGCCCCGCACAAAGTGAAAAACGGGAGGCGATCTTTGATCCGGca                                                                                                     >  2:1890952/4‑149 (MQ=255)
                                                  aataatAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCg                                                                                                 >  1:593471/1‑150 (MQ=255)
                                                    taataaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCggg                                                                                               <  1:681330/150‑1 (MQ=255)
                                                     aataaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGc                                                                                              >  2:3886463/1‑150 (MQ=255)
                                                     aataaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGc                                                                                              >  2:50266/1‑150 (MQ=255)
                                                        aaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCacg                                                                                           >  2:2933216/1‑148 (MQ=255)
                                                        aaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAg                                                                                           >  1:2166352/1‑150 (MQ=255)
                                                          ggAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGaa                                                                                         >  2:4333917/1‑150 (MQ=255)
                                                              aaGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTg                                                                                     >  2:1072869/1‑150 (MQ=255)
                                                                 tgatgaCAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGt                                                                                  >  2:3361536/1‑150 (MQ=255)
                                                                   atgaCAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCt                                                                                >  2:841835/1‑150 (MQ=255)
                                                                       cAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTActg                                                                            >  2:1071271/1‑150 (MQ=255)
                                                                        agagCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTActgc                                                                           >  1:102320/1‑150 (MQ=255)
                                                                          agCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTActgctg                                                                         >  1:2428840/1‑150 (MQ=255)
                                                                           gCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGa                                                                        <  2:1795024/150‑1 (MQ=255)
                                                                            ccGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGAt                                                                       <  2:2776878/150‑1 (MQ=255)
                                                                                  gTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAAc                                                                 <  1:1028483/150‑1 (MQ=255)
                                                                                  gTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAAc                                                                 >  1:779472/1‑150 (MQ=255)
                                                                                     gCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGaa                                                              <  2:2825711/150‑1 (MQ=255)
                                                                                          aTTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGa                                                         <  2:377590/150‑1 (MQ=255)
                                                                                             aaGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGt                                                      <  1:2727932/150‑1 (MQ=255)
                                                                                             aaGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGt                                                      >  2:3120280/1‑150 (MQ=255)
                                                                                               gCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAg                                                    >  2:4374539/1‑150 (MQ=255)
                                                                                               gCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAg                                                    >  2:3386514/1‑150 (MQ=255)
                                                                                                ccAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGa                                                   >  2:1485428/1‑150 (MQ=255)
                                                                                                 cAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGAt                                                  >  2:4447993/1‑150 (MQ=255)
                                                                                                   ggTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATaa                                                >  2:3023858/1‑150 (MQ=255)
                                                                                                         aCAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGt                                          >  1:3867198/1‑150 (MQ=255)
                                                                                                           aGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCa                                        <  2:3711192/150‑1 (MQ=255)
                                                                                                           aGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCa                                        >  1:1470244/1‑150 (MQ=255)
                                                                                                            gTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCAt                                       <  2:713468/150‑1 (MQ=255)
                                                                                                            gTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCAt                                       <  1:2955718/150‑1 (MQ=255)
                                                                                                            gTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCAt                                       >  2:3064660/1‑150 (MQ=255)
                                                                                                               tctcGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTc                                    <  2:3271882/150‑1 (MQ=255)
                                                                                                                  cGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGc                                 >  1:2708316/1‑150 (MQ=255)
                                                                                                                  cGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCAGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGATGATGCCAACGAAGCGGAGGTAGATAACGCAGTCATGTCTGc                                 >  1:3194701/1‑150 (MQ=255)
                                                                                                                   gATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGcc                                >  2:1343943/1‑150 (MQ=255)
                                                                                                                    aTCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCt                               >  1:3384825/1‑150 (MQ=255)
                                                                                                                     tCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTg                              <  2:238869/150‑1 (MQ=255)
                                                                                                                      cGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTgg                             >  1:889670/1‑150 (MQ=255)
                                                                                                                      cGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTgg                             >  1:889817/1‑150 (MQ=255)
                                                                                                                       gTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGc                            >  1:558929/1‑150 (MQ=255)
                                                                                                                       gTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGCTGTAGATAACGCAGTCATGTCTGCCTGGc                            >  2:3804568/1‑150 (MQ=255)
                                                                                                                         cAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCgg                          <  2:1531354/150‑1 (MQ=255)
                                                                                                                         cAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCgg                          <  2:1530906/150‑1 (MQ=255)
                                                                                                                          aaCACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCggg                         >  2:2439704/1‑150 (MQ=255)
                                                                                                                          aaCACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCggg                         <  1:2479276/150‑1 (MQ=255)
                                                                                                                           acacGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGc                        >  1:4343443/1‑150 (MQ=255)
                                                                                                                            cacGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGcc                       >  1:3507946/1‑150 (MQ=255)
                                                                                                                                 tCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGCTGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGt                  >  2:3517006/1‑150 (MQ=255)
                                                                                                                                   tttATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGACTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGtct                >  1:4436637/1‑150 (MQ=255)
                                                                                                                                       tatTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTcgcg            <  1:2522592/150‑1 (MQ=255)
                                                                                                                                          tGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTTCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCg         <  1:4521009/150‑1 (MQ=255)
                                                                                                                                          tGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCg         <  1:3028362/150‑1 (MQ=255)
                                                                                                                                           gATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGc        >  2:1585588/1‑150 (MQ=255)
                                                                                                                                            aTGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCt       <  1:3322576/150‑1 (MQ=255)
                                                                                                                                             tGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTg      >  2:2468089/1‑150 (MQ=255)
                                                                                                                                              ggTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTgg     >  1:3642327/1‑150 (MQ=255)
                                                                                                                                              ggTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTgg     >  1:3798688/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GAAAAAAGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCC  >  CP009273/1441709‑1442003

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: