Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,468,165 A→G 100% intergenic (‑37/‑125) argI ← / → rraB ornithine carbamoyltransferase 1/protein inhibitor of RNase E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,468,1650AG99.3% 1132.7 / ‑3.3 293intergenic (‑37/‑125)argI/rraBornithine carbamoyltransferase 1/protein inhibitor of RNase E
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (165/126);  total (165/128)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 293 total aligned reads displayed.

ACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCAAACCCGGAACAACTG  >  CP009273/4468018‑4468310
                                                                                                                                                   |                                                                                                                                                 
aCCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGtt                                                                                                                                                 >  1:2066101/1‑150 (MQ=255)
 ccGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTc                                                                                                                                                <  2:1799833/150‑1 (MQ=255)
   gCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCaa                                                                                                                                              <  2:3550474/150‑1 (MQ=255)
   gCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCaa                                                                                                                                              >  1:4461973/1‑150 (MQ=255)
       ttCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTta                                                                                                                                          >  1:4720590/1‑150 (MQ=255)
        tCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTtat                                                                                                                                         <  1:3676904/150‑1 (MQ=255)
         cTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAAGTtata                                                                                                                                        >  2:766389/1‑150 (MQ=255)
           tATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATg                                                                                                                                      <  2:4604698/150‑1 (MQ=255)
           tATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATg                                                                                                                                      >  2:264739/1‑150 (MQ=255)
           tATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATg                                                                                                                                      <  2:2034619/150‑1 (MQ=255)
           tATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATg                                                                                                                                      >  1:2343634/1‑150 (MQ=255)
             tCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGga                                                                                                                                    >  2:2019598/1‑150 (MQ=255)
              cGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGgat                                                                                                                                   >  1:723128/1‑150 (MQ=255)
               ggCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGgatg                                                                                                                                  <  1:2018372/150‑1 (MQ=255)
               ggCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGgatg                                                                                                                                  >  2:2276891/1‑150 (MQ=255)
               ggCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGgatg                                                                                                                                  <  2:1870997/150‑1 (MQ=255)
                   ttCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGatta                                                                                                                              >  1:2208310/1‑150 (MQ=255)
                   ttCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGatta                                                                                                                              >  1:1970560/1‑150 (MQ=255)
                     cAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGattatt                                                                                                                            <  2:3292138/150‑1 (MQ=255)
                     cAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGattatt                                                                                                                            <  2:3772370/150‑1 (MQ=255)
                       gCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCa                                                                                                                          >  2:3642504/1‑150 (MQ=255)
                           cgcgGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTg                                                                                                                      >  1:2644004/1‑150 (MQ=255)
                            gcgGCTAACTGCAGCAGGCTGGTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCTATTTATATGGATGATTATTCATTTGc                                                                                                                     <  1:4118798/150‑1 (MQ=255)
                             cgGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCa                                                                                                                    >  1:3011514/1‑150 (MQ=255)
                             cgGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCa                                                                                                                    >  2:723015/1‑150 (MQ=255)
                              ggCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCaa                                                                                                                   <  2:3449201/150‑1 (MQ=255)
                              ggCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCaa                                                                                                                   >  1:158695/1‑150 (MQ=255)
                               gCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGATTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCaat                                                                                                                  <  1:2534617/150‑2 (MQ=255)
                                cTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGt                                                                                                                 >  2:1702367/1‑150 (MQ=255)
                                cTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGt                                                                                                                 >  1:1422941/1‑150 (MQ=255)
                                 tAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTc                                                                                                                <  1:3318354/150‑1 (MQ=255)
                                   aCTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCCTTTGCAAGTCTa                                                                                                              >  2:4876490/1‑150 (MQ=255)
                                     tGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTaaa                                                                                                            >  2:2875861/1‑150 (MQ=255)
                                         gcagGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCAt                                                                                                        >  1:1529355/1‑150 (MQ=255)
                                              cTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATc                                                                                                   >  2:952125/1‑150 (MQ=255)
                                                    aGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTc                                                                                             >  2:3353952/1‑150 (MQ=255)
                                                     gTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTca                                                                                            <  1:2727606/150‑1 (MQ=255)
                                                            tGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGt                                                                                     <  2:4125963/150‑1 (MQ=255)
                                                             ggCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTg                                                                                    >  2:4735062/1‑150 (MQ=255)
                                                               cGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGa                                                                                  >  1:4095298/1‑150 (MQ=255)
                                                               cGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGa                                                                                  <  2:2566383/150‑1 (MQ=255)
                                                                gTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAg                                                                                 >  1:2928340/1‑150 (MQ=255)
                                                                gTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAg                                                                                 <  1:1546385/150‑1 (MQ=255)
                                                                 tGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCATCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAgg                                                                                <  1:3526542/150‑1 (MQ=255)
                                                                    aaTCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                                                             >  1:2393795/1‑150 (MQ=255)
                                                                     aTCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCACGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAAt                                                                            >  2:3181704/1‑150 (MQ=255)
                                                                       cGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGt                                                                          >  2:1616401/1‑150 (MQ=255)
                                                                       cGAGAAATTTCAGGAAATGCTAATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATATATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAAAGt                                                                          >  2:397892/1‑150 (MQ=255)
                                                                          gTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGGGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAg                                                                       <  2:4312766/150‑1 (MQ=255)
                                                                           tAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgt                                                                      <  1:4289751/150‑1 (MQ=255)
                                                                            aaTTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtg                                                                     <  1:4994846/150‑1 (MQ=255)
                                                                              tttCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtggt                                                                   <  1:4886092/150‑1 (MQ=255)
                                                                              tttCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtggt                                                                   <  2:2373410/150‑1 (MQ=255)
                                                                                 cAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTTTCACAAAGGTGGAGGCAATGTCAGTGgtgtg                                                                >  2:1789370/1‑150 (MQ=255)
                                                                                  aGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGgtgtgt                                                               <  2:1129757/150‑1 (MQ=255)
                                                                                  aGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGgtgtgt                                                               <  2:2066846/150‑1 (MQ=255)
                                                                                    gAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGa                                                             <  1:1228729/150‑1 (MQ=255)
                                                                                     aaaTGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGAc                                                            >  2:2175857/1‑150 (MQ=255)
                                                                                      aaTGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGACTTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACa                                                           >  2:563989/1‑150 (MQ=255)
                                                                                        tGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAAt                                                         >  2:1441310/1‑150 (MQ=255)
                                                                                         gCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATa                                                        <  2:3199408/150‑1 (MQ=255)
                                                                                            tATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAaga                                                     >  1:4811529/1‑150 (MQ=255)
                                                                                             aTGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAagag                                                    >  2:4051971/1‑150 (MQ=255)
                                                                                              tGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAATGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGt                                                   >  1:3633465/1‑150 (MQ=255)
                                                                                              tGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGt                                                   >  1:1901099/1‑150 (MQ=255)
                                                                                               gATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTa                                                  >  1:4375413/1‑150 (MQ=255)
                                                                                                aTAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTAt                                                 >  1:424339/1‑150 (MQ=255)
                                                                                                 tAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATc                                                >  1:3756321/1‑150 (MQ=255)
                                                                                                 tAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATc                                                >  1:5050976/1‑150 (MQ=255)
                                                                                                 tAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATc                                                >  1:1037021/1‑150 (MQ=255)
                                                                                                  aaaaCCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCg                                               <  1:1071109/150‑1 (MQ=255)
                                                                                                  aaaaCCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCg                                               >  2:2329052/1‑150 (MQ=255)
                                                                                                  aaaaCCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCg                                               <  2:1020765/150‑1 (MQ=255)
                                                                                                    aaCCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGc                                             <  2:2643941/150‑1 (MQ=255)
                                                                                                      cccGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAg                                           >  2:4550204/1‑150 (MQ=255)
                                                                                                       ccGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAgg                                          <  2:4487936/150‑1 (MQ=255)
                                                                                                         ggACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGAc                                        <  1:5131136/150‑1 (MQ=255)
                                                                                                             aTAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTa                                    >  2:1305184/1‑150 (MQ=255)
                                                                                                              tAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTaa                                   <  2:8370/150‑1 (MQ=255)
                                                                                                              tAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTaa                                   >  1:3284165/1‑150 (MQ=255)
                                                                                                               aGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAag                                  >  2:4268923/1‑150 (MQ=255)
                                                                                                               aGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAAGTTATATGGATGATTCTTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAag                                  >  2:1959073/1‑150 (MQ=255)
                                                                                                                 aTCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAagag                                >  1:3800779/1‑150 (MQ=255)
                                                                                                                 aTCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAagag                                >  2:4666820/1‑150 (MQ=255)
                                                                                                                   cccTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGa                              <  1:1332844/150‑1 (MQ=255)
                                                                                                                     ctcctGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAAt                            >  1:4528838/1‑150 (MQ=255)
                                                                                                                         tgtgGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGc                        <  1:610632/150‑1 (MQ=255)
                                                                                                                          gtgGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTAGTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGcc                       <  2:4720590/150‑1 (MQ=255)
                                                                                                                             gCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATg                    <  1:2876833/150‑1 (MQ=255)
                                                                                                                               tAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGc                  >  2:1289480/1‑150 (MQ=255)
                                                                                                                                aaCGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCa                 >  1:4517262/1‑150 (MQ=255)
                                                                                                                                aaCGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCa                 <  1:1749988/150‑1 (MQ=255)
                                                                                                                                aaCGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCa                 >  2:3320224/1‑150 (MQ=255)
                                                                                                                                   gCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCAAAc              >  1:2587706/1‑150 (MQ=255)
                                                                                                                                    ccTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCAAAcc             >  2:369759/1‑150 (MQ=255)
                                                                                                                                      tCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCAAACCCg           <  2:3428456/150‑1 (MQ=255)
                                                                                                                                          tGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCAAACCCGGaac       >  1:846665/1‑150 (MQ=255)
                                                                                                                                           gAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCAAACCCGGaaca      >  2:1499666/1‑150 (MQ=255)
                                                                                                                                           gAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCAAACCCGGaaca      <  2:5032224/150‑1 (MQ=255)
                                                                                                                                               tAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCAAACCCGGAACAACTg  >  2:2994242/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCAAACCCGGAACAACTG  >  CP009273/4468018‑4468310

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: