Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,136,421 A→G 7.8% D3D (GAT→GAC udk ← uridine/cytidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,136,4210AG 7.8% 1179.1 / 66.9 357D3D (GAT→GACudkuridine/cytidine kinase
Reads supporting (aligned to +/- strand):  ref base A (186/143);  new base G (15/13);  total (201/156)
Fisher's exact test for biased strand distribution p-value = 8.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 359 total aligned reads displayed.

TGATCTTTGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGACAAAGAATTCC  >  CP009273/2136272‑2136570
                                                                                                                                                     |                                                                                                                                                     
  aTCTTTGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATc                                                                                                                                                     <  2:4193/150‑1 (MQ=255)
       tGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCa                                                                                                                                                <  1:732443/150‑1 (MQ=255)
        gTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCat                                                                                                                                               <  1:4162535/150‑1 (MQ=255)
         tAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCata                                                                                                                                              >  2:4297686/1‑150 (MQ=255)
             aGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATAtt                                                                                                                                          <  1:1329465/150‑1 (MQ=255)
                aGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAg                                                                                                                                       >  1:1472304/1‑150 (MQ=255)
                aGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAg                                                                                                                                       >  1:1973807/1‑150 (MQ=255)
                    ttCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGAc                                                                                                                                   >  2:625957/1‑150 (MQ=255)
                      cGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCt                                                                                                                                 <  2:1070750/150‑1 (MQ=255)
                        ggAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGa                                                                                                                               >  2:906485/1‑150 (MQ=255)
                         gAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGAt                                                                                                                              <  1:2376934/150‑1 (MQ=255)
                               cGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCt                                                                                                                        >  2:4962991/1‑150 (MQ=255)
                                 ccGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTgg                                                                                                                      <  1:3632067/150‑1 (MQ=255)
                                   gATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGAt                                                                                                                    >  2:4584064/1‑150 (MQ=255)
                                    aTGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATg                                                                                                                   <  1:4842071/150‑1 (MQ=255)
                                     tgtgTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATGCCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGt                                                                                                                  >  1:1157150/1‑150 (MQ=255)
                                     tgtgTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGt                                                                                                                  <  2:3901493/150‑1 (MQ=255)
                                       tgtTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTa                                                                                                                >  1:1001083/1‑150 (MQ=255)
                                       tgtTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTa                                                                                                                <  2:4055514/150‑1 (MQ=255)
                                         ttCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATCATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAgg                                                                                                              >  2:3335764/1‑150 (MQ=255)
                                         ttCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAgg                                                                                                              <  1:2467419/150‑1 (MQ=255)
                                          tcatcaCCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGa                                                                                                             >  1:4571561/1‑150 (MQ=255)
                                           catcaCCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGaa                                                                                                            >  2:5016706/1‑150 (MQ=255)
                                             tcaCCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAgg                                                                                                          >  1:5132942/1‑150 (MQ=255)
                                             tcaCCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAgg                                                                                                          >  2:4271399/1‑150 (MQ=255)
                                             tcaCCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCCGATTAACCTGGATGTTAGGAAgg                                                                                                          <  2:4749876/150‑1 (MQ=255)
                                             tcaCCGACGTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAgg                                                                                                          <  2:698094/150‑1 (MQ=255)
                                              caCCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAggg                                                                                                         >  2:4011178/1‑150 (MQ=255)
                                                ccGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAggggg                                                                                                       <  1:1718911/150‑1 (MQ=255)
                                                 cGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGc                                                                                                      >  1:1915396/1‑150 (MQ=255)
                                                  gACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGcg                                                                                                     <  1:1766929/150‑1 (MQ=255)
                                                   aCTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGcga                                                                                                    >  2:3152133/1‑150 (MQ=255)
                                                   aCTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGcga                                                                                                    >  2:3868087/1‑150 (MQ=255)
                                                       gCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGaa                                                                                                >  1:3166274/1‑150 (MQ=255)
                                                          cACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAgcg                                                                                             >  1:2190603/1‑150 (MQ=255)
                                                           aCGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAgcgc                                                                                            <  1:2556834/150‑1 (MQ=255)
                                                            cGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGcc                                                                                           >  1:531334/1‑150 (MQ=255)
                                                             gCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCa                                                                                          <  1:4916912/150‑1 (MQ=255)
                                                              cAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCaa                                                                                         <  1:3837044/150‑1 (MQ=255)
                                                                aTTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAAc                                                                                       >  2:783965/1‑150 (MQ=255)
                                                                    aCGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCgg                                                                                   >  2:1912982/1‑150 (MQ=255)
                                                                    aCGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCgg                                                                                   >  2:125174/1‑150 (MQ=255)
                                                                          aaGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATta                                                                             >  2:1171286/1‑150 (MQ=255)
                                                                           aGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATtat                                                                            <  1:1050250/150‑1 (MQ=255)
                                                                             ggTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAg                                                                          >  1:4757145/1‑150 (MQ=255)
                                                                              gTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAgg                                                                         >  1:4411459/1‑150 (MQ=255)
                                                                                 cTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGAt                                                                      <  2:872971/150‑1 (MQ=255)
                                                                                 cTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGAt                                                                      <  2:2146522/150‑1 (MQ=255)
                                                                                  tGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGAtt                                                                     <  1:881936/150‑1 (MQ=255)
                                                                                  tGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTCAGTCATATATATAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGAtt                                                                     <  1:3914245/150‑1 (MQ=255)
                                                                                  tGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGAtt                                                                     <  1:3781183/150‑1 (MQ=255)
                                                                                  tGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAAGGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGAtt                                                                     <  2:4851610/150‑1 (MQ=255)
                                                                                   ggCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGAttt                                                                    <  1:3515043/150‑1 (MQ=255)
                                                                                    gCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTc                                                                   <  1:3335873/150‑1 (MQ=255)
                                                                                     cAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCa                                                                  >  1:2767792/1‑150 (MQ=255)
                                                                                     cAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCa                                                                  >  1:4239911/1‑150 (MQ=255)
                                                                                       aTAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATc                                                                >  1:4976933/1‑150 (MQ=255)
                                                                                        tAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATcc                                                               >  2:3274329/1‑150 (MQ=255)
                                                                                        tAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATcc                                                               <  2:4524967/150‑1 (MQ=255)
                                                                                          aGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGt                                                             >  1:3736946/1‑150 (MQ=255)
                                                                                           gACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTc                                                            >  2:5172258/1‑150 (MQ=255)
                                                                                            aCTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATACTGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCt                                                           <  2:2609777/150‑1 (MQ=255)
                                                                                            aCTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCt                                                           >  1:660584/1‑150 (MQ=255)
                                                                                            aCTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCt                                                           >  2:4133306/1‑150 (MQ=255)
                                                                                              tctTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGa                                                         <  1:1434011/150‑1 (MQ=255)
                                                                                                   ccGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCa                                                    >  2:177643/1‑150 (MQ=255)
                                                                                                    cGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAg                                                   >  1:428369/1‑150 (MQ=255)
                                                                                                     ggAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGt                                                  <  1:3387983/150‑1 (MQ=255)
                                                                                                      gAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTc                                                 >  2:3497429/1‑150 (MQ=255)
                                                                                                          ccGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAAt                                             <  1:2892345/150‑1 (MQ=255)
                                                                                                            gATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAg                                           <  2:1570588/150‑1 (MQ=255)
                                                                                                            gATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAg                                           >  2:347148/1‑150 (MQ=255)
                                                                                                              tGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCg                                         <  1:408942/150‑1 (MQ=255)
                                                                                                              tGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCg                                         <  1:2513417/150‑1 (MQ=255)
                                                                                                                cgcCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGtt                                       <  2:3429634/150‑1 (MQ=255)
                                                                                                                cgcCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGtt                                       <  2:4680267/150‑1 (MQ=255)
                                                                                                                     gCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCcgc                                  <  1:2337594/150‑1 (MQ=255)
                                                                                                                      cGATACCGATAATGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCcgcg                                 <  1:3016547/150‑1 (MQ=255)
                                                                                                                        aTACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCcgcgcg                               >  1:2414784/1‑150 (MQ=255)
                                                                                                                        aTACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCcgcgcg                               >  2:2315024/1‑150 (MQ=255)
                                                                                                                          aCCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCt                             <  1:3936785/150‑1 (MQ=255)
                                                                                                                            cGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCa                           >  1:808485/1‑150 (MQ=255)
                                                                                                                               tAATGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCagag                        <  2:2396591/150‑1 (MQ=255)
                                                                                                                               tAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCagag                        >  2:2013761/1‑150 (MQ=255)
                                                                                                                                 aTGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGtt                      >  2:2513942/1‑150 (MQ=255)
                                                                                                                                 aTGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGtt                      >  2:367384/1‑150 (MQ=255)
                                                                                                                                  tGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTa                     >  1:5038274/1‑150 (MQ=255)
                                                                                                                                  tGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATGAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTa                     <  2:3115184/150‑1 (MQ=255)
                                                                                                                                        aCTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAAttgtt               >  1:3145422/1‑150 (MQ=255)
                                                                                                                                         cTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAAttgttg              >  2:1145920/1‑150 (MQ=255)
                                                                                                                                          tgatgaGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGa             <  1:4700074/150‑1 (MQ=255)
                                                                                                                                           gatgaGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGAc            >  2:2960631/1‑150 (MQ=255)
                                                                                                                                           gatgaGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATAGTTGAc            >  2:2667172/1‑150 (MQ=255)
                                                                                                                                            atgaGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGACa           >  2:3978312/1‑150 (MQ=255)
                                                                                                                                             tgaGACTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGACaa          >  2:1729400/1‑150 (MQ=255)
                                                                                                                                              gagaCTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGACaaa         >  2:2592858/1‑150 (MQ=255)
                                                                                                                                                gaCTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGACAAAGa       <  1:1595873/150‑1 (MQ=255)
                                                                                                                                                  cTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGACAAAGAAt     <  1:1806327/150‑1 (MQ=255)
                                                                                                                                                   tGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGACAAAGAAtt    >  2:1851278/1‑150 (MQ=255)
                                                                                                                                                    gATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGACAAAGAATTc   >  2:2368481/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGATCTTTGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGACAAAGAATTCC  >  CP009273/2136272‑2136570

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: