Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,959,004 G→A 100% I131I (ATC→ATT lgt ← phosphatidylglycerol‑prolipoprotein diacylglyceryl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,959,0040GA95.8% 824.7 / 17.8 239I131I (ATC→ATTlgtphosphatidylglycerol‑prolipoprotein diacylglyceryl transferase
Reads supporting (aligned to +/- strand):  ref base G (8/2);  new base A (109/120);  total (117/122)
Fisher's exact test for biased strand distribution p-value = 5.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 239 total aligned reads displayed.

GTTGGTTTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCGATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCGGCCATAAACTGCG  >  CP009273/2958857‑2959149
                                                                                                                                                   |                                                                                                                                                 
gttggttTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAAt                                                                                                                                                 <  1:3534793/150‑1 (MQ=255)
gttggttTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAAt                                                                                                                                                 <  2:606456/150‑1 (MQ=255)
gttggttTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAAt                                                                                                                                                 >  2:578977/1‑150 (MQ=255)
 ttggttTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCGATa                                                                                                                                                <  1:2721862/150‑1 (MQ=255)
  tggtgTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATaa                                                                                                                                               <  2:1038706/150‑1 (MQ=255)
    gttTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATaaaa                                                                                                                                             >  2:3637245/1‑150 (MQ=255)
    gttTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATaaaa                                                                                                                                             >  1:2804743/1‑150 (MQ=255)
    gttTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATaaaa                                                                                                                                             <  2:681387/150‑1 (MQ=255)
       tGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCa                                                                                                                                          <  2:2995093/150‑1 (MQ=255)
         cagcagcaAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCaga                                                                                                                                        <  1:3995133/150‑1 (MQ=255)
         cagcagcaAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCaga                                                                                                                                        <  2:328889/150‑1 (MQ=255)
          agcagcaAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCagag                                                                                                                                       >  2:1420199/1‑150 (MQ=255)
          agcagcaAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCagag                                                                                                                                       >  2:3194955/1‑150 (MQ=255)
             agcaAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGAcc                                                                                                                                    >  2:2494019/1‑150 (MQ=255)
             agcaAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGAcc                                                                                                                                    >  2:2909334/1‑150 (MQ=255)
              gcaAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCt                                                                                                                                   <  1:3561465/150‑1 (MQ=255)
                aaaaTATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTgg                                                                                                                                 <  2:2911301/150‑1 (MQ=255)
                aaaaTATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTgg                                                                                                                                 <  2:2911166/150‑1 (MQ=255)
                aaaaTATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTgg                                                                                                                                 >  1:2910620/1‑150 (MQ=255)
                aaaaTATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTgg                                                                                                                                 >  2:735319/1‑150 (MQ=255)
                  aaTATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGgaa                                                                                                                               >  2:3122170/1‑150 (MQ=255)
                  aaTATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGgaa                                                                                                                               <  1:1589742/150‑1 (MQ=255)
                    tatCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGgaaga                                                                                                                             >  1:3340917/1‑150 (MQ=255)
                     atCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGgaagaa                                                                                                                            <  1:380225/150‑1 (MQ=255)
                     atCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGgaagaa                                                                                                                            >  2:954962/1‑150 (MQ=255)
                     atCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGgaagaa                                                                                                                            <  1:3302602/150‑1 (MQ=255)
                          ctGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAAc                                                                                                                       >  2:1744030/1‑150 (MQ=255)
                           tGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACg                                                                                                                      <  1:720853/150‑1 (MQ=255)
                              aCGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCGATAAAATCAGAGACCTGGAAGAAGGACCGttt                                                                                                                   >  1:2492701/1‑150 (MQ=255)
                              aCGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGttt                                                                                                                   <  2:2255335/150‑1 (MQ=255)
                                gggAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAg                                                                                                                 <  1:4042613/150‑1 (MQ=255)
                                gggAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAg                                                                                                                 >  1:1728519/1‑150 (MQ=255)
                                   aGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTAc                                                                                                              <  1:4230295/150‑1 (MQ=255)
                                      cAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTAcggc                                                                                                           >  1:4634410/1‑150 (MQ=255)
                                         ggAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACggcggg                                                                                                        <  1:1501666/150‑1 (MQ=255)
                                          gAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACggcgggc                                                                                                       <  1:3178394/150‑1 (MQ=255)
                                             cAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACATTTAGTACGGCGGGCGaa                                                                                                    <  1:1869237/150‑1 (MQ=255)
                                                cATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGAt                                                                                                 <  2:1775822/150‑1 (MQ=255)
                                                cATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGAt                                                                                                 <  2:1776214/150‑1 (MQ=255)
                                                 aTGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATa                                                                                                >  2:4062885/1‑150 (MQ=255)
                                                         cGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCGATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATcac                                                                                        <  1:3094755/150‑1 (MQ=255)
                                                           ggAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATcacca                                                                                      <  2:4171169/150‑1 (MQ=255)
                                                             aaGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACg                                                                                    >  2:655556/1‑150 (MQ=255)
                                                                ttCTGGTCATCGCGGCACCACTATTCACCGTTAATAAAGTTGCCCATACGCCCGGCACCAAGACCAAACGGTATGTTTGGTTCAATAACATCATAGACCTGGAATAAGGAACGTTTAGTACGGCGGGCTAAGATAATCATAAACACGAtt                                                                                 <  2:3895736/150‑2 (MQ=255)
                                                                  cGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAAc                                                                               <  1:2513401/150‑1 (MQ=255)
                                                                        aaCGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAAt                                                                         <  2:4313999/150‑1 (MQ=255)
                                                                          cgcgGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCa                                                                       >  2:1274760/1‑150 (MQ=255)
                                                                           gcgGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAg                                                                      <  2:2859722/150‑1 (MQ=255)
                                                                           gcgGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAg                                                                      <  2:375728/150‑1 (MQ=255)
                                                                           gcgGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAg                                                                      <  1:791158/150‑1 (MQ=255)
                                                                            cgGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAgg                                                                     <  1:472914/150‑1 (MQ=255)
                                                                             ggCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGgc                                                                    >  2:520955/1‑150 (MQ=255)
                                                                               ccccACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCGATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGgccg                                                                  >  1:3171817/1‑150 (MQ=255)
                                                                               ccccACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGgccg                                                                  <  2:1632935/150‑1 (MQ=255)
                                                                                 ccACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGgccgcc                                                                <  2:2047205/150‑1 (MQ=255)
                                                                                  cacaATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGccgccg                                                               >  1:189686/1‑150 (MQ=255)
                                                                                     aaTTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTgg                                                            <  1:4049650/150‑1 (MQ=255)
                                                                                     aaTTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTgg                                                            >  1:3455254/1‑150 (MQ=255)
                                                                                     aaTTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTgg                                                            >  1:932952/1‑150 (MQ=255)
                                                                                              ttAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATg                                                   <  1:1997648/150‑1 (MQ=255)
                                                                                              ttAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATg                                                   >  1:3364009/1‑150 (MQ=255)
                                                                                               tactaaAGTTTCCCAGACGCCCGGCACCAAGACCAAACTGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACTCCAATCAGGCCGCCGTGGAAAGACATgc                                                  <  2:1474337/147‑1 (MQ=255)
                                                                                                aataaAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCGATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATgcc                                                 >  2:3917137/1‑150 (MQ=255)
                                                                                                aataaAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATgcc                                                 >  2:3805629/1‑150 (MQ=255)
                                                                                                     aGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGt                                            <  1:1070621/150‑1 (MQ=255)
                                                                                                     aGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGt                                            <  1:1070494/150‑1 (MQ=255)
                                                                                                      gTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTc                                           >  2:542890/1‑150 (MQ=255)
                                                                                                       ttGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTcc                                          >  2:656604/1‑150 (MQ=255)
                                                                                                        tGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTccc                                         >  2:1834434/1‑150 (MQ=255)
                                                                                                        tGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTccc                                         >  2:1503507/1‑150 (MQ=255)
                                                                                                        tGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTccc                                         <  1:3169965/150‑1 (MQ=255)
                                                                                                          cccAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAg                                       <  1:1218093/150‑1 (MQ=255)
                                                                                                           ccAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGa                                      >  2:3176082/1‑150 (MQ=255)
                                                                                                           ccAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGa                                      <  2:486396/150‑1 (MQ=255)
                                                                                                           ccAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGa                                      >  1:21276/1‑150 (MQ=255)
                                                                                                            cAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGac                                     <  1:4559799/150‑1 (MQ=255)
                                                                                                              gACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGacac                                   <  2:1354743/150‑1 (MQ=255)
                                                                                                                 gCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGa                                <  2:4258775/150‑1 (MQ=255)
                                                                                                                    cGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACa                             <  1:3613425/150‑1 (MQ=255)
                                                                                                                     ggCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAg                            >  2:3474027/1‑150 (MQ=255)
                                                                                                                     ggCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAg                            >  1:3884728/1‑150 (MQ=255)
                                                                                                                     ggCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAg                            >  2:3876405/1‑150 (MQ=255)
                                                                                                                      gCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGa                           <  1:2512220/150‑1 (MQ=255)
                                                                                                                          cAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGTAAAGACATGCCGCCGTCCCAGACACGGAACAGATACa                       >  1:3052208/1‑150 (MQ=255)
                                                                                                                          cAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACa                       >  2:384961/1‑150 (MQ=255)
                                                                                                                          cAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACa                       <  2:1852625/150‑1 (MQ=255)
                                                                                                                          cAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACa                       >  2:2575919/1‑150 (MQ=255)
                                                                                                                           aaGACCAAACGGACTGAGTGGTGCAAGAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTCCGGCGGGCGAAGATAATCATCACCACGATAACGCCCATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACCCGGACCCGATACAg                      >  2:3544631/1‑150 (MQ=255)
                                                                                                                           aaGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAg                      >  2:4166991/1‑150 (MQ=255)
                                                                                                                            aGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGc                     <  1:4300128/150‑1 (MQ=255)
                                                                                                                            aGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGc                     >  1:1621460/1‑150 (MQ=255)
                                                                                                                                cAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGAt                 <  2:2798216/150‑1 (MQ=255)
                                                                                                                                 aaaCGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATc                <  1:915621/150‑1 (MQ=255)
                                                                                                                                  aaCGGAATGAGTGGTGCGATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCg               >  2:2005485/1‑150 (MQ=255)
                                                                                                                                    cGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCGGc             >  2:2604110/1‑150 (MQ=255)
                                                                                                                                       aaTGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCGGCCAt          <  2:1697527/150‑1 (MQ=255)
                                                                                                                                          gAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCGGCCATaaa       >  1:1779144/1‑150 (MQ=255)
                                                                                                                                           aGTGGTGCGATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCGGCCATAAAc      >  2:3656059/1‑150 (MQ=255)
                                                                                                                                              ggtgCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCGGCCATAAACTGc   >  1:2337789/1‑150 (MQ=255)
                                                                                                                                               gtgCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCGGCCATAAACTGCg  >  1:870219/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTTGGTTTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCGATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCGGCCATAAACTGCG  >  CP009273/2958857‑2959149

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: