Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 828,374 G→T 100% V11V (GTC→GTA ybiA ← DUF1768 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273828,3740GT95.3% 796.6 / 25.8 249V11V (GTC→GTAybiADUF1768 family protein
Reads supporting (aligned to +/- strand):  ref base G (7/5);  new base T (115/122);  total (122/127)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 249 total aligned reads displayed.

TTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTTTT  >  CP009273/828227‑828522
                                                                                                                                                   |                                                                                                                                                    
tttttCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTAc                                                                                                                                                    <  1:2649483/150‑1 (MQ=255)
   ttCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACgtg                                                                                                                                                 >  1:337669/1‑150 (MQ=255)
      atcaAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTg                                                                                                                                              <  2:35012/150‑1 (MQ=255)
       tcaAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTgg                                                                                                                                             >  1:4251062/1‑150 (MQ=255)
         aaGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGAt                                                                                                                                           >  1:306520/1‑150 (MQ=255)
         aaGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGAt                                                                                                                                           >  2:677099/1‑150 (MQ=255)
         aaGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGAt                                                                                                                                           >  1:1649410/1‑150 (MQ=255)
                 tttGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTtgtg                                                                                                                                   >  2:4223628/1‑150 (MQ=255)
                    gtgCCTGAAATTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTc                                                                                                                                <  1:3880152/150‑1 (MQ=255)
                        cTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAAc                                                                                                                            <  2:3204997/150‑1 (MQ=255)
                        cTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAAc                                                                                                                            >  1:1422159/1‑150 (MQ=255)
                        cTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAAc                                                                                                                            <  2:1788285/150‑1 (MQ=255)
                        cTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAAc                                                                                                                            >  1:2575525/1‑150 (MQ=255)
                         tGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACg                                                                                                                           <  2:4509039/150‑1 (MQ=255)
                          gAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACgg                                                                                                                          <  2:3940212/150‑1 (MQ=255)
                          gAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACgg                                                                                                                          >  1:3083254/1‑150 (MQ=255)
                            aaaTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGc                                                                                                                        >  1:4356209/1‑150 (MQ=255)
                              aTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCAt                                                                                                                      <  1:1204099/150‑1 (MQ=255)
                                aGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTa                                                                                                                    >  2:1935258/1‑150 (MQ=255)
                                         aGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCAtt                                                                                                           <  2:3468494/150‑1 (MQ=255)
                                           gTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTAt                                                                                                         >  2:392048/1‑150 (MQ=255)
                                            tAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATg                                                                                                        >  2:2802499/1‑150 (MQ=255)
                                             aGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGt                                                                                                       <  2:3341840/150‑1 (MQ=255)
                                                 cAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtggt                                                                                                   <  2:1315495/150‑1 (MQ=255)
                                                  aGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtggtg                                                                                                  >  1:2997608/1‑150 (MQ=255)
                                                   ggTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtggtgg                                                                                                 <  2:2170969/150‑1 (MQ=255)
                                                   ggTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtggtgg                                                                                                 <  2:794477/150‑1 (MQ=255)
                                                      tttCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTa                                                                                              <  2:4368437/150‑1 (MQ=255)
                                                        tCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTAtt                                                                                            <  2:4087670/150‑1 (MQ=255)
                                                         cccGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTg                                                                                           <  1:3536504/150‑1 (MQ=255)
                                                          ccGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGc                                                                                          <  2:4153435/150‑1 (MQ=255)
                                                          ccGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGc                                                                                          <  2:2358050/150‑1 (MQ=255)
                                                                        gggCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcag                                                                            >  2:3069399/1‑150 (MQ=255)
                                                                         ggCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcaga                                                                           >  1:3572872/1‑150 (MQ=255)
                                                                           ccAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaa                                                                         >  1:1055698/1‑150 (MQ=255)
                                                                           ccAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaa                                                                         <  2:919806/150‑1 (MQ=255)
                                                                             aGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaaca                                                                       <  1:2357064/150‑1 (MQ=255)
                                                                             aGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaaca                                                                       <  1:1223570/150‑1 (MQ=255)
                                                                             aGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaaca                                                                       <  2:3036142/150‑1 (MQ=255)
                                                                              ggcggcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaacag                                                                      <  2:631144/150‑1 (MQ=255)
                                                                               gcggcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaccaga                                                                     <  2:1259761/150‑1 (MQ=255)
                                                                               gcggcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGGAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaacaga                                                                     <  2:2976517/150‑1 (MQ=255)
                                                                               gaggcAAAGTTAGATTAATCACCGTAGTCATCGCTGGTGCTGTAGAATTTTATTATGGTGTCTTTCTTTACGTGTTGGATCCTTTGTGCTAGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaatcaga                                                                     <  2:1402570/148‑1 (MQ=255)
                                                                                cggcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaacagaa                                                                    <  2:2422363/150‑1 (MQ=255)
                                                                                cggcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGCGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaacagaa                                                                    >  2:3418010/1‑150 (MQ=255)
                                                                                 ggcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaacagaaa                                                                   >  2:2069429/1‑150 (MQ=255)
                                                                                   cAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGaaaaa                                                                 >  2:221976/1‑150 (MQ=255)
                                                                                    aaaGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAAc                                                                >  2:2411983/1‑150 (MQ=255)
                                                                                    aaaGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAAc                                                                <  1:1139760/150‑1 (MQ=255)
                                                                                         tAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTa                                                           >  2:2660432/1‑150 (MQ=255)
                                                                                          aGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTAc                                                          >  1:4473295/1‑150 (MQ=255)
                                                                                          aGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTAc                                                          >  2:2899941/1‑150 (MQ=255)
                                                                                           gAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTAcc                                                         >  2:4264422/1‑150 (MQ=255)
                                                                                              aaaTCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCt                                                      >  1:2052010/1‑150 (MQ=255)
                                                                                               aaTCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTg                                                     <  1:744276/150‑1 (MQ=255)
                                                                                                aTCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGa                                                    >  1:4136150/1‑150 (MQ=255)
                                                                                                aTCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGa                                                    >  1:3713110/1‑150 (MQ=255)
                                                                                                aTCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGa                                                    >  1:3662707/1‑150 (MQ=255)
                                                                                                 tCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGaa                                                   <  1:4058627/150‑1 (MQ=255)
                                                                                                  cACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTGTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGaaa                                                  <  1:2201574/150‑1 (MQ=255)
                                                                                                     cGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAAcc                                               >  2:535567/1‑150 (MQ=255)
                                                                                                     cGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAAcc                                               <  2:1173802/150‑1 (MQ=255)
                                                                                                     cGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAAcc                                               >  1:2577134/1‑150 (MQ=255)
                                                                                                     cGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAAcc                                               <  1:3544043/150‑1 (MQ=255)
                                                                                                       tAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGa                                             <  1:1067238/150‑1 (MQ=255)
                                                                                                       tAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGa                                             <  1:4428543/150‑1 (MQ=255)
                                                                                                         gTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGaaa                                           >  2:1040776/1‑150 (MQ=255)
                                                                                                          tCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGaaaa                                          <  2:1271256/150‑1 (MQ=255)
                                                                                                          tCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGaaaa                                          <  2:3511929/150‑1 (MQ=255)
                                                                                                           cATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGaaaaa                                         <  2:4593828/150‑1 (MQ=255)
                                                                                                            aTCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAAt                                        >  1:3438136/1‑150 (MQ=255)
                                                                                                              cGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGc                                      >  1:2984057/1‑150 (MQ=255)
                                                                                                                cTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCca                                    >  1:3489838/1‑150 (MQ=255)
                                                                                                                 tGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCcac                                   >  2:2869594/1‑150 (MQ=255)
                                                                                                                    tGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAat                                <  1:2836811/150‑1 (MQ=255)
                                                                                                                     gCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAata                               <  2:985549/150‑1 (MQ=255)
                                                                                                                         tAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTgg                           <  2:1316462/150‑1 (MQ=255)
                                                                                                                            aaaTTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTg                        >  1:2485982/1‑150 (MQ=255)
                                                                                                                             aaTTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATAATGTTTGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGt                       >  2:2507078/1‑150 (MQ=255)
                                                                                                                              aTTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGtt                      >  1:4493035/1‑150 (MQ=255)
                                                                                                                              aTTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGtt                      >  1:471006/1‑150 (MQ=255)
                                                                                                                               tttATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGttt                     <  2:2615103/150‑1 (MQ=255)
                                                                                                                                ttATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTta                    <  1:611262/150‑1 (MQ=255)
                                                                                                                                ttATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTta                    <  1:2739519/150‑1 (MQ=255)
                                                                                                                                ttATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCACGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTta                    >  2:1270349/1‑150 (MQ=255)
                                                                                                                                   tgatgGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATAc                 <  2:3604543/150‑1 (MQ=255)
                                                                                                                                     atgGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTCTACAg               >  2:1726601/1‑150 (MQ=255)
                                                                                                                                      tgGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGt              >  1:1991772/1‑150 (MQ=255)
                                                                                                                                      tgGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGt              <  2:1452840/150‑1 (MQ=255)
                                                                                                                                       ggTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCCCAATATTGGCTGTTTATACAGTa             >  2:4133633/1‑150 (MQ=255)
                                                                                                                                       ggTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTa             >  1:1788903/1‑150 (MQ=255)
                                                                                                                                        gtgtCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTAt            >  2:1101878/1‑150 (MQ=255)
                                                                                                                                        gtgtCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTAt            <  1:3930919/150‑1 (MQ=255)
                                                                                                                                        gtgtCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTAt            <  1:1606872/150‑1 (MQ=255)
                                                                                                                                         tgtCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTAtt           >  2:4233380/1‑150 (MQ=255)
                                                                                                                                         tgtCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTAtt           >  1:4215860/1‑150 (MQ=255)
                                                                                                                                           tCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTc         <  1:4308657/150‑1 (MQ=255)
                                                                                                                                              tGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAgg      <  1:1049098/150‑1 (MQ=255)
                                                                                                                                              tGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAgg      <  1:4223450/150‑1 (MQ=255)
                                                                                                                                                cATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGtt    >  1:1351909/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTTTT  >  CP009273/828227‑828522

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: