Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 777,175 A→T 8.6% 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,1750AT 8.6% 337.6 / 5.6 139intergenic (+67/‑366)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base A (97/26);  new base T (2/13);  total (99/40)
Fisher's exact test for biased strand distribution p-value = 9.40e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

CCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTG  >  CP009273/777026‑777324
                                                                                                                                                     |                                                                                                                                                     
  gCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAgg                                                                                                                                                     <  1:332150/150‑1 (MQ=35)
  gCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAgg                                                                                                                                                     >  1:1732141/1‑150 (MQ=35)
    acgtggtggtTAGATCAGTTGGTAGAGCAGTTGACTTTTAATTAATTGGTCGCAGGTTCGTAGGGTGCTCGACCCTCCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATATGGGATACTCAAAGTTTCGATCCTGCTGGAGGtt                                                                                                                                                   <  1:1551258/146‑1 (MQ=11)
    aaGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGtt                                                                                                                                                   <  1:4145633/150‑1 (MQ=34)
      tgggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttg                                                                                                                                                 <  2:1180626/149‑1 (MQ=17)
      ggggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTAGCAGGGTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCGAAAGGCGCTTTTTTGCTATCTTCGATACTCAAAGATTAGAACCTGCAGAAGGTttg                                                                                                                                                 >  2:1487782/1‑150 (MQ=11)
        ggTCTTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag                                                                                                                                               <  1:4421619/150‑1 (MQ=17)
         gTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagt                                                                                                                                              >  2:3096650/1‑150 (MQ=32)
         gTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagt                                                                                                                                              <  1:915116/150‑1 (MQ=32)
            gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                                                                           >  1:1216369/1‑150 (MQ=32)
             ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                                                                          >  2:5115707/1‑150 (MQ=32)
             ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                                                                          >  2:5115584/1‑150 (MQ=32)
             ttAGCTCAGTTGGTAGAGCAGTAGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                                                                                                                          >  2:2538156/1‑150 (MQ=17)
              tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgc                                                                                                                                         <  2:4612510/150‑1 (MQ=32)
              tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgc                                                                                                                                         <  2:4612174/150‑1 (MQ=32)
                gCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTCAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcgc                                                                                                                                       <  1:2680879/150‑1 (MQ=17)
                    aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCg                                                                                                                                   >  2:2940151/1‑150 (MQ=32)
                    aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCAGCGATACTCAAAGATAAGAAACTGCAGCTGTTTTGAGTTGAgcgccgag                                                                                                                                   >  2:327605/1‑146 (MQ=11)
                     gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGTTTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGa                                                                                                                                  >  2:550585/1‑150 (MQ=17)
                     gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGa                                                                                                                                  >  1:4276485/1‑150 (MQ=32)
                     gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATACGATCCTGCAGCAGGTGTGAGTTGATCGCAGCGa                                                                                                                                  >  2:175794/1‑150 (MQ=255)
                      ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaa                                                                                                                                 >  2:1116481/1‑150 (MQ=32)
                       tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaca                                                                                                                                >  2:3811284/1‑148 (MQ=17)
                       tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTGTGAGGTGAGCGCAGCGaaa                                                                                                                                >  2:4882014/1‑150 (MQ=17)
                        ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACACAAAGATTCGAACCTGAAGCAGGTTTGAGGTGAGCGCAGAGAaac                                                                                                                               >  1:3291462/1‑150 (MQ=255)
                          tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaa                                                                                                                             >  2:1857080/1‑150 (MQ=32)
                           agagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaac                                                                                                                            <  1:1477503/150‑1 (MQ=32)
                           agagCAGTTGACTTTTAATCAATTGGTCGCAGGTTAGAATACTGCACGACCATCCAATGTAATAAAGCGACCTACAGGCGCTTTTTTGCTATCTGATATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaac                                                                                                                            <  1:4073576/150‑1 (MQ=11)
                                   tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCAATCTGCGATAATCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGc                                                                                                                    >  1:4840028/1‑150 (MQ=17)
                                        tttAATCATTGGGTCGCTGGTTCTAATCCGGCACGCCCCACCAATTTAAAAAAGCGCCCTAACGGCGCTTTTTTGCTATGTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGTGTTGAGCGCAGCGAAACAACGGAGCCGCTcgcg                                                                                                               <  2:2480792/150‑1 (MQ=11)
                                         ttAATCAATTGGTCGAAGTTTCTAATCCTGCACGACCCACCACTTTAAAAAAGCGCCCTAACGGCGATTTTTTGCTATATGCGATAATCAAAGATTCGAACCTGCAGCAGTTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcgg                                                                                                              <  2:4855959/150‑1 (MQ=11)
                                          tAATCAAGGGGTCGCAGGTTCGAATCCGGCACGCCCCACCAAGGTAAAAAAGCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGACCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggc                                                                                                             <  2:4149394/150‑1 (MQ=11)
                                           aaTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggcg                                                                                                            >  1:3071445/1‑150 (MQ=32)
                                               aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                                                                        >  1:983157/1‑150 (MQ=32)
                                                aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGc                                                                                                       >  2:2610086/1‑150 (MQ=32)
                                                  tGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGGTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGccc                                                                                                     >  2:4039656/1‑150 (MQ=17)
                                                   ggTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCGCGCGGCGCCGGCCCg                                                                                                    >  1:3521380/1‑150 (MQ=17)
                                                    gTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGa                                                                                                   >  2:3852563/1‑150 (MQ=32)
                                                       gCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAgg                                                                                                >  1:2986189/1‑150 (MQ=32)
                                                       gCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAgg                                                                                                >  2:2212159/1‑150 (MQ=32)
                                                        cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg                                                                                               >  2:3848656/1‑150 (MQ=32)
                                                        cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg                                                                                               >  2:1252261/1‑150 (MQ=32)
                                                          ggTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcg                                                                                             >  1:2245665/1‑150 (MQ=32)
                                                           gTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcga                                                                                            >  1:2272233/1‑150 (MQ=32)
                                                            ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                                                                           >  1:1331146/1‑150 (MQ=32)
                                                                attCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCTAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaa                                                                                       <  2:3403491/148‑1 (MQ=17)
                                                                aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcgac                                                                                       >  2:451615/1‑149 (MQ=17)
                                                                 aTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACACAGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaag                                                                                      >  1:880219/1‑150 (MQ=17)
                                                                     tGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAg                                                                                  >  1:2672009/1‑150 (MQ=32)
                                                                     tGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGGTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAg                                                                                  >  1:54818/1‑150 (MQ=17)
                                                                       cACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTc                                                                                >  2:2933271/1‑150 (MQ=32)
                                                                       cACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATAATCAAAGATTCGAACCGGCAGCAGGTTTGAGGTGAGCGCAGCGAGACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTc                                                                                >  1:1789892/1‑150 (MQ=255)
                                                                         cgccccACAAATTAAAAAACGCGCCATAACGGCTTTTTTTTGATATCTCCTATCCACAACTATACGAACCTGCAGCTGTTTTGAGTTGTGCGCAGCGAAAAAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCAt                                                                              <  1:3582938/147‑1 (MQ=2)
                                                                          gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATc                                                                             >  2:1870098/1‑150 (MQ=32)
                                                                            cccACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCt                                                                           >  1:850185/1‑150 (MQ=32)
                                                                             ccaccaATGTAAAAAAGCGCCCTACAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTg                                                                          <  2:5171167/150‑1 (MQ=17)
                                                                                ccaTTGTAAAAAAGCGCCCTACAGGCTTTTTTTTGCTATCTGCGATACCCAAAGTTTCGACCGTGCTGCTGTTTTTAGTTGAGGTCAGCGACACTTCGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCAc                                                                       <  2:2426663/150‑1 (MQ=11)
                                                                                ccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCAc                                                                       >  2:3907111/1‑150 (MQ=32)
                                                                                 caATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCCTCCTGCACg                                                                      >  2:3648058/1‑150 (MQ=17)
                                                                                   ttGTAAAAAAGCGCCCTAACGGCGATTTTTTGCTATCTGCGATACTCAAATATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAc                                                                    <  1:3874982/149‑1 (MQ=255)
                                                                                    tGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAcc                                                                   >  1:5111070/1‑150 (MQ=32)
                                                                                    tGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTCGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAcc                                                                   >  1:5110385/1‑150 (MQ=17)
                                                                                     gTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAccc                                                                  >  1:3103447/1‑150 (MQ=32)
                                                                                     gTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAACGAAGCGAGTCATCCTGCACGAccc                                                                  >  1:622686/1‑150 (MQ=17)
                                                                                       aaaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac                                                                >  2:3444928/1‑150 (MQ=32)
                                                                                        aaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc                                                               >  1:1185617/1‑150 (MQ=32)
                                                                                        aaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc                                                               >  1:1185654/1‑150 (MQ=32)
                                                                                        aaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc                                                               >  2:1921060/1‑150 (MQ=32)
                                                                                        aaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc                                                               >  2:1775732/1‑150 (MQ=32)
                                                                                           aaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCCCCAAt                                                            >  1:2820520/1‑150 (MQ=17)
                                                                                           aaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt                                                            >  2:606104/1‑150 (MQ=32)
                                                                                           aaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt                                                            >  1:4319222/1‑150 (MQ=32)
                                                                                           aaGCGCCCAAAAGGCGCTTTTTTGCTATCAGCCAAACTCAAAGATTCGAACCTGCAGCAGTTATGAGTTGAGCGCCGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt                                                            <  2:2758978/150‑1 (MQ=11)
                                                                                            aGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCCTCCTGCACGACCCACCAATg                                                           >  1:1961616/1‑150 (MQ=17)
                                                                                            aGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATg                                                           >  1:1961576/1‑150 (MQ=32)
                                                                                               gcCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaa                                                        >  1:492301/1‑150 (MQ=32)
                                                                                                 ccTAACGTCGATTTTTTGCTATGTGCTATACACAAATATCCTAACCTTCAGGTGTTTTTTGTTTAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGACGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaa                                                      <  1:4147198/150‑1 (MQ=11)
                                                                                                 ccTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaa                                                      >  1:1698248/1‑150 (MQ=32)
                                                                                                          gcTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCCATGTAAAAAAGCGCCCt                                             >  2:2908312/1‑150 (MQ=17)
                                                                                                          gcTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCt                                             >  2:4156466/1‑150 (MQ=32)
                                                                                                          gcTTTTTTGCTATCTGCGATACTAAATGTTGCGACCCTGCAGCTGGTTTGCGTTGAGCGCAGCTAAACAACGGAGCGGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCt                                             <  2:2229209/150‑1 (MQ=11)
                                                                                                          gcTTTTTTGCTATATGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCt                                             <  2:4208722/150‑1 (MQ=17)
                                                                                                           cTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTa                                            >  1:121258/1‑150 (MQ=32)
                                                                                                           cTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTa                                            >  1:3300826/1‑150 (MQ=32)
                                                                                                                 tGCTATCTGGTATACTCAAAGATTCTAACCTGCAGCAGGTTTGAGTTGAGCGCTGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                      <  1:5115707/150‑1 (MQ=11)
                                                                                                                 tGCTATCTGCTATACTAAAATTTTCCAACCTGATGATGGTTTGTGTTGCGCGAAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                      <  1:2538156/150‑1 (MQ=11)
                                                                                                                              acacAATGATTCTAACCTGCTGCCGGTTTTAGTTGCGCGCCCCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAtt                         <  1:1058402/147‑1 (MQ=11)
                                                                                                                                 cAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTTCTATTCAg                      >  1:4578157/1‑150 (MQ=17)
                                                                                                                                 cAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAg                      >  1:3315393/1‑150 (MQ=32)
                                                                                                                                  aaaTATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCTAAACAACGGAGCCGCACGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATTTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAgg                     <  2:3187613/150‑1 (MQ=255)
                                                                                                                                   acgagtCGCACCTGCCCCAGGTTTGCGTTTAGCGCCGCGAATCAACGGAGCCGCTCGCGGCGAGGGCCAGACGGGCGAGCGCAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGc                    <  2:3678500/145‑1 (MQ=11)
                                                                                                                                    aGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCa                   >  2:56879/1‑150 (MQ=35)
                                                                                                                                      tttCGACCCAGCAGCATGTTTGCGTTGAGGGCAGGTAACCACCGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATc                 <  2:2015512/149‑1 (MQ=11)
                                                                                                                                         cGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTc              >  2:3172164/1‑150 (MQ=39)
                                                                                                                                                gcagctGGTTTTTGTTGCGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCa       <  1:4007494/150‑1 (MQ=14)
                                                                                                                                                 cagcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCAc      >  2:655559/1‑150 (MQ=255)
                                                                                                                                                 cagcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCAc      >  2:2771629/1‑150 (MQ=255)
                                                                                                                                                    cagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACttt   >  2:5189361/1‑150 (MQ=255)
                                                                                                                                                     agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTg  >  1:2098041/1‑150 (MQ=255)
                                                                                                                                                     agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTg  >  1:2098493/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTG  >  CP009273/777026‑777324

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