Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,315 G→A 56.2% Y79Y (TAC→TAT USA300HOU_0026 ← IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,3150GA56.2% 231.2 / inf 1187Y79Y (TAC→TATUSA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base G (244/264);  new base A (298/361);  total (556/631)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

Only 100 of 1200 total aligned reads displayed.

TTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTC  >  CP000730/36175‑36455
                                                                                                                                            |                                                                                                                                            
ttACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATg                                                                                                                                              <  2:1101521/141‑1 (MQ=32)
 tACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAATTAGCTCCATTTTCCTTTTATTTTGATGt                                                                                                                                             >  2:736489/1‑141 (MQ=17)
  aCCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTa                                                                                                                                            <  1:309796/141‑1 (MQ=32)
      aaTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGtc                                                                                                                                        <  2:62015/141‑1 (MQ=12)
        ttGTTTAATTAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTTCGCAACCAATTATCTAATGTATGTTCCTCTCCATCACTGGCACGATATAAATAGCTCCATTTTCCTTTTATTCTGATGTACGtctc                                                                                                                                      <  2:467487/141‑1 (MQ=11)
          gTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCtcat                                                                                                                                    >  1:835750/1‑141 (MQ=17)
             tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCaa                                                                                                                                 >  1:622761/1‑141 (MQ=17)
             tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGGTGTACGTCTCATCaa                                                                                                                                 >  2:188214/1‑141 (MQ=17)
             tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCaa                                                                                                                                 >  1:394981/1‑141 (MQ=32)
             tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGGCtcctcat                                                                                                                                 >  1:37822/1‑140 (MQ=17)
              actGAGACGTTTGATACACGCATCTGCTGAATGATTATCTCGTTGATTACGCAACCAAATATATAATGTATGCCCCTCTGTATCAATGGCACGATATAAATAGCACCATTTTCCTTTTATTTTGATATTCGTCTCATCAAt                                                                                                                                >  2:316178/3‑141 (MQ=11)
              aaTGAGACGTTTGATAAACGCATATGCTGAATGATCATCCCGTTGCTGACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAAt                                                                                                                                <  1:852298/141‑1 (MQ=12)
               aTGTGACGTTTGATAAACGCATATGTTGAATGATTCTCTCGTTGCTTACGCAACCAAATATCTACTGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATa                                                                                                                               <  2:930022/141‑1 (MQ=11)
               aTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATa                                                                                                                               <  1:177143/141‑1 (MQ=12)
                 gagaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACg                                                                                                                             <  1:949352/141‑1 (MQ=12)
                 gagaCGATTGAAAAACGCATATACTAAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACg                                                                                                                             <  1:963197/141‑1 (MQ=11)
                  ggactgttGAGAAGCGCATATGCTGAATGATTATCTGTTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGGATCAATGGCACGATATAAAGAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGc                                                                                                                            <  2:602737/135‑1 (MQ=11)
                  agaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGc                                                                                                                            <  1:892819/141‑1 (MQ=32)
                       tttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCAttt                                                                                                                       >  2:164307/1‑141 (MQ=32)
                        ttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTg                                                                                                                      >  2:635207/1‑141 (MQ=17)
                        ttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTAATTTCGATCTAATGTTTTTCAATACGCCACTTg                                                                                                                      >  2:225405/1‑141 (MQ=11)
                         tGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGt                                                                                                                     >  2:1089376/1‑141 (MQ=17)
                           ggacagcgctatGCGGAATGATTCGCTCGTTGTCTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGtaa                                                                                                                   <  2:959655/132‑1 (MQ=11)
                           aTAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTACATTTTCCTTTTATTTTGATATACGTCTAATCAATACGCCATTTGtaa                                                                                                                   >  2:578137/1‑141 (MQ=255)
                              aaCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGtaataa                                                                                                                <  1:158538/141‑1 (MQ=18)
                               aCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGTTTTACGTCTCATCACTACGCCATTTGTAATAAg                                                                                                               >  1:90386/1‑141 (MQ=255)
                                cgtatttGCTGAATGATTATGTTGTTGCTTACGCAACAAAATATCTAAGGTATGTCCCTCAGTATCGCGTGCTGGAGATTTATCGCTCGATTTGCCTTTTATTTTAATATACGTCTCATCAATACGCCATTTGCAATAAGc                                                                                                              <  1:1091942/135‑1 (MQ=2)
                                   atatGCTGAATGATTATCTCGTTGCTTACGCAACCATATATCTAATGCATGTCCCTCTGTATCAATGGCACGAGATTAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAACACGCCATTTGTAAGAAGCttt                                                                                                           <  1:695236/141‑1 (MQ=11)
                                    tatGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCtttt                                                                                                          <  1:352099/141‑1 (MQ=21)
                                    tatGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTTTTTTGTTTTACGTCTACTCAATCTGCCATTTGTACGACGCtttt                                                                                                          >  1:135331/1‑141 (MQ=11)
                                      tGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATTCGTCTCATCAATACGCCATTTGTATTTAACtttttt                                                                                                        >  1:243865/1‑141 (MQ=11)
                                      tGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCTTCAATGGCTCGTTATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGCCTCATCAATACGCCATTTGTAATAAGCtttttt                                                                                                        >  1:78107/1‑141 (MQ=255)
                                       gCTGAATGATTATCTCTTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCTTTTGTAATAAGCttttttt                                                                                                       >  1:551681/1‑140 (MQ=255)
                                       gCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCAATTTGTAATAAGCTTTTTTa                                                                                                       <  1:669237/141‑1 (MQ=38)
                                        cTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGATCCATTTTCCTTTTATTTTGTTTTTCGTCTTCTCAATACGTCTTTTGTTATAAGCtttttttt                                                                                                      >  1:743331/1‑139 (MQ=11)
                                        cTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTAt                                                                                                      <  2:977924/141‑1 (MQ=38)
                                          gcaTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATTCATGGCCCGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                                                                    <  2:1060485/139‑1 (MQ=14)
                                          gAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                                                                    >  2:782287/1‑141 (MQ=255)
                                           aaTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTTTGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCt                                                                                                   <  2:883362/141‑1 (MQ=255)
                                            aTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCtt                                                                                                  >  2:1078894/1‑141 (MQ=40)
                                              gATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCtttt                                                                                                >  2:746228/1‑141 (MQ=255)
                                               aTTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCttttt                                                                                               >  2:892978/1‑141 (MQ=40)
                                                ttATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGGTATACGTCTCTTAAATATACCATTTGTAGTTAGATTTTTTATGGtttttt                                                                                              >  1:282703/1‑140 (MQ=11)
                                                ttATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTTATCAATAAGCCATTTGTAATAATCTTTTTTATGCTTTttc                                                                                              >  1:1029352/1‑141 (MQ=14)
                                                ttATCGCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttc                                                                                              <  1:719676/141‑1 (MQ=255)
                                                  ttctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCACTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCGTTTGTAATAAGCTTTTTTATGCTTTctctt                                                                                            <  2:671721/140‑1 (MQ=14)
                                                  aTCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttctt                                                                                            <  1:653798/141‑1 (MQ=40)
                                                   tctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATGTGTAATAAGCTTTTTTATGCTTTttcttc                                                                                           >  2:864915/1‑141 (MQ=255)
                                                    ctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTcc                                                                                          <  1:871182/141‑1 (MQ=255)
                                                         tGCGTAAGCAACGAGATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTCATAAGCTTTTTTATGCTTTTTCTTCCAAAtt                                                                                     >  2:838133/1‑141 (MQ=14)
                                                           cgtaagCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAAtttt                                                                                   >  2:717009/6‑140 (MQ=255)
                                                           cTTACGCAACCAAATATATAATGTATGTCACTCTGCATCAATGCCACGATATAAATAGCACCATTTTCCTTTTATTTTGATGTACGTCTGATCAATACGCCATTTTTAATAAGCTTTTTTATGCTTTTTCTTTCAAttttt                                                                                   >  2:1002335/1‑140 (MQ=14)
                                                              aCGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATa                                                                                >  2:590383/1‑141 (MQ=255)
                                                                  aaCCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCGAAATTTGATATaaa                                                                            <  1:575817/141‑1 (MQ=11)
                                                                   aCCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATCCaaaa                                                                           >  2:39121/1‑141 (MQ=255)
                                                                     cAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCCTTTGTAATAAGCTTTTTTATGTCTTTTCTTCAATTTTTGATACaacctt                                                                         >  2:1031824/1‑137 (MQ=14)
                                                                       aaTTTCCAATGTATGTTCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgg                                                                       <  2:761138/141‑1 (MQ=11)
                                                                        atatCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTTATGTACGTCTCTTCATTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTggg                                                                      >  1:422972/1‑141 (MQ=255)
                                                                         tatCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgggg                                                                     <  1:200023/141‑1 (MQ=255)
                                                                         tatCTAATCTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATCTACGTCTCATCAATACGCCATTTGTAATAAGCTTGTTTATGCTTTTTCTTCCAAATTTTATATACAAttagcg                                                                     >  2:210165/1‑137 (MQ=11)
                                                                           tCTAATGTATGTCCCTCTGTATCACTGGCACGATATCAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTACCAATTTTAATATaaccttgagcta                                                                   >  2:428258/1‑131 (MQ=11)
                                                                           tCTAATGTATCTCCATCTGCATCACTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCa                                                                   <  1:112689/141‑1 (MQ=255)
                                                                            cTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTTGGGCat                                                                  >  2:441653/1‑141 (MQ=255)
                                                                                tgtatgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTc                                                              <  1:220196/141‑1 (MQ=255)
                                                                                   atgtCCCTCTGTATCAAGGGCACGATATAAATAGCTCCCTTTTCCTTTTATTTTGATATACGGCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATAGTCTTg                                                           >  2:251546/1‑141 (MQ=11)
                                                                                        ccTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTTAAccc                                                      >  1:993370/1‑141 (MQ=21)
                                                                                          tctGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATCTAAAATTTGGGCATATTTTTGATACCaa                                                    >  1:480575/1‑141 (MQ=11)
                                                                                           ctGTATCAATGGCACGATATAAATCGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATCCTTTTTCTTTCAAATTTGATATAAAATTGGGGCCTGTTCCTgggatccga                                                   >  1:771372/1‑133 (MQ=11)
                                                                                            tGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAAGAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACg                                                  >  2:799415/1‑141 (MQ=21)
                                                                                            ggtgtCAATGGCACGATATAAAGAGCTCCATTTTCCTTTTATTTTGATAGACGGCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTGTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACg                                                  <  2:555541/137‑1 (MQ=11)
                                                                                             gtaTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACgg                                                 <  1:922816/139‑1 (MQ=255)
                                                                                                tCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACTTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTTTTCCAATTTTGATATAAAATTGGGGCATTTTCTTTAACCCAAAGGTAg                                              >  2:70008/1‑141 (MQ=11)
                                                                                                       cACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTg                                       <  2:537936/141‑1 (MQ=255)
                                                                                                          gATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAat                                    >  1:256977/1‑141 (MQ=255)
                                                                                                           atatAAATAGCTTCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCGTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatg                                   <  1:824597/141‑1 (MQ=255)
                                                                                                           atatAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatg                                   <  2:293196/141‑1 (MQ=18)
                                                                                                           atatAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatg                                   <  2:825158/141‑1 (MQ=21)
                                                                                                           agataAATAGCTCCATTTTCCTTTTATGTTGATGTACGTCTCATCAATACGCCATGTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatg                                   <  2:150735/139‑1 (MQ=255)
                                                                                                            tataAATAGCTGCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatga                                  <  1:459574/141‑1 (MQ=21)
                                                                                                               aaaTAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAAtgatga                               >  1:501211/1‑141 (MQ=255)
                                                                                                                aaTAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTTTTTTTTTTTTTATAAAAAATTGGGGCAAATTTTTGAATCCAGCGGTGGTCCGTTATAtgaggat                              >  2:514763/1‑140 (MQ=11)
                                                                                                                   atcTCCCTTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGt                           <  2:333412/139‑1 (MQ=11)
                                                                                                                   aGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGCATGTTGAACGt                           >  2:478778/1‑141 (MQ=11)
                                                                                                                      tcctttttccttttATTTTTATATACGTCTCATCAATACTCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTTATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTa                        <  1:1055908/141‑1 (MQ=11)
                                                                                                                       cctttttCCTTTTATTTTGATATACGCCTCGTCAATACGACATTTGTAATCAGGTCTTTTCTGCTTTTTCTTCCAAATTTGATAGAAAATTGGGGCATAGTCTTGAACCCGACGGTAGATCGTTGAATGATGAACGTTTac                       <  2:369748/138‑1 (MQ=11)
                                                                                                                        cATTTTCCTTTTATTTCGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGACCGTTTaca                      <  1:338358/141‑1 (MQ=255)
                                                                                                                         aTTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTacac                     <  1:1079110/141‑1 (MQ=255)
                                                                                                                          ttttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAGTGATGAACGTTTAcacc                    <  1:443985/141‑1 (MQ=16)
                                                                                                                            ttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcaccac                  >  1:933024/1‑141 (MQ=40)
                                                                                                                             tCCTTGTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACg                 <  1:408131/141‑1 (MQ=255)
                                                                                                                              ccTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGGACGTTTACACCACGt                <  2:124817/141‑1 (MQ=18)
                                                                                                                              ccTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGt                <  2:948602/141‑1 (MQ=21)
                                                                                                                               cTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAGTTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGtt               >  1:526894/1‑141 (MQ=25)
                                                                                                                               cTTTTATTTTGATAAACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGtt               >  1:591537/1‑141 (MQ=255)
                                                                                                                                  ttATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTccc            <  2:401856/141‑1 (MQ=21)
                                                                                                                                  ttATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTccc            >  1:1000234/1‑141 (MQ=21)
                                                                                                                                    aTTTTGATGTACGTCTCATCAATACGCCATGTGTAATAAGCTTTCTTCAACATGGGCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt          <  1:799415/141‑1 (MQ=11)
                                                                                                                                    aTTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt          <  1:150291/141‑1 (MQ=255)
                                                                                                                                       ttGATGTACTTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAat       <  2:487250/141‑1 (MQ=255)
                                                                                                                                           tGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATAttt   >  1:600797/1‑141 (MQ=21)
                                                                                                                                            |                                                                                                                                            
TTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTC  >  CP000730/36175‑36455

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: