Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 777,175 A→T 11.7% intergenic (+67/‑366) lysQ → / → nadA tRNA‑Lys/quinolinate synthase, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273777,1750AT11.7% 345.0 / 4.0 165intergenic (+67/‑366)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base A (94/42);  new base T (3/22);  total (97/68)
Fisher's exact test for biased strand distribution p-value = 9.39e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

CAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTG  >  CP009273/777029‑777324
                                                                                                                                                   |                                                                                                                                                     
cAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGt                                                                                                                                                    <  1:2031892/150‑1 (MQ=35)
cAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTGTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGt                                                                                                                                                    >  2:2969067/1‑150 (MQ=25)
 aaGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGtt                                                                                                                                                   <  2:661716/150‑1 (MQ=34)
  aGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                                                                                  <  1:3923450/150‑1 (MQ=33)
   ggggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttg                                                                                                                                                 >  2:1554189/1‑150 (MQ=32)
      gTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagt                                                                                                                                              >  2:3459874/1‑150 (MQ=32)
         gTTAGCTCCGTTGGTAGAGCAGTTGACTATTACTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                                                                           <  2:2627393/150‑1 (MQ=12)
            aGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcg                                                                                                                                        >  1:847050/1‑150 (MQ=32)
            aGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAACTGAAGCAGGTTTGAGTTGAgcg                                                                                                                                        >  1:2678469/1‑150 (MQ=17)
             gCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcgc                                                                                                                                       >  1:3146542/1‑150 (MQ=32)
                 aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATATGCGATACTCAAAGATACGATGCTGCAGCAGGTTTGAGTTGAGCGCAGCg                                                                                                                                   >  2:1992747/1‑150 (MQ=255)
                    tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                                                                                                                >  1:2163054/1‑150 (MQ=32)
                    tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACAAAAAGATACGAACCTGAAGCAGGTTAGAGGTGAGCTCAGCGaaa                                                                                                                                >  1:3762815/1‑150 (MQ=11)
                     ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTCTCTGCGATACACAAAGATACGAACCTACAGCAGTTTTGAGATGAGCGCAGCGCaac                                                                                                                               >  2:1878478/1‑150 (MQ=11)
                         gagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGATGAGCGCAGCGAACCAACg                                                                                                                           >  2:2105184/1‑150 (MQ=17)
                         gagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCAGCGATACTGAAATTTGAGAAACTGCAGCAGGTTAGAGATGAGCGCATCGaagaaccg                                                                                                                           >  2:2510195/1‑144 (MQ=11)
                            cAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGTTTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAg                                                                                                                        >  2:1052222/1‑150 (MQ=17)
                            cAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAg                                                                                                                        >  1:4351909/1‑150 (MQ=32)
                             aGTTGACTTTTAATCAATTGGTCGCAGGTTCTAATCCTGCACGACCCACCAATTTAAAAAAGCGCCATAAAGGCGATTTTTTGCTATCTGCTATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGc                                                                                                                       <  2:473846/150‑1 (MQ=11)
                                            aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                                                                        >  1:1001582/1‑150 (MQ=32)
                                               tGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGccc                                                                                                     >  2:2024413/1‑150 (MQ=32)
                                                ggTCGCTGGTTCGAATCCTGCCCGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTGATATCTGCTATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCg                                                                                                    <  2:3132687/150‑1 (MQ=11)
                                                ggTCGCTGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCg                                                                                                    <  1:3369991/150‑1 (MQ=17)
                                                ggTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGTACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCg                                                                                                    >  1:678822/1‑150 (MQ=17)
                                                ggTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCg                                                                                                    >  2:2357174/1‑150 (MQ=32)
                                                  tCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGaa                                                                                                  >  2:2852626/1‑150 (MQ=32)
                                                    gCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAgg                                                                                                >  1:2886818/1‑150 (MQ=32)
                                                     cAGGTTCTAATCCTGCACGACCCACCAAGTTAAAAACGCGCCCTAAAGGCGATTTTTTGCTTTCTGGTATACTCAAAGATTCGACCCTTCAGCAGTTTTGAGTTGAGCGCAGCGAAAAAACGGAGCCGCTCGCGGCGACGGCCCGAAggg                                                                                               <  2:3833495/150‑1 (MQ=11)
                                                     cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTGAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg                                                                                               >  1:1409318/1‑150 (MQ=17)
                                                       ggTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcg                                                                                             >  1:3504998/1‑150 (MQ=32)
                                                        gtactaaTCCTGCACGACCCACCAATGAAAAAAAGCGCCCTAATGGCTCTTTTTTGCTATCTGCGATACTCAACGATTCGAACCGGCAGCTGGTTTTAGTTGCGCGCAGCGAAACCACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcga                                                                                            <  2:3124604/145‑1 (MQ=11)
                                                         gtccaccccTGCACCACCCACCAATTAAAAAACGCCCCCTACCGCCGATTTTTTGCTATCTGCGATACTCAAATATCCTACCCGGCCGCTGTTTTGTTTTGCGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                                                                           <  1:2736790/143‑1 (MQ=2)
                                                           ctaaTCCTGCACGACCCACCAATGTAAAAACGCGCCCTACCGGCGCTTTTTTGCTATCTGCGATACTCAATGATTCTACCCTGCAGCAGGTTTGAGTTGCGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcg                                                                                         <  1:645049/148‑1 (MQ=11)
                                                           cGAATCCTGCACGACCCACCAATGTAAAAAAGCTCCCTAAAGGCGCTTTTATGCTATCTGCGATACTCAAAGAGTCGAACCTGCAGAAGCTTTGAGTTGAGCGCAGCGAAGAAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcg                                                                                         >  2:4528303/1‑150 (MQ=11)
                                                           cGAAGCCTGCACGACCCCCCAATGTAAAAAAGCGCCCTCAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGCGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcg                                                                                         <  2:2761404/150‑1 (MQ=11)
                                                           cGAACCCTGCACGACCCACAATTGTAAAAAAGCGCCCTACAGGCGTTTTTTTGGTATCTGCTATACTCAAATATTCGAACCTGCAGCTGGTTTGTGTTGCGCGCCGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcg                                                                                         <  1:4428647/150‑1 (MQ=11)
                                                            taaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCAAAAGGCGATTTTTTGCTATCTGCGATACTCAAATATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcga                                                                                        <  1:357661/149‑1 (MQ=12)
                                                            gAATCCTGCCCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACGACAAGTTGCGAACCTCCAGCTGGTTTGAGTTGACCGCAGAGAACCGCCGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcga                                                                                        >  2:4580328/1‑150 (MQ=11)
                                                            gAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCCAAGATTCGAACCTGCAGCAGGTTTGAGGTGAGCGCATCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcga                                                                                        >  2:1513486/1‑150 (MQ=12)
                                                            gAATCCTGCACGACCCACCAAGGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAGAGATTCGAAGCTGCAGCAGGTTTGAGGTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcga                                                                                        >  2:4580310/1‑150 (MQ=255)
                                                            gAATCCTGCACGACCCACAAATGTAAAAAAGCGCCCTAAAGGCGATTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcga                                                                                        <  1:357591/150‑1 (MQ=17)
                                                             aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaa                                                                                       <  1:4169368/150‑1 (MQ=32)
                                                               tCCTGCACGACCCACCAATTTAAAAATGCGCCCTAACGGCTCTTTTTTGCTATATGCTATACTCAACGATTCTACCCTGAAGCCGGTTTTAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagc                                                                                     <  1:2486711/150‑1 (MQ=11)
                                                                ccTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcg                                                                                    >  1:4579482/1‑150 (MQ=32)
                                                                 cTGCACGACCCACAAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcga                                                                                   <  2:1039998/150‑1 (MQ=17)
                                                                  tGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAg                                                                                  >  2:3577005/1‑150 (MQ=32)
                                                                    cACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTc                                                                                <  2:4062122/150‑1 (MQ=32)
                                                                        aCCCACCAATTAAAAAAAGCGCCCCAACGGCGTTTTTTTGCTATCTGCGATACTAAAATATCCCAACCTGCAGCTGTTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATcc                                                                            <  1:1738877/150‑1 (MQ=11)
                                                                        aCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCCAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATcc                                                                            >  1:511838/1‑150 (MQ=17)
                                                                         cccACCAATTTAAAAAACCGCCCTAACGGCGCTTTTTTGCTATCTGCTATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCt                                                                           <  1:1776730/150‑1 (MQ=11)
                                                                         ccaaccaATTAAAAAATCCGCCCAAACGGGGCTTTTTTGCTATCTGCCATACACAAATATTCCACCCGGCAGCTGTTTTGTGTGGACCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCt                                                                           <  1:4271909/147‑1 (MQ=11)
                                                                          cccccaATTTAAAAAACCCGCCCTAACGCCGATTTTTTGCTTTCTGCTATACTCAATTATTCCACCCTGCAGCGGTTTTGTGTTGAGCGCCGCCAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCt                                                                           <  1:4271880/147‑1 (MQ=11)
                                                                           caccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGc                                                                         >  2:4216445/1‑150 (MQ=32)
                                                                           caccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCAGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGc                                                                         >  2:2426804/1‑150 (MQ=17)
                                                                             ccaCTGTAAAAAAGCGCCCTACCGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCAc                                                                       <  2:3676168/150‑1 (MQ=255)
                                                                              caATTTAAAAAAGCGCCCTAACGGCGATTTTTTGCTATCTGCGATACTCAAATATTCGAACCTGCAGCTGTTTTGAGTTGAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACg                                                                      <  1:2282086/150‑1 (MQ=11)
                                                                               aaTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGa                                                                     >  1:1993373/1‑150 (MQ=32)
                                                                                aTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAc                                                                    >  1:1031795/1‑150 (MQ=32)
                                                                                  gTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAccc                                                                  >  2:3796283/1‑150 (MQ=32)
                                                                                    aaaaatGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac                                                                >  2:3108774/1‑150 (MQ=12)
                                                                                       aaaCCGCCCAAACCGCGCTTTTTTGCTATCTGCGATACCAAAAGATTCCAACCTGCAGCTGTTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa                                                             <  1:3164438/150‑1 (MQ=11)
                                                                                        aaGCGCCCTACCCGCGCTTTTTTGCTATCTGCCACACTCAAAGATTCGAACGTGCAGCTGGTTTTAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt                                                            <  1:1068233/150‑1 (MQ=11)
                                                                                        aaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt                                                            >  1:4246063/1‑150 (MQ=32)
                                                                                            cccccTAAAGGCGCTTTTTGGCTATCTGCGATACCAACAGTTTCGAACCTCCAGAAGTTTGGAGTTGCGCGCTGCGAATCTACTGAGCGGCTCTCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTa                                                         <  2:2549928/148‑1 (MQ=11)
                                                                                             gcCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaa                                                        >  1:4083273/1‑150 (MQ=32)
                                                                                             gcCCTAAAGGCGCTTTTTTGCTATATGCGATACTCAAATATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCCCCAATGTaa                                                        >  1:1358093/1‑150 (MQ=12)
                                                                                               ccaccccgcgcTTTTTTGCTCTCTGCGATACTCAACGATTCGAACCTGCTGGAGGTTTGAGTTGAGCGCAGCGAAACAACGTAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaa                                                      <  2:1832153/143‑1 (MQ=11)
                                                                                               ccTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaa                                                      >  2:2450804/1‑150 (MQ=32)
                                                                                                cTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaa                                                     >  1:1583083/1‑150 (MQ=32)
                                                                                                   atcgcgcTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAgc                                                  <  1:3836812/147‑1 (MQ=17)
                                                                                                       cgcTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGccc                                              >  2:1313249/1‑150 (MQ=32)
                                                                                                        ttttttttGCTTTTTTCTATACTCAAAGATACCAAATTTCAGCTGTTTTTAGTTTTTCCCAGCTAAACATATTAGCGGCTGGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCt                                             <  1:2340754/148‑1 (MQ=2)
                                                                                                        tcTTTTTTGCTTTCTGCGATACTCAAAGATTCGAACCTGCACCTGGTTTTAGTTGCTCTCAGCGAAACTACGGAGCCGCTCGCGGCTACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCt                                             <  2:2180712/149‑1 (MQ=11)
                                                                                                        gcTTTCTTGCTCCCTGCTATACTAAAAGCTTCGATCCTGCAGCTTGTTTTAGTTGCTGTCAGATAAACAACGGAGCGGCTCGCGGCGAGGGCCATAAGGGCGTTCGAAGCTTGTAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCt                                             <  2:2327500/150‑1 (MQ=2)
                                                                                                         cTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTa                                            >  1:1112265/1‑150 (MQ=32)
                                                                                                         cTTTATTGCTATCTGCGATACTAAAAGTTTCGAACCTGCAGCAGTTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTa                                            <  1:4608042/150‑1 (MQ=255)
                                                                                                          ttttttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaa                                           >  1:2416585/1‑150 (MQ=32)
                                                                                                          ttttttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaa                                           >  1:2081804/1‑150 (MQ=32)
                                                                                                           tttttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaac                                          >  1:2816201/1‑149 (MQ=17)
                                                                                                             tttGCTCTCTGCGATACTCAAATATTCTAACCTGCTGCAGTTTTTTGTTTCGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGTCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAgg                                        <  2:3472993/150‑1 (MQ=11)
                                                                                                                              tCAAAGATTCGATCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCa                       <  2:4258837/150‑1 (MQ=17)
                                                                                                                                 aaGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCGAAAGGCGCTTTTTTGCTATTCAGGc                    >  2:4367341/1‑150 (MQ=21)
                                                                                                                                  ctaTTCGAACCTGAAGCAGGTTTGAGTTGAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCa                   <  1:3608210/148‑1 (MQ=255)
                                                                                                                                  ataTTCGAACCTGCTGCAGGTTTGATTTGAGCGCAGCTAACCAACGGAGCCGCTCGCGGCGACGGCCCGACGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCa                   <  1:1270945/148‑1 (MQ=11)
                                                                                                                                  aGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCa                   >  2:4024307/1‑150 (MQ=35)
                                                                                                                                  aGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCa                   >  2:4024203/1‑150 (MQ=35)
                                                                                                                                     ttCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGTGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACAAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATcc                <  1:1802592/150‑1 (MQ=255)
                                                                                                                                     atCGAACCTGCAGCAGGATTGAGTTGAGCGCAGTGAACCAAAGGAGCCGCTCGCGGCGACGGACCTTACGGCGAGCGAAGCGAGTCATCCTGCACGACCCACAAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATcc                <  1:1802593/149‑1 (MQ=11)
                                                                                                                                      tCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCt               >  2:3181973/1‑150 (MQ=38)
                                                                                                                                       cGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTc              <  1:1427622/150‑1 (MQ=39)
                                                                                                                                        gAACCTGCAGCCGGTTTGAGTTGAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCa             <  1:2726495/150‑1 (MQ=255)
                                                                                                                                        gAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCa             >  1:916375/1‑150 (MQ=40)
                                                                                                                                         acccTGCAGCAGGTTTGAGTTGAGCGCAGCGTAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCaa            <  2:1092502/148‑1 (MQ=255)
                                                                                                                                         aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCaa            >  1:3332964/1‑150 (MQ=40)
                                                                                                                                          aCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAAt           >  1:2112897/1‑150 (MQ=255)
                                                                                                                                           ccTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAAtt          <  1:1124522/150‑1 (MQ=255)
                                                                                                                                                agcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACt     >  2:4020820/1‑150 (MQ=255)
                                                                                                                                                agcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACt     >  1:1351669/1‑150 (MQ=255)
                                                                                                                                                   agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGACGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTg  <  1:4222247/150‑1 (MQ=255)
                                                                                                                                                   agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTg  >  2:1263597/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTG  >  CP009273/777029‑777324

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: