Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,590,906 C→A 100% A5A (GCG→GCT)  livK ← leucine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,590,9060CA98.7% 993.6 / 2.9 307A5A (GCG→GCT) livKleucine transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (134/169);  total (136/171)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 307 total aligned reads displayed.

AATTGCCTGACGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAAT  >  CP009273/3590762‑3591050
                                                                                                                                                |                                                                                                                                                
aaTTGCCTGACGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCa                                                                                                                                               <  1:343520/148‑1 (MQ=255)
  ttGCCTGACGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCAtt                                                                                                                                             >  2:369018/1‑148 (MQ=255)
  ttGCCTGACGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCAtt                                                                                                                                             <  1:3500745/148‑1 (MQ=255)
   tGCCTGACGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTc                                                                                                                                            >  2:2752051/1‑148 (MQ=255)
     ccTGACGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCg                                                                                                                                          <  1:1047707/148‑1 (MQ=255)
        gACGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGttt                                                                                                                                       <  1:1329236/148‑1 (MQ=255)
        gACGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGttt                                                                                                                                       <  1:3815306/148‑1 (MQ=255)
        gACGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGttt                                                                                                                                       <  2:3107492/148‑1 (MQ=255)
        gACGCGCGCCGTTAAAGTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGttt                                                                                                                                       <  2:2317010/148‑1 (MQ=255)
         aCGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTc                                                                                                                                      <  2:1520263/148‑1 (MQ=255)
           gcgcgcCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCAt                                                                                                                                    <  1:3301065/148‑1 (MQ=255)
                ccGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCa                                                                                                                               >  2:1311328/1‑148 (MQ=255)
                  gTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCaaa                                                                                                                             <  2:1624982/148‑1 (MQ=255)
                   ttAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCaaaa                                                                                                                            >  1:1079551/1‑148 (MQ=255)
                   ttAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCaaaa                                                                                                                            >  1:698222/1‑148 (MQ=255)
                   ttAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCaaaa                                                                                                                            <  2:2858556/148‑1 (MQ=255)
                   ttAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCaaaa                                                                                                                            >  1:18805/1‑148 (MQ=255)
                    tAAATTCCATATCGCCCCACTGGGCCATCGGGCCGTACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGGGTGTGAAAATGCCAGTGCAATCATCCCTGAGATGATAGTATTAGCATACCGTTTCATAGTCaaaaa                                                                                                                           >  2:3278035/1‑148 (MQ=255)
                      aaTTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATc                                                                                                                         >  1:1724764/1‑148 (MQ=255)
                       aTTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATcc                                                                                                                        >  2:2345078/1‑148 (MQ=255)
                        ttCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATccc                                                                                                                       >  1:919461/1‑148 (MQ=255)
                        ttCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATccc                                                                                                                       <  2:3027482/148‑1 (MQ=255)
                         tCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATcccc                                                                                                                      >  1:459796/1‑148 (MQ=255)
                         tCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATcccc                                                                                                                      <  2:2276122/148‑1 (MQ=255)
                           cATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCAt                                                                                                                    <  1:51720/148‑1 (MQ=255)
                                cGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTg                                                                                                               <  2:1791429/148‑1 (MQ=255)
                                cGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTg                                                                                                               <  1:753752/148‑1 (MQ=255)
                                  ccccACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAt                                                                                                             >  2:286591/1‑148 (MQ=255)
                                    ccACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgt                                                                                                           <  1:883065/148‑1 (MQ=255)
                                      aCTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttg                                                                                                         <  1:1393912/148‑1 (MQ=255)
                                      aCTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttg                                                                                                         <  1:1402089/148‑1 (MQ=255)
                                      aCTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttg                                                                                                         <  2:62518/148‑1 (MQ=255)
                                        tGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttgtg                                                                                                       <  2:3656858/148‑1 (MQ=255)
                                         gggCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttgtgt                                                                                                      <  1:928152/148‑1 (MQ=255)
                                          ggCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttgtgtt                                                                                                     >  2:1663320/1‑148 (MQ=255)
                                           gCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttgtgttg                                                                                                    <  1:1122469/148‑1 (MQ=255)
                                              aTCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCtt                                                                                                 >  1:3306648/1‑148 (MQ=255)
                                              aTCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCtt                                                                                                 <  1:2656350/148‑1 (MQ=255)
                                              aTCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCtt                                                                                                 <  2:2915974/148‑1 (MQ=255)
                                                     cGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGttttt                                                                                          <  1:2703949/148‑1 (MQ=255)
                                                     cGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGttttt                                                                                          <  1:2469720/148‑1 (MQ=255)
                                                        aCATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTAtg                                                                                       >  2:2469210/1‑148 (MQ=255)
                                                        aCATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTAtg                                                                                       >  1:2059851/1‑148 (MQ=255)
                                                           tCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATgtgt                                                                                    >  2:1755078/1‑148 (MQ=255)
                                                           tCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATgtgt                                                                                    >  1:126621/1‑148 (MQ=255)
                                                            cgcgCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGtt                                                                                   >  2:3374249/1‑148 (MQ=255)
                                                             gcgcCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTa                                                                                  <  1:797254/148‑1 (MQ=255)
                                                              cgcCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTaa                                                                                 <  1:2638630/148‑1 (MQ=255)
                                                                       cGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGa                                                                        <  1:1867940/148‑1 (MQ=255)
                                                                       cGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGa                                                                        <  2:1165719/148‑1 (MQ=255)
                                                                       cGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGa                                                                        <  1:3210782/148‑1 (MQ=255)
                                                                       cGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGa                                                                        >  2:3186832/1‑148 (MQ=255)
                                                                       cGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTACGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGa                                                                        <  1:351752/148‑1 (MQ=255)
                                                                         gCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACt                                                                      >  1:1768594/1‑148 (MQ=255)
                                                                          cGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTg                                                                     <  2:1079551/148‑1 (MQ=255)
                                                                          cGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTg                                                                     <  1:1506141/148‑1 (MQ=255)
                                                                           gACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGt                                                                    <  2:1244091/148‑1 (MQ=255)
                                                                             cTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTc                                                                  <  2:708516/148‑1 (MQ=255)
                                                                             cTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTc                                                                  >  2:2297820/1‑148 (MQ=255)
                                                                              tttAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCt                                                                 <  2:1282408/148‑1 (MQ=255)
                                                                              tttAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCt                                                                 >  2:685672/1‑148 (MQ=255)
                                                                                tAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCttt                                                               >  1:3014446/1‑148 (MQ=255)
                                                                                tAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCttt                                                               <  2:3700083/148‑1 (MQ=255)
                                                                                  atatCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCttttt                                                             >  2:2517678/1‑148 (MQ=255)
                                                                                   tatCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCtttttt                                                            <  1:1895159/148‑1 (MQ=255)
                                                                                      cGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTtat                                                         >  2:1555776/1‑148 (MQ=255)
                                                                                           gCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGc                                                    >  1:2725605/1‑148 (MQ=255)
                                                                                            ccATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCa                                                   >  1:2925853/1‑148 (MQ=255)
                                                                                            ccATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCa                                                   <  1:2353178/148‑1 (MQ=255)
                                                                                            ccATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCa                                                   >  1:1318695/1‑148 (MQ=255)
                                                                                            ccATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCa                                                   >  2:1824602/1‑148 (MQ=255)
                                                                                              aTAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACt                                                 >  2:2590950/1‑148 (MQ=255)
                                                                                               tAGCGGTGTGTGAAATTGCCAGTGCAATCATCACTGAGATGATAGTTATAGCATTCCGTTTCATAGTCAAAATTCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTATCATATCAGAGTGTTCTTTTTTTATACGGCACTg                                                <  2:1915162/148‑1 (MQ=255)
                                                                                                       tgtgAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGATTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGcc                                        <  2:2200473/148‑1 (MQ=255)
                                                                                                          gAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGt                                     >  1:3587359/1‑148 (MQ=255)
                                                                                                          gAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGt                                     <  2:211201/148‑1 (MQ=255)
                                                                                                          gAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGt                                     <  2:210595/148‑1 (MQ=255)
                                                                                                            aaTTGCCAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCt                                   <  1:3121285/148‑1 (MQ=255)
                                                                                                                  cAGTGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTaa                             <  2:506108/148‑1 (MQ=255)
                                                                                                                   aGTGCAATCATCCCTGCGATGATTGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAg                            <  2:3145411/148‑1 (MQ=255)
                                                                                                                     tGCAATCATCCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAggg                          >  1:1484970/1‑148 (MQ=255)
                                                                                                                         atcatcCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTa                      >  1:2269900/1‑148 (MQ=255)
                                                                                                                         atcatcCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTa                      >  2:3082116/1‑148 (MQ=255)
                                                                                                                          tcatcCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAg                     >  2:1882847/1‑148 (MQ=255)
                                                                                                                          tcatcCCTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAg                     >  1:2168388/1‑148 (MQ=255)
                                                                                                                              cccTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGTTGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCa                 >  2:865447/1‑148 (MQ=255)
                                                                                                                              cccTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCa                 >  2:2780759/1‑148 (MQ=255)
                                                                                                                               ccTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAg                <  1:2039095/148‑1 (MQ=255)
                                                                                                                                cTGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGt               <  2:1707798/148‑1 (MQ=255)
                                                                                                                                 tGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTa              >  2:3480833/1‑148 (MQ=255)
                                                                                                                                 tGCGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTa              <  2:2679116/148‑1 (MQ=255)
                                                                                                                                   cGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTAtt            >  1:1619184/1‑148 (MQ=255)
                                                                                                                                   cGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTAtt            >  1:2386654/1‑148 (MQ=255)
                                                                                                                                   cGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTAtt            <  1:1473359/148‑1 (MQ=255)
                                                                                                                                   cGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTAtt            <  1:3375195/148‑1 (MQ=255)
                                                                                                                                   cGATGATAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTAtt            >  1:1298469/1‑148 (MQ=255)
                                                                                                                                    gatgatAGTTTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTAttt           <  1:3753108/148‑1 (MQ=255)
                                                                                                                                         tAGATTTAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGa      <  2:1968660/148‑1 (MQ=255)
                                                                                                                                            ttttAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGataa   >  2:1433260/1‑148 (MQ=255)
                                                                                                                                            ttttAGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTATGATTAAGGGGTTAGCGCAGTATTTTGTGataa   >  2:581580/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AATTGCCTGACGCGCGCCGTTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAAT  >  CP009273/3590762‑3591050

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: ‑