Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,570,897 G→T 10.5% S128Y (TCT→TAT)  hsdS ← specificity determinant for hsdM and hsdR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,570,8970GT10.5% 977.7 / 55.8 281S128Y (TCT→TAT) hsdSspecificity determinant for hsdM and hsdR
Reads supporting (aligned to +/- strand):  ref base G (129/122);  new base T (12/18);  total (141/140)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 297 total aligned reads displayed.

CGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAAC  >  CP009273/4570749‑4571045
                                                                                                                                                    |                                                                                                                                                    
cgcCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGagaag                                                                                                                                                      <  2:771684/149‑1 (MQ=255)
  ccAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGAt                                                                                                                                                    <  1:3555158/149‑1 (MQ=255)
  ccAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGAt                                                                                                                                                    <  1:1561885/149‑1 (MQ=255)
     gcagcGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAATAtttt                                                                                                                                                 <  2:2378840/149‑1 (MQ=255)
       agcGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTtgt                                                                                                                                               <  2:196508/149‑1 (MQ=255)
            aTCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAATATTTTGTGAAAt                                                                                                                                          >  1:1527312/1‑149 (MQ=255)
                 gTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCa                                                                                                                                     >  1:2229296/1‑149 (MQ=255)
                  tttttCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCaa                                                                                                                                    <  1:4051433/149‑1 (MQ=255)
                   ttttCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAAt                                                                                                                                   <  1:1903746/149‑1 (MQ=255)
                      tCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATaaa                                                                                                                                <  1:4305488/149‑1 (MQ=255)
                        aGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAATATTTTGTGAAATGAGCAATAAAAc                                                                                                                              <  1:2687426/149‑1 (MQ=255)
                         gCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAAcc                                                                                                                             >  2:1125279/1‑149 (MQ=255)
                         gCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAAcc                                                                                                                             >  1:321984/1‑149 (MQ=255)
                           gatgatTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAg                                                                                                                           >  1:4554338/1‑149 (MQ=255)
                           gatgatTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAg                                                                                                                           >  2:2158187/1‑149 (MQ=255)
                            atgatTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGa                                                                                                                          >  1:4318586/1‑149 (MQ=255)
                               atTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGaaaa                                                                                                                       <  1:3296252/149‑1 (MQ=255)
                                 tttttGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAata                                                                                                                     >  1:654327/1‑149 (MQ=255)
                                  ttttGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAatat                                                                                                                    <  2:2786892/149‑1 (MQ=255)
                                  ttttGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAatat                                                                                                                    <  2:4439195/149‑1 (MQ=255)
                                          ggCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAATATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGttttt                                                                                                            <  1:3719819/149‑1 (MQ=255)
                                          ggCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGttttt                                                                                                            <  2:4684628/149‑1 (MQ=255)
                                          ggCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGttttt                                                                                                            >  2:2297418/1‑149 (MQ=255)
                                          ggCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGttttt                                                                                                            <  1:2195240/149‑1 (MQ=255)
                                            cAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCa                                                                                                          <  2:3801682/149‑1 (MQ=255)
                                             aaGTGGTGTGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAg                                                                                                         <  2:2927109/149‑1 (MQ=255)
                                             aaGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAg                                                                                                         <  1:4165021/149‑1 (MQ=255)
                                             aaGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAg                                                                                                         <  2:503567/149‑1 (MQ=255)
                                             aaGTGGTGGGATTGGTATATTGATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAg                                                                                                         <  1:4460507/149‑1 (MQ=255)
                                              aGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAgg                                                                                                        >  2:1272167/1‑149 (MQ=255)
                                               gtggtgGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGa                                                                                                       <  2:3978692/149‑1 (MQ=255)
                                                  gtggGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGt                                                                                                    >  1:4495369/1‑149 (MQ=255)
                                                  gtggGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGt                                                                                                    >  2:1600671/1‑149 (MQ=255)
                                                  gtggGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGt                                                                                                    >  1:3539922/1‑149 (MQ=255)
                                                     ggATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAAt                                                                                                 <  1:699796/149‑1 (MQ=255)
                                                      gATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATa                                                                                                >  1:3381612/1‑149 (MQ=255)
                                                        ttGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATaca                                                                                              >  1:3025957/1‑149 (MQ=255)
                                                             atatTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCa                                                                                         >  1:2870252/1‑149 (MQ=255)
                                                             atatTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTCAGGACGTAATACACCGCaa                                                                                        >  2:3483085/1‑149 (MQ=255)
                                                              tatTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCaa                                                                                        >  1:3709367/1‑149 (MQ=255)
                                                              tatTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCaa                                                                                        >  1:8524/1‑149 (MQ=255)
                                                              tatTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCaa                                                                                        >  2:3621966/1‑149 (MQ=255)
                                                                     caaatcaaaGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATgcgc                                                                                 <  1:2017532/149‑1 (MQ=255)
                                                                     caaatcaaaGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATgcgc                                                                                 <  2:2852899/149‑1 (MQ=255)
                                                                        atcaaaGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGa                                                                              >  1:576239/1‑149 (MQ=255)
                                                                         tcaaaGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGaa                                                                             >  2:3054560/1‑149 (MQ=255)
                                                                            aaGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACt                                                                          >  1:108214/1‑149 (MQ=255)
                                                                              gCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTAc                                                                        >  1:1706402/1‑149 (MQ=255)
                                                                              gCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTAc                                                                        <  2:4282714/149‑1 (MQ=255)
                                                                               cTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACa                                                                       <  1:1000491/149‑1 (MQ=255)
                                                                                ttGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACAt                                                                      >  2:171973/1‑149 (MQ=255)
                                                                                 tGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACAtt                                                                     <  1:1220328/149‑1 (MQ=255)
                                                                                 tGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACAtt                                                                     <  2:1420776/149‑1 (MQ=255)
                                                                                      gCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAAt                                                                >  2:1772012/1‑149 (MQ=255)
                                                                                       cTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATg                                                               <  2:2741278/149‑1 (MQ=255)
                                                                                        ttAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATgg                                                              <  1:3592908/149‑1 (MQ=255)
                                                                                        ttAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATgg                                                              <  1:522971/149‑1 (MQ=255)
                                                                                         tAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGt                                                             <  1:1865248/149‑1 (MQ=255)
                                                                                          aaTATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGATCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTa                                                            >  1:278881/1‑149 (MQ=255)
                                                                                           atatTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGGGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAg                                                           >  1:758166/1‑149 (MQ=255)
                                                                                           atatTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAg                                                           >  1:2060586/1‑149 (MQ=255)
                                                                                               tattaATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGc                                                       <  1:1757368/149‑1 (MQ=255)
                                                                                                 ttaATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAATATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTg                                                     <  1:131237/149‑1 (MQ=255)
                                                                                                 ttaATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTg                                                     >  2:388886/1‑149 (MQ=255)
                                                                                                 ttaATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTg                                                     >  2:2005276/1‑149 (MQ=255)
                                                                                                   aaTATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTGATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGAt                                                   <  2:4243583/149‑1 (MQ=255)
                                                                                                    atatTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAATATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGAtg                                                  >  1:1821557/1‑149 (MQ=255)
                                                                                                    atatTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGAtg                                                  >  2:870285/1‑149 (MQ=255)
                                                                                                       tttGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGc                                               <  1:3982396/149‑1 (MQ=255)
                                                                                                         tGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCgg                                             <  1:1603203/149‑1 (MQ=255)
                                                                                                          gCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGa                                            <  1:3905995/149‑1 (MQ=255)
                                                                                                            aCCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGAtt                                          <  2:3025957/149‑1 (MQ=255)
                                                                                                               agcagAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTAc                                       >  2:150758/1‑149 (MQ=255)
                                                                                                                gcagAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTAcc                                      >  1:1821434/1‑149 (MQ=255)
                                                                                                                    aaaGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACt                                  >  1:3139570/1‑149 (MQ=255)
                                                                                                                    aaaGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACAAACt                                  >  2:1831420/1‑149 (MQ=255)
                                                                                                                    aaaGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATGTCCCAACt                                  >  2:98879/1‑149 (MQ=255)
                                                                                                                        tgatgaAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACgg                              <  1:1231278/149‑1 (MQ=255)
                                                                                                                          atgaAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAt                            <  1:492762/149‑1 (MQ=255)
                                                                                                                          atgaAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAt                            <  2:697743/149‑1 (MQ=255)
                                                                                                                          atgaAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAt                            <  2:576239/149‑1 (MQ=255)
                                                                                                                          atgaAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAt                            >  2:2224187/1‑149 (MQ=255)
                                                                                                                           tgaAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAtt                           <  2:892149/149‑1 (MQ=255)
                                                                                                                           tgaAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAtt                           >  2:1667997/1‑149 (MQ=255)
                                                                                                                            gaAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAttt                          <  1:132353/149‑1 (MQ=255)
                                                                                                                             aaaTTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTg                         >  2:2648982/1‑149 (MQ=255)
                                                                                                                                ttttGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTc                      <  1:4750503/149‑1 (MQ=255)
                                                                                                                                  ttGTTTCGATAAAGAGAATATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTccc                    <  1:816446/149‑1 (MQ=255)
                                                                                                                                  ttGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTccc                    <  2:1296090/149‑1 (MQ=255)
                                                                                                                                      ttCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCtgat                >  2:4470153/1‑149 (MQ=255)
                                                                                                                                         gATAAAGAGAATATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGAc             <  1:1118541/149‑1 (MQ=255)
                                                                                                                                          aTAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACa            <  1:2619208/149‑1 (MQ=255)
                                                                                                                                          aTAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACa            >  2:204731/1‑149 (MQ=255)
                                                                                                                                          aTAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACa            >  1:159682/1‑149 (MQ=255)
                                                                                                                                            aaaGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACAtt          <  2:2515810/149‑1 (MQ=255)
                                                                                                                                            aaaGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACAtt          <  2:1475131/149‑1 (MQ=255)
                                                                                                                                             aaGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTg         >  2:3168941/1‑149 (MQ=255)
                                                                                                                                                agaagaTTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCaa      >  2:1279208/1‑149 (MQ=255)
                                                                                                                                                agaagaTTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCaa      <  1:3628763/149‑1 (MQ=255)
                                                                                                                                                   agaTTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATaa   >  2:3685243/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAAC  >  CP009273/4570749‑4571045

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: