Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,477,570 G→A 7.4% D85N (GAC→AAC)  hrpA → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,477,5700GA 7.4% 956.3 / 42.9 289D85N (GAC→AAC) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (128/135);  new base A (11/10);  total (139/150)
Fisher's exact test for biased strand distribution p-value = 8.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 296 total aligned reads displayed.

GAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGC  >  CP009273/1477422‑1477718
                                                                                                                                                     |                                                                                                                                                    
gaagGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAgg                                                                                                                                                      >  2:4012657/1‑149 (MQ=255)
 aagGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGa                                                                                                                                                     >  1:1440395/1‑149 (MQ=255)
   ggTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACa                                                                                                                                                   <  2:2618808/149‑1 (MQ=255)
      tAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtc                                                                                                                                                <  1:4068855/149‑1 (MQ=255)
       aaaaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtct                                                                                                                                               <  1:1849325/149‑1 (MQ=255)
        aaaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtctc                                                                                                                                              <  2:3788794/149‑1 (MQ=255)
          aaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGa                                                                                                                                            <  1:4234337/149‑1 (MQ=255)
          aaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGa                                                                                                                                            >  1:2091247/1‑149 (MQ=255)
          aaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGa                                                                                                                                            <  1:1597530/149‑1 (MQ=255)
           aTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGaa                                                                                                                                           <  1:3700373/149‑1 (MQ=255)
               tGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGa                                                                                                                                       >  1:1963141/1‑149 (MQ=255)
                gATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGAt                                                                                                                                      <  2:3494383/149‑1 (MQ=255)
                 aTGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGAtt                                                                                                                                     <  1:527537/149‑1 (MQ=255)
                 aTGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGAtt                                                                                                                                     <  1:2708832/149‑1 (MQ=255)
                 aTGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGAtt                                                                                                                                     <  1:4746641/149‑1 (MQ=255)
                  tGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTc                                                                                                                                    >  1:3427216/1‑149 (MQ=255)
                  tGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTc                                                                                                                                    <  1:1360308/149‑1 (MQ=255)
                     acaacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGACTCGTg                                                                                                                                 >  1:1795001/1‑149 (MQ=255)
                      caacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGa                                                                                                                                >  2:3662508/1‑149 (MQ=255)
                      caacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGa                                                                                                                                >  2:3418337/1‑149 (MQ=255)
                        acaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATc                                                                                                                              <  1:3055206/149‑1 (MQ=255)
                          aGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCAc                                                                                                                            <  2:3170101/149‑1 (MQ=255)
                           ggCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCAcc                                                                                                                           >  1:4787936/1‑149 (MQ=255)
                               aTTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAggt                                                                                                                       <  2:3142417/149‑1 (MQ=255)
                                ttttCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAggtg                                                                                                                      <  1:3217324/149‑1 (MQ=255)
                                 tttCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAggtgg                                                                                                                     >  1:400694/1‑149 (MQ=255)
                                  ttCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAggtggt                                                                                                                    <  1:4427662/149‑1 (MQ=255)
                                   tCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGgtggtg                                                                                                                   <  2:1201995/149‑1 (MQ=255)
                                    ccAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGa                                                                                                                  <  2:1174110/149‑1 (MQ=255)
                                     cAGGAGATGGCGAAATAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAt                                                                                                                 <  1:35230/149‑1 (MQ=255)
                                        gagaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcgt                                                                                                              <  2:2341716/149‑1 (MQ=255)
                                        gagaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcgt                                                                                                              <  2:979594/149‑1 (MQ=255)
                                        gagaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcgt                                                                                                              >  1:1463881/1‑149 (MQ=255)
                                         agaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcgtc                                                                                                             >  2:2785771/1‑149 (MQ=255)
                                         agaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcgtc                                                                                                             >  2:2968930/1‑149 (MQ=255)
                                            tGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGcc                                                                                                          >  1:1833384/1‑149 (MQ=255)
                                             ggCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCg                                                                                                         >  1:1117335/1‑149 (MQ=255)
                                                gAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGAGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCgggg                                                                                                      <  2:2694285/149‑1 (MQ=255)
                                                 aaaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACACTTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGa                                                                                                     >  1:3503032/1‑149 (MQ=255)
                                                 aaaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGa                                                                                                     >  1:4398467/1‑149 (MQ=255)
                                                  aaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGaa                                                                                                    >  2:2467066/1‑149 (MQ=255)
                                                     agaTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACg                                                                                                 >  2:4716359/1‑149 (MQ=255)
                                                         tGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGtt                                                                                             <  2:114849/149‑1 (MQ=255)
                                                          gACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTc                                                                                            <  2:1885574/149‑1 (MQ=255)
                                                          gACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTGACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTc                                                                                            <  2:4275503/149‑1 (MQ=255)
                                                             cAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTgg                                                                                         <  1:3824868/149‑1 (MQ=255)
                                                              aGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGt                                                                                        <  1:2341039/149‑1 (MQ=255)
                                                                gcggcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTaa                                                                                      <  2:3356455/149‑1 (MQ=255)
                                                                  ggcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTaaaa                                                                                    <  1:1008689/149‑1 (MQ=255)
                                                                  ggcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTaaaa                                                                                    <  2:497065/149‑1 (MQ=255)
                                                                     aGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGa                                                                                 >  2:2836689/1‑149 (MQ=255)
                                                                     aGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGa                                                                                 >  2:4064794/1‑149 (MQ=255)
                                                                        gAAAGTCCTGCTGCGTGAAGCGGCACGACCGTAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTc                                                                              <  2:4192425/149‑1 (MQ=255)
                                                                             tCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAAGTGACCGGTTAGTCAGAACAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTa                                                                         >  1:326515/1‑149 (MQ=255)
                                                                               ctgctgCTTTAAGCGGCACGACCGGAAATTACTTATCATGACATTTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTTATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGTTAAAACGACTCAGTTAcc                                                                       <  2:4583730/149‑1 (MQ=255)
                                                                               ctgctgCGTGATGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGTAAACGGGTTCTGGTAAAACGACTCAGTTAcc                                                                       <  2:4591694/149‑1 (MQ=255)
                                                                               ctgatgCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTAcc                                                                       >  2:4606463/1‑149 (MQ=255)
                                                                                tgctgcGTGAAGCGGCACTACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTTATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCg                                                                      <  2:1978591/149‑1 (MQ=255)
                                                                                    gcGTTAAGCGGCACGACCGGAAATTACTTATCCTGAAAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGaaaa                                                                  <  2:2213174/149‑1 (MQ=255)
                                                                                    gcGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGaaaa                                                                  >  1:1670315/1‑149 (MQ=255)
                                                                                    gcGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGaaaa                                                                  <  2:2370073/149‑1 (MQ=255)
                                                                                      gTGACGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACTGGACATTCTCGAAGCGATTCGTTATCACCAGGTGGGGATCGTCGCCGGGGAAACGGGTTCGGGTAAAACGACTCAGTTACCGAAAATc                                                                <  2:4145833/149‑1 (MQ=255)
                                                                                         aaGCGGCACGCCCGGAAATTACTGCTCCTGACATTTTACCGGTTTGTCAGACAACACTGGACATTCTCGACGCGCTTCGTGATCACCTGGGGGTGATCGGCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGt                                                             <  2:3610267/149‑1 (MQ=255)
                                                                                         aaGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGt                                                             >  1:2228710/1‑149 (MQ=255)
                                                                                            cGGCACGACCGGAAATTACTTAACCGGACAATTTACCGGTTAGTCAGAAAAAACTGGACATTCACGACGCGTTTCGTGATCACCCGGTGGTGATCGTCGCCGGGGAAAGGGGTTCGGGTAAAACGACTCAGTTACCGAAAATCTGTATg                                                          <  1:650902/149‑1 (MQ=255)
                                                                                             ggCCCGACCGGAAATTACTTATCCGGACAAGTTCCCGGTGTGTCAGAACACACCGGACCGTCGGGATGCGATTCGTGCTCACCAGGGGGTGAGCGTCGCCGGGGATCCGGGTTCGGGTAAAACGACTCAGTTACCGAAACTCTGTATgg                                                         <  2:784566/149‑1 (MQ=255)
                                                                                                   aCCGTAAATTACTTATCCTTAAATTTTCCCGGTTAGTCATAAAAAACTGTACTTTCTCGACGCGATTCGTGATCACCCGGTGGTTATCGTCGCCGGGGAAAGGGGTTCTGTTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTgg                                                   <  1:872132/149‑1 (MQ=255)
                                                                                                         aaTTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGcgcg                                             >  1:868946/1‑149 (MQ=255)
                                                                                                          tttACTTATCCTGAAAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCCGGTGGTGATCGTCGCCGGGTAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCgg                                            <  1:4577399/148‑1 (MQ=255)
                                                                                                          aTTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCgg                                            >  2:3012358/1‑149 (MQ=255)
                                                                                                              tttATCATGAAAATTTACCGGTTATTCAGAAAAAACAGGACATGCTCTAAGCGATTCGTGATCACCAGGTGGTGTTGGTCGCCGGGGAAAGGGGTTCTGTTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGAtt                                        <  2:4492143/148‑1 (MQ=255)
                                                                                                               ttATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTa                                       >  1:964934/1‑149 (MQ=255)
                                                                                                                 aTCCTGACAATTTACCGGTTAGTCAGAAAAAACTGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCGGGGGAAAGGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTaaa                                     <  1:1909000/149‑1 (MQ=255)
                                                                                                                 aTCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTaaa                                     >  1:4659573/1‑149 (MQ=255)
                                                                                                                 aTCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTaaa                                     >  2:3957504/1‑149 (MQ=255)
                                                                                                                 aTCCTGAAAATTTACCGTTTAGTCATAAAAAACAGGACATTCTCGAAGCTATTCGTGATCACCTGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTaaa                                     <  1:2489811/149‑1 (MQ=255)
                                                                                                                  tCCTGACATTTTACCGGTTAGTCAGAAAAAACCGGACATTCTCGAAGCGATTCGTGATCACCCGGTGGTGATCGTCGCCGGGGAAAGGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAg                                    <  2:2029876/149‑1 (MQ=255)
                                                                                                                  tCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAg                                    >  2:455678/1‑149 (MQ=255)
                                                                                                                  tCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAg                                    >  1:2579211/1‑149 (MQ=255)
                                                                                                                     tGAAAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGCTTCGTGATCACCAGGTGGGGATCGTCGCCGGGTAACCGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGAc                                 <  2:3262200/149‑1 (MQ=255)
                                                                                                                      gaaaaTTTACGGGTTAGACAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCTGGTGGTTATCGTCGCCGGGTAAAGGGGTTCTGTAAACACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACt                                <  1:1760109/146‑1 (MQ=255)
                                                                                                                         attttACCGGTTCGTCAGAAAAAACAGGCCATTCTCGAAGCGATTCGTGATCACCCGGTGGGGATCGTCGCCGGGGAACCGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGAt                             <  1:288337/147‑1 (MQ=255)
                                                                                                                         aaTTTACCGTTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCGGTATGGAGCTGGGGCGCGGGATTAAAGGACTGAt                             <  2:61462/149‑1 (MQ=255)
                                                                                                                          aTTTTCCGGTTTAGTCATAAAAAACCGGAACTGCTCTACGCGAGTCGTGATCACCTGGGGGTGATGGGCGCCGGGGAACGGGGTTGTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGAt                             <  1:2042489/149‑1 (MQ=255)
                                                                                                                          aTTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATc                            >  1:3221444/1‑149 (MQ=255)
                                                                                                                           tttACCGGTTATTCACAAAAAACAGGACATTCCCGATGCGATTCGTGATCACCAGGTGGTGATCGTCGCGGGGTAAACGGGTAGTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCg                           <  1:533716/149‑1 (MQ=255)
                                                                                                                           tttACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCg                           >  1:1054858/1‑149 (MQ=255)
                                                                                                                                      aGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAg                >  2:216009/1‑149 (MQ=255)
                                                                                                                                        tCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGc               >  2:3126035/1‑149 (MQ=255)
                                                                                                                                           gAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCcgc            >  1:1380158/1‑149 (MQ=255)
                                                                                                                                            aaaaaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCcgcg           >  1:2182315/1‑149 (MQ=255)
                                                                                                                                             aaaaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgt          >  2:1365024/1‑149 (MQ=255)
                                                                                                                                              aaaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTAACGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgtc         >  1:3932975/1‑149 (MQ=255)
                                                                                                                                                aaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgtcgt       >  2:3457870/1‑149 (MQ=255)
                                                                                                                                                 aCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCgtcgtc      <  2:1827839/149‑1 (MQ=255)
                                                                                                                                                  cAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCt     >  2:292213/1‑149 (MQ=255)
                                                                                                                                                  cAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCt     >  2:1745703/1‑149 (MQ=255)
                                                                                                                                                   aGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTg    >  2:3575528/1‑149 (MQ=255)
                                                                                                                                                   aGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTg    >  2:1132629/1‑149 (MQ=255)
                                                                                                                                                     gACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTggc  >  2:4651221/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGC  >  CP009273/1477422‑1477718

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: