Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 733,570 A→G 5.5% E8G (GAA→GGA)  ybfD → H repeat‑associated putative transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273733,5700AG 5.5% 213.1 / 3.1 73E8G (GAA→GGA) ybfDH repeat‑associated putative transposase
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base G (2/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

CTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTA  >  CP009273/733299‑733818
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        
cgaTTAATACATTAACGAAGGTTCAGGTATTAAATTCCCGGTTAAGACGTAATATGCGGTAAGCGGTCTTGGATAGAAACTATCTCGAGTGAGTCCTTGCAGGGACCGCGCTTTCGTACTCTTTCCCTGCTTTTTGTACTAAATAAGTCGCATTTTATGTGTAAAAATTGAGATATTCCTAATTACCTGAAGCTGTTTTTTATTGCTTATACCTGATCAAATACGTCGTACATAATTAAGGAGAACAACATGGAACTTAAAAAATGGATCGac                                                                                                                                                                                                                                                         <  2:245445/271‑2 (MQ=18)
           ttAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAATTAGAACAAAATGGAACTTAAAAAATTGATAGAAAATTTTTttatta                                                                                                                                                                                                                                            >  2:573926/1‑275 (MQ=255)
                             ttCAAGGCAGAGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAGTTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTATTGCTTATACATGATCAAATACGTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACaa                                                                                                                                                                                                                           >  1:722802/1‑274 (MQ=255)
                                                  aaGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATCTTTCCCTGCTTTTTGTACTNAATAATTCGCATTTTATGTGTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATaaa                                                                                                                                                                                                      <  1:573939/274‑1 (MQ=255)
                                                       tcGAATAGCAGTCATGGTTTGAATCTATCTCGCGTCTGTTCTGGCGGGGATCTGTCTTTCGTAATCTTTCCCTGCTTTTTTNACTAAATANTTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGATAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATAACAGACAAGCCTGGAAAGTGGAACATAAATTATc                                                                                                                                                                                                 <  1:497671/273‑1 (MQ=17)
                                                        cGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAGTTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCgg                                                                                                                                                                                               <  2:722788/275‑1 (MQ=255)
                                                        cGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTGTCGTAATCTTTCCCTGCTTTTTGTACTACATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTAACCACGATTACAGACAACAATGGAAAGTGGAACAAAAATTATCgg                                                                                                                                                                                               >  2:674691/1‑275 (MQ=255)
                                                            aaGCAGCCATTGATATAAACTATCTCGCGTCAGTCCTTGCGGTGAGCGGGCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGANATTCCTTATTAGCTGAAGCTGTTTTTTATTGCTTATACATTATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGa                                                                                                                                                                                              <  1:674705/272‑1 (MQ=255)
                                                                 gcCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTCTTGCTTATACATGATCTACTACTTCTTACATAATTCAGTAGTACAAACTGGAACTTCACAACTTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAGGTGGACCATAAATTCTCGGCTATTCTAc                                                                                                                                                                                      >  2:678472/1‑275 (MQ=18)
                                                                       gATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTATTGCTTATACATGATACAATACTTCTTNCATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGa                                                                                                                                                                                >  1:679849/1‑275 (MQ=255)
                                                                                                        aGCGGGCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGc                                                                                                                                                <  1:210555/274‑1 (MQ=255)
                                                                                                           gggCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGa                                                                                                                                                                                                                                                          <  1:687634/165‑1 (MQ=255)
                                                                                                                  cGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAAAATAAATTATCGGATATTCTACTg                                                                                                                                                                                    >  1:432074/1‑228 (MQ=255)
                                                                                                                           tatttGCTTTGTGTACTAAATAATTCTTATTTTTTGTTTAAAAATTTAGTTTTTCCTTATTACCTGAATGTGTTTTTTATTGATTATACATGATCAAATTCTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGAttt                                                                                                                            <  2:619110/271‑1 (MQ=11)
                                                                                                                           tCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTCACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAACCTGGAAAGTGGAACATAAATTATCGTATAATCTACTGTTGACTATTTGTGCAGTTATTTCTGTTTCAGAAGGTGGGGAAGATATAGAGGAttt                                                                                                                            >  2:322798/1‑275 (MQ=255)
                                                                                                                                                                            atatTCCTTATTACCTGAAGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTAGTTGTGCCGTTATTTCTGGTGCAGa                                                                                                                                                    >  1:413778/1‑202 (MQ=21)
                                                                                                                                                                                tCCTTATTACCTGAAGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTa                                                                                                                                                                                                                                  >  2:536735/1‑120 (MQ=35)
                                                                                                                                                                                tCCTTATTACCTGAAGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGc                                                                                                                                                <  2:413771/202‑1 (MQ=35)
                                                                                                                                                                                    tattaCCTGAAGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGa                                                                                                                                                                                                                              <  1:536745/120‑1 (MQ=35)
                                                                                                                                                                                         ccTGAAGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTa                                                                                                                                                                                                                                  >  2:641371/1‑111 (MQ=35)
                                                                                                                                                                                            gAAGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATcccc                                                                                                                                                                                                                                       >  1:607615/1‑103 (MQ=35)
                                                                                                                                                                                             aaGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGa                                                                                                                                                                                                                              <  1:641385/111‑1 (MQ=35)
                                                                                                                                                                                                cTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGAtt                                                                                                                                                                                                                                   <  2:607602/103‑1 (MQ=35)
                                                                                                                                                                                                   ttttttATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCATAAGGCTGGGAAGa                                                                                                                                        >  1:401445/1‑191 (MQ=14)
                                                                                                                                                                                                       ttATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGAtata                                                                                                                                    <  2:401438/191‑1 (MQ=14)
                                                                                                                                                                                                        tATTGCTTATACATTATCAAATACTTCTCACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCa                                               <  1:322804/275‑1 (MQ=25)
                                                                                                                                                                                                                   cATGATCAAAGACTTCTGACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTANCGGATATTCTAGTGTTGACTATGTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATcc                                    <  1:678486/275‑1 (MQ=255)
                                                                                                                                                                                                                       aTCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGa                                                                                                                                                                                                             >  1:697698/1‑102 (MQ=40)
                                                                                                                                                                                                                       aTCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGa                                                                                                                                                                                              >  2:325552/1‑117 (MQ=39)
                                                                                                                                                                                                                       aTCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTTTa                                >  1:594441/1‑275 (MQ=25)
                                                                                                                                                                                                                           aaTACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACAt                                                                                                                                                                                                         <  2:697684/102‑1 (MQ=40)
                                                                                                                                                                                                                           aaTACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATAtt                                                                                                                                                                                          <  1:325558/117‑1 (MQ=39)
                                                                                                                                                                                                                                ttcttACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGaa                                                                                                                                                                                                                                                         <  2:188977/49‑1 (MQ=40)
                                                                                                                                                                                                                                ttcttACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAgg                                                                                                                                >  1:394645/1‑170 (MQ=37)
                                                                                                                                                                                                                                  cttACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTAt                                                                                                                                                                                                  >  1:14561/1‑102 (MQ=35)
                                                                                                                                                                                                                                    tACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTAc                                                                                                                                                                                      >  2:674651/1‑112 (MQ=35)
                                                                                                                                                                                                                                    tACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCACCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGAttt                                                                                                                            <  2:394638/170‑1 (MQ=25)
                                                                                                                                                                                                                                      cATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGa                                                                                                                                                                                              <  2:14561/102‑1 (MQ=35)
                                                                                                                                                                                                                                        tAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGtt                                                                                                                                                                                  <  1:674665/112‑1 (MQ=32)
                                                                                                                                                                                                                                         aaTTAAGGAGAACAAACTGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACt                                                                                                                                                                              >  1:47406/1‑115 (MQ=17)
                                                                                                                                                                                                                                         aaTTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTAc                                                                                                                                                                                      >  2:564498/1‑107 (MQ=17)
                                                                                                                                                                                                                                         aaTTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTAc                                                                                                                                                                                      >  2:640961/1‑107 (MQ=17)
                                                                                                                                                                                                                                          aTTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAAttttt             >  2:618895/1‑273 (MQ=25)
                                                                                                                                                                                                                                             aaGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGtt                                                                                                                                                                                  <  1:640975/107‑1 (MQ=17)
                                                                                                                                                                                                                                             aaGGAGAACAAAATGGAACTTAAAAAATTGATGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGtt                                                                                                                                                                                  <  1:564510/107‑1 (MQ=17)
                                                                                                                                                                                                                                             aaGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTAttt                                                                                                                                                                          <  2:47405/115‑1 (MQ=32)
                                                                                                                                                                                                                                              aGGATAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGNCAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTACTGTTCACGATACGATTGCCAGAGTTGTa                                       <  1:503906/245‑1 (MQ=18)
                                                                                                                                                                                                                                              aGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGCGAAACACATCCCGATTTTTTGAAGc                                                                                               >  1:563284/1‑189 (MQ=25)
                                                                                                                                                                                                                                               ggAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATgaga        <  1:618908/275‑1 (MQ=38)
                                                                                                                                                                                                                                                  gaACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGATTTTTTGAAGCAata                                                                                           <  2:563272/189‑1 (MQ=35)
                                                                                                                                                                                                                                                   accAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATCAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGCTGTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCATGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTGCTTGTAAATTTTCTgagatgtg    >  2:176629/3‑271 (MQ=11)
                                                                                                                                                                                                                                                   aaCAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAAt                                                                                                                                                                                                     >  2:366684/1‑82 (MQ=33)
                                                                                                                                                                                                                                                   aaCAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTAt                                                                                                                                                                                                  >  2:352184/1‑85 (MQ=33)
                                                                                                                                                                                                                                                   aaCAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTAt                                                                                                                                                                                                  >  1:398125/1‑85 (MQ=33)
                                                                                                                                                                                                                                                   aaCAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTAt                                                                                                                                                                                                  >  2:454613/1‑85 (MQ=33)
                                                                                                                                                                                                                                                   aaCAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCg                                                                                                                                                                                                >  2:157020/1‑87 (MQ=33)
                                                                                                                                                                                                                                                   aaCAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGatat                                                                                                                                                                                           >  1:189410/1‑92 (MQ=33)
                                                                                                                                                                                                                                                   aaCAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGc      >  1:71420/1‑273 (MQ=38)
                                                                                                                                                                                                                                                     cAAAATGGAACTTAAGAAATTGATAGAAGATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAAt                                                                                                                                                                                                     >  1:510184/1‑80 (MQ=11)
                                                                                                                                                                                                                                                     cAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTAt                                                                                                                                                                                                  >  1:188074/1‑83 (MQ=33)
                                                                                                                                                                                                                                                     cAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTtg                                                                                                                                                                         >  2:575465/1‑108 (MQ=32)
                                                                                                                                                                                                                                                     cAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGa                                                                                                                                                    >  2:289928/1‑129 (MQ=32)
                                                                                                                                                                                                                                                     cAAAATGGAACTTAAAAAATTGATGGAACATAGTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCACGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTACTGTTCACGATACCATTGCCAGAGTTGTATCCTTTATCTTTCCTGCGAAATTGCATTAAAACtttt  >  2:586019/1‑274 (MQ=11)
                                                                                                                                                                                                                                                       aaaTGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATc                                                                                                                                                                                                 <  1:366691/82‑1 (MQ=33)
                                                                                                                                                                                                                                                       aaaTGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGata                                                                                                                                                                                            <  1:157022/87‑1 (MQ=33)
                                                                                                                                                                                                                                                       aaaTGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGa                                                                                                                                                                                              <  1:454620/85‑1 (MQ=33)
                                                                                                                                                                                                                                                       aaaTGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGa                                                                                                                                                                                              <  1:352190/85‑1 (MQ=33)
                                                                                                                                                                                                                                                       aaaTGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGa                                                                                                                                                                                              <  2:398118/85‑1 (MQ=33)
                                                                                                                                                                                                                                                       aaaTGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTa                                                                                                                                                                                       <  2:189408/92‑1 (MQ=33)
                                                                                                                                                                                                                                                        aaTGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCag                                              >  1:40068/1‑228 (MQ=34)
                                                                                                                                                                                                                                                         aTGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATc                                                                                                                                                                                                 <  2:510175/80‑1 (MQ=34)
                                                                                                                                                                                                                                                         aTGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGa                                                                                                                                                                                              <  2:188072/83‑1 (MQ=33)
                                                                                                                                                                                                                                                         aTGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGcc                                                                                                                                                                     <  1:575478/108‑1 (MQ=32)
                                                                                                                                                                                                                                                         aTGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGc                                                                                                                                                <  1:289933/129‑1 (MQ=33)
                                                                                                                                                                                                                                                                   aaaaaTTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCAtt                                                   >  1:332729/1‑212 (MQ=34)
                                                                                                                                                                                                                                                                       aTTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCa                                               <  2:332723/212‑1 (MQ=34)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        
CTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTATTGCTTATACATGATCAAATACTTCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAATTGATGGAACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAGTGGAACATAAATTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCTGTATCTGTCCTGCGAAATTTCATGAGAGCTTTA  >  CP009273/733299‑733818

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