Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,798,582 C→A 9.3% Q243H (CAG→CAT waaP ← kinase that phosphorylates core heptose of lipopolysaccharide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,798,5820CA 9.3% 767.9 / 36.0 225Q243H (CAG→CATwaaPkinase that phosphorylates core heptose of lipopolysaccharide
Reads supporting (aligned to +/- strand):  ref base C (120/84);  new base A (8/13);  total (128/97)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

Only 100 of 226 total aligned reads displayed.

TAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGC  >  CP009273/3798435‑3798728
                                                                                                                                                   |                                                                                                                                                  
tAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTc                                                                                                                                                    >  2:2953515/1‑148 (MQ=255)
  aTTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTg                                                                                                                                                  >  2:2761955/1‑148 (MQ=255)
        aaTCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAag                                                                                                                                            >  2:2578435/1‑148 (MQ=255)
        aaTCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCGTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAag                                                                                                                                            >  2:2579335/1‑148 (MQ=255)
          tCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAagag                                                                                                                                          >  1:454403/1‑148 (MQ=255)
              cAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACaa                                                                                                                                      >  2:3146076/1‑148 (MQ=255)
               aaCCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAat                                                                                                                                     >  1:2266797/1‑148 (MQ=255)
                 cccAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCCTGTTCAATAAATTCTGTTCAAGAGACAatat                                                                                                                                   >  2:2290670/1‑148 (MQ=255)
                 cccAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAatat                                                                                                                                   >  2:596282/1‑148 (MQ=255)
                   cAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATAttt                                                                                                                                 >  1:2919991/1‑148 (MQ=255)
                    aaTTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATAtttt                                                                                                                                <  2:3208996/148‑1 (MQ=255)
                          taaaataaaTAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTaa                                                                                                                          >  2:592020/1‑148 (MQ=255)
                           aaaataaaTAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTaaa                                                                                                                         >  2:3327045/1‑148 (MQ=255)
                           aaaataaaTAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTaaa                                                                                                                         >  1:1418884/1‑148 (MQ=255)
                               taaataGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCa                                                                                                                     <  1:3429913/148‑1 (MQ=255)
                               taaataGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCa                                                                                                                     <  1:2359768/148‑1 (MQ=255)
                                aaataGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCAt                                                                                                                    >  2:1952729/1‑148 (MQ=255)
                                  ataGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATcc                                                                                                                  >  1:2435020/1‑148 (MQ=255)
                                    aGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCa                                                                                                                <  1:1544069/148‑1 (MQ=255)
                                    aGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCa                                                                                                                <  2:1575913/148‑1 (MQ=255)
                                      tCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCaaa                                                                                                              <  2:3122661/148‑1 (MQ=255)
                                      tCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCaaa                                                                                                              >  2:181406/1‑148 (MQ=255)
                                       cATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCaaaa                                                                                                             >  2:3489102/1‑148 (MQ=255)
                                         taataaAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATa                                                                                                           >  1:2489326/1‑148 (MQ=255)
                                         taataaAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATa                                                                                                           <  1:434254/148‑1 (MQ=255)
                                         taataaAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATa                                                                                                           >  2:1435606/1‑148 (MQ=255)
                                         taataaAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATa                                                                                                           >  2:2243706/1‑148 (MQ=255)
                                         taataaAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATa                                                                                                           >  2:456931/1‑148 (MQ=255)
                                         taataaAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATa                                                                                                           >  1:2243554/1‑148 (MQ=255)
                                           ataaAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATaaa                                                                                                         <  1:723035/148‑1 (MQ=255)
                                            taaAGTTAGTTCCAGTACATACTTATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAAc                                                                                                        <  2:2334472/148‑1 (MQ=255)
                                             aaaGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAAcc                                                                                                       <  1:2012488/148‑1 (MQ=255)
                                             aaaGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAAcc                                                                                                       >  2:1228589/1‑148 (MQ=255)
                                                 ttagttCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCa                                                                                                   >  1:46807/1‑148 (MQ=255)
                                                 ttagttCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCa                                                                                                   <  1:3103613/148‑1 (MQ=255)
                                                    gttCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATaa                                                                                                >  2:3440425/1‑148 (MQ=255)
                                                    gttCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATaa                                                                                                >  1:1162239/1‑148 (MQ=255)
                                                    gttCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATaa                                                                                                >  2:306728/1‑148 (MQ=255)
                                                    gttCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATaa                                                                                                >  2:2731773/1‑148 (MQ=255)
                                                      tCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAAt                                                                                              >  2:426381/1‑148 (MQ=255)
                                                      tCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAAt                                                                                              <  1:2871591/148‑1 (MQ=255)
                                                            aCATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAg                                                                                        >  1:1684554/1‑148 (MQ=255)
                                                                  taataaATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATAtct                                                                                  <  2:3500588/148‑1 (MQ=255)
                                                                     taaataTTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCtt                                                                               <  2:61251/148‑1 (MQ=255)
                                                                      aaataTTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCttt                                                                              >  1:21809/1‑148 (MQ=255)
                                                                       aataTTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTg                                                                             >  1:634208/1‑148 (MQ=255)
                                                                            ttttATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCa                                                                        >  1:509834/1‑148 (MQ=255)
                                                                              ttATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAgg                                                                      <  2:3180779/148‑1 (MQ=255)
                                                                               tataATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGc                                                                     >  2:892910/1‑148 (MQ=255)
                                                                                  aaTCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCaa                                                                  >  1:985409/1‑148 (MQ=255)
                                                                                    tCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAata                                                                >  1:946800/1‑148 (MQ=255)
                                                                                    tCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAata                                                                >  2:192569/1‑148 (MQ=255)
                                                                                     ccTTTGCTTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAatat                                                               <  1:2953515/148‑1 (MQ=255)
                                                                                     ccTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAatat                                                               >  1:1067030/1‑148 (MQ=255)
                                                                                     ccTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAatat                                                               >  1:382468/1‑148 (MQ=255)
                                                                                     ccTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAatat                                                               >  1:381563/1‑148 (MQ=255)
                                                                                      cTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATAtt                                                              >  1:356309/1‑148 (MQ=255)
                                                                                       tttGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTc                                                             >  2:303836/1‑148 (MQ=255)
                                                                                        ttGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCa                                                            >  1:1694058/1‑148 (MQ=255)
                                                                                        ttGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCa                                                            <  2:1709250/148‑1 (MQ=255)
                                                                                           cGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATtg                                                         <  1:1603769/148‑1 (MQ=255)
                                                                                              tgtgTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTc                                                             <  1:3232131/141‑1 (MQ=255)
                                                                                              tgtgTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATtgatg                                                      <  2:201793/148‑1 (MQ=255)
                                                                                                     gCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAata                                               >  1:2546385/1‑148 (MQ=255)
                                                                                                        ctTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATaa                                            <  2:3380981/148‑1 (MQ=255)
                                                                                                        ctTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATaa                                            <  1:662880/148‑1 (MQ=255)
                                                                                                        ctTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATaa                                            >  1:3452629/1‑148 (MQ=255)
                                                                                                            aaTACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCa                                        <  1:3145431/148‑1 (MQ=255)
                                                                                                             aTACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCaa                                       <  1:229389/148‑1 (MQ=255)
                                                                                                              tACGCTCGGCCTTAACTCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAAt                                      >  1:1342183/1‑148 (MQ=255)
                                                                                                                 gCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAg                                   <  1:1013539/148‑1 (MQ=255)
                                                                                                                  cTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAgg                                  <  2:931815/148‑1 (MQ=255)
                                                                                                                   tCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGt                                 <  1:3329234/148‑1 (MQ=255)
                                                                                                                     ggCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCt                               <  1:1940776/148‑1 (MQ=255)
                                                                                                                       ccTTAACTCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCttt                             <  2:2782659/148‑1 (MQ=255)
                                                                                                                       ccTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCttt                             <  1:3171529/148‑1 (MQ=255)
                                                                                                                       ccTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCttt                             >  1:32699/1‑148 (MQ=255)
                                                                                                                          tAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATc                          >  2:801354/1‑148 (MQ=255)
                                                                                                                           aaCGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATcg                         <  1:596282/148‑1 (MQ=255)
                                                                                                                           aaCGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATcg                         <  1:2731773/148‑1 (MQ=255)
                                                                                                                            aCGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATcgc                        >  2:1086662/1‑148 (MQ=255)
                                                                                                                             cGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATcgcg                       <  2:208156/148‑1 (MQ=255)
                                                                                                                              gCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCgcgc                      <  1:963738/148‑1 (MQ=255)
                                                                                                                              gCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCgcgc                      <  2:2484180/148‑1 (MQ=255)
                                                                                                                              gCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCAAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCgcgc                      <  1:1165234/148‑1 (MQ=255)
                                                                                                                              gCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCgcgc                      <  2:2243554/148‑1 (MQ=255)
                                                                                                                                 tGCCATTTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAg                   <  2:3375135/148‑1 (MQ=255)
                                                                                                                                   ccATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAgcg                 >  2:689758/1‑148 (MQ=255)
                                                                                                                                     aTGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAgcggc               <  1:851714/148‑1 (MQ=255)
                                                                                                                                      tGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAgcggcg              >  2:1832323/1‑148 (MQ=255)
                                                                                                                                       gTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGgcgc             >  2:2299600/1‑148 (MQ=255)
                                                                                                                                        ttCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGcgcg            >  1:3442763/1‑148 (MQ=255)
                                                                                                                                         tCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCgg           >  2:147445/1‑148 (MQ=255)
                                                                                                                                          cAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGt          >  2:3060236/1‑148 (MQ=255)
                                                                                                                                             tAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGTAAACTATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACt       <  1:3489102/148‑1 (MQ=255)
                                                                                                                                             tAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACt       >  1:147924/1‑148 (MQ=255)
                                                                                                                                              aaaTTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACtt      >  1:3007163/1‑148 (MQ=255)
                                                                                                                                                  tCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGc  >  1:2043278/1‑148 (MQ=255)
                                                                                                                                                  tCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGc  >  1:195048/1‑148 (MQ=255)
                                                                                                                                                  tCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGc  >  1:1780487/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGC  >  CP009273/3798435‑3798728

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: