Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 170,868 C→A 13.3% V168F (GTT→TTT)  hemL ← glutamate‑1‑semialdehyde aminotransferase (aminomutase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273170,8680CA13.3% 844.2 / 82.7 248V168F (GTT→TTT) hemLglutamate‑1‑semialdehyde aminotransferase (aminomutase)
Reads supporting (aligned to +/- strand):  ref base C (105/110);  new base A (14/19);  total (119/129)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 249 total aligned reads displayed.

GGCAGCGGCGGAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCATGGT  >  CP009273/170720‑171015
                                                                                                                                                    |                                                                                                                                                   
   aGCGGCGGAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGcc                                                                                                                                                  <  2:3318461/149‑1 (MQ=255)
     cggcggAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCg                                                                                                                                                <  1:4033728/149‑1 (MQ=255)
      ggcggAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGTCAATCTCTTGCGGGTATTGCTCAAATGCGGCGAGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCAGGAACGCCCgg                                                                                                                                               >  1:4005760/1‑149 (MQ=255)
        cggAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCg                                                                                                                                             <  1:3241068/149‑1 (MQ=255)
         ggAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGa                                                                                                                                            >  1:476762/1‑149 (MQ=255)
          gAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAg                                                                                                                                           <  2:2422757/149‑1 (MQ=255)
           aaCACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGt                                                                                                                                          >  1:4364563/1‑149 (MQ=255)
            acacAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGtt                                                                                                                                         <  2:2853682/149‑1 (MQ=255)
             cacaGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGttt                                                                                                                                        >  2:2789850/1‑149 (MQ=255)
                aGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTtggc                                                                                                                                     <  2:4269757/149‑1 (MQ=255)
                     atatTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGcc                                                                                                                                <  2:4534156/149‑1 (MQ=255)
                     atagTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGcc                                                                                                                                <  1:1802805/149‑1 (MQ=255)
                      tatTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCt                                                                                                                               <  2:2384943/149‑1 (MQ=255)
                       atTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTa                                                                                                                              >  2:1459085/1‑149 (MQ=255)
                        ttGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTaa                                                                                                                             <  2:4337001/149‑1 (MQ=255)
                        ttGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTaa                                                                                                                             >  1:1750067/1‑149 (MQ=255)
                             tgccACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGa                                                                                                                        <  2:888887/149‑1 (MQ=255)
                              gccACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAg                                                                                                                       >  1:1109589/1‑149 (MQ=255)
                                 aCCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTCTTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTgc                                                                                                                    >  2:129152/1‑149 (MQ=255)
                                   cGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTgcgc                                                                                                                  <  2:2551240/149‑1 (MQ=255)
                                    ggCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGcc                                                                                                                 >  2:1674041/1‑149 (MQ=255)
                                     gCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCa                                                                                                                >  2:1241430/1‑149 (MQ=255)
                                     gCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCa                                                                                                                >  1:2112448/1‑149 (MQ=255)
                                         gacgaTAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAAc                                                                                                            >  2:41445/1‑149 (MQ=255)
                                         gacgaTAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAAc                                                                                                            <  2:93830/149‑1 (MQ=255)
                                            gaTAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCgg                                                                                                         >  1:3109/1‑149 (MQ=255)
                                            gaTAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCgg                                                                                                         >  2:1871722/1‑149 (MQ=255)
                                               aataCAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCtt                                                                                                      <  2:4012292/149‑1 (MQ=255)
                                                ataCAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCttt                                                                                                     >  2:2357423/1‑149 (MQ=255)
                                                  aCAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCa                                                                                                   <  1:2967021/149‑1 (MQ=255)
                                                     ggCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACca                                                                                                <  1:2506951/149‑1 (MQ=255)
                                                      gCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACcag                                                                                               >  2:3121257/1‑149 (MQ=255)
                                                      gCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACcag                                                                                               >  2:1496843/1‑149 (MQ=255)
                                                               tGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTc                                                                                      <  1:2049738/149‑1 (MQ=255)
                                                                 cGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAg                                                                                    >  1:3728220/1‑149 (MQ=255)
                                                                  gggTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGc                                                                                   <  1:1816937/149‑1 (MQ=255)
                                                                   ggTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCg                                                                                  <  1:3890237/149‑1 (MQ=255)
                                                                      aTTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGa                                                                               >  2:1276559/1‑149 (MQ=255)
                                                                        tGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGAcc                                                                             >  1:2842855/1‑149 (MQ=255)
                                                                        tGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGAcc                                                                             >  1:837681/1‑149 (MQ=255)
                                                                          cTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGTCAGTCAGCGTGACCAt                                                                           <  2:1311113/149‑1 (MQ=255)
                                                                          cTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCAt                                                                           <  1:2175806/149‑1 (MQ=255)
                                                                          cTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCAt                                                                           <  1:3293614/149‑1 (MQ=255)
                                                                              aaTGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCGGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTa                                                                       <  1:648791/149‑1 (MQ=255)
                                                                              aaTGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTa                                                                       >  2:657740/1‑149 (MQ=255)
                                                                              aaTGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATAGTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTa                                                                       <  1:2988704/149‑1 (MQ=255)
                                                                                   ggcgCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACAcc                                                                  >  2:1390720/1‑149 (MQ=255)
                                                                                     cgcgTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCt                                                                <  1:150817/149‑1 (MQ=255)
                                                                                       cgTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTc                                                              <  2:2165316/149‑1 (MQ=255)
                                                                                       cgTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTc                                                              <  1:3608104/149‑1 (MQ=255)
                                                                                          aCAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCaaa                                                           <  2:1750067/149‑1 (MQ=255)
                                                                                            aGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAAtt                                                         <  2:3602450/149‑1 (MQ=255)
                                                                                            aGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAAtt                                                         >  1:4272535/1‑149 (MQ=255)
                                                                                            aGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAAtt                                                         <  2:1827984/149‑1 (MQ=255)
                                                                                                gCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTtaat                                                     <  2:178705/149‑1 (MQ=255)
                                                                                                 ccAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTtaata                                                    <  1:3806421/149‑1 (MQ=255)
                                                                                                  cAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTtaataa                                                   >  2:408042/1‑149 (MQ=255)
                                                                                                   aGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTaataat                                                  >  1:4628674/1‑149 (MQ=255)
                                                                                                    gATCATTATAAGTACATGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAAtt                                                 <  2:4005760/149‑1 (MQ=255)
                                                                                                     aTCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAAttt                                                >  2:300120/1‑149 (MQ=255)
                                                                                                     aTCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAAttt                                                >  1:3366453/1‑149 (MQ=255)
                                                                                                     aTCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAAttt                                                >  2:3690450/1‑149 (MQ=255)
                                                                                                        aTTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGt                                             <  1:1579354/149‑1 (MQ=255)
                                                                                                         ttATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTc                                            <  2:3503542/149‑1 (MQ=255)
                                                                                                         ttATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTc                                            >  2:1432749/1‑149 (MQ=255)
                                                                                                         ttATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTc                                            >  1:544179/1‑149 (MQ=255)
                                                                                                         ttATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTc                                            <  2:2502288/149‑1 (MQ=255)
                                                                                                            taAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTcgcg                                         >  1:1693561/1‑149 (MQ=255)
                                                                                                            taAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTcgcg                                         >  2:3547922/1‑149 (MQ=255)
                                                                                                            taAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTcgcg                                         >  1:1008236/1‑149 (MQ=255)
                                                                                                                 aCAGGTTAAGGTATATTTGGCGAAATCTGCTGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCgg                                    <  1:2326863/149‑1 (MQ=255)
                                                                                                                 aCAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCgg                                    >  1:1199020/1‑149 (MQ=255)
                                                                                                                  cAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGt                                   <  1:4146543/149‑1 (MQ=255)
                                                                                                                  cAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGt                                   <  1:4393099/149‑1 (MQ=255)
                                                                                                                   atgTTAAGGTATTTTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCTGTTCGGGCTTTGTCCAGCAGGCAGTCAGCGTGTCCATGGTATCATCGTTCATATTTTATAATTTTTTCTCGACTGGtt                                  <  2:1313742/147‑2 (MQ=255)
                                                                                                                    ggTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCAGGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTaa                                 >  2:2016926/1‑149 (MQ=255)
                                                                                                                     gTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTaaa                                >  2:121856/1‑149 (MQ=255)
                                                                                                                        aaGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAAcc                             >  1:731228/1‑149 (MQ=255)
                                                                                                                          ggTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCAc                           <  1:2682916/149‑1 (MQ=255)
                                                                                                                            tataTTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACgg                         <  1:1269354/149‑1 (MQ=255)
                                                                                                                             atatTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACggg                        >  1:470412/1‑149 (MQ=255)
                                                                                                                             atatTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTGAATAATTTTGTCGCGACCGGTACAACCACggg                        >  1:1347209/1‑149 (MQ=255)
                                                                                                                               atTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGcc                      >  2:153139/1‑149 (MQ=255)
                                                                                                                               atTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGcc                      <  1:3915684/149‑1 (MQ=255)
                                                                                                                                tttGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCa                     <  2:1004509/149‑1 (MQ=255)
                                                                                                                                tttGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCa                     >  2:1921800/1‑149 (MQ=255)
                                                                                                                                  tGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCACATTTAATCATTTTGTCGCGACCGGTAAAACCACGGGCCAgg                   >  2:4400147/1‑149 (MQ=255)
                                                                                                                                  tGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGACAGCGTGACCATGGTCACACCCTTCAAAGTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAgg                   >  1:2633269/1‑149 (MQ=255)
                                                                                                                                   ggCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCAGTAAAACCACGGGCCAGGc                  <  1:2800187/149‑1 (MQ=255)
                                                                                                                                      gAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGa               <  1:4563547/149‑1 (MQ=255)
                                                                                                                                        aaTCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATg             >  1:3561710/1‑149 (MQ=255)
                                                                                                                                          tCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGgc           >  2:3179409/1‑149 (MQ=255)
                                                                                                                                            tGCCGGACCGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGgcgc         <  2:2309383/149‑1 (MQ=255)
                                                                                                                                            tGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCTCGACCGGTAAAACCACGGGCCAGGCGGATGgcgc         >  2:3220108/1‑149 (MQ=255)
                                                                                                                                            tGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGgcgc         <  1:2334840/149‑1 (MQ=255)
                                                                                                                                            tGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGgcgc         <  1:2513715/149‑1 (MQ=255)
                                                                                                                                            tGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGgcgc         >  2:1568370/1‑149 (MQ=255)
                                                                                                                                               cGGAACGCCCGGCGAGTTTGGCTGGCCTCACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCa      <  2:474130/149‑1 (MQ=255)
                                                                                                                                               cGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCGTGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCa      <  2:3840950/149‑1 (MQ=255)
                                                                                                                                                   aCGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCATGGt  >  1:2208919/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCAGCGGCGGAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCATGGT  >  CP009273/170720‑171015

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: