Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,468,165 A→G 90.6% 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,1650AG90.6% 843.6 / 54.0 256intergenic (‑37/‑125)argI/rraBornithine carbamoyltransferase 1/protein inhibitor of RNase E
Reads supporting (aligned to +/- strand):  ref base A (13/11);  new base G (118/114);  total (131/125)
Fisher's exact test for biased strand distribution p-value = 8.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 256 total aligned reads displayed.

ACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCAAACCCGGAACAACTG  >  CP009273/4468018‑4468310
                                                                                                                                                   |                                                                                                                                                 
aCCGCTTTTCTTATCGTCTTTCAGCTTCGCGGCTAACTGCATCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGtt                                                                                                                                                 <  1:4152167/150‑1 (MQ=255)
aCCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGACTTAAAGtt                                                                                                                                                 >  1:1074760/1‑150 (MQ=255)
aCCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGtt                                                                                                                                                 >  2:3831424/1‑150 (MQ=255)
aCCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGtt                                                                                                                                                 <  2:277320/150‑1 (MQ=255)
aCCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGtt                                                                                                                                                 >  2:2254665/1‑150 (MQ=255)
aCCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGtt                                                                                                                                                 <  1:3941833/150‑1 (MQ=255)
 ccGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTc                                                                                                                                                <  2:1267605/150‑1 (MQ=255)
    cTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAAt                                                                                                                                             >  1:1800641/1‑150 (MQ=255)
      tttCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAAttt                                                                                                                                           >  1:1040505/1‑150 (MQ=255)
          ttATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTatat                                                                                                                                       >  2:4284740/1‑150 (MQ=255)
          ttATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTatat                                                                                                                                       <  2:3586042/150‑1 (MQ=255)
           tATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATg                                                                                                                                      >  1:4518034/1‑150 (MQ=255)
            aTCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATgg                                                                                                                                     <  1:713140/150‑1 (MQ=255)
            aTCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATgg                                                                                                                                     <  2:2823673/150‑1 (MQ=255)
              cGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGgat                                                                                                                                   <  2:4506228/150‑1 (MQ=255)
              cGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGgat                                                                                                                                   >  2:2684559/1‑150 (MQ=255)
               ggCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGgatg                                                                                                                                  >  2:1688640/1‑150 (MQ=255)
                 cTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGgatgat                                                                                                                                >  1:276039/1‑150 (MQ=255)
                  tttCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGatt                                                                                                                               <  1:2787770/150‑1 (MQ=255)
                   ttCAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGatta                                                                                                                              <  1:2156513/150‑1 (MQ=255)
                    tCAGCTTCGCGGCTAACTGCAGCTGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGattat                                                                                                                             <  2:2351007/150‑1 (MQ=255)
                     cAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGattatt                                                                                                                            >  1:1129691/1‑150 (MQ=255)
                     cAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGattatt                                                                                                                            <  2:1314298/150‑1 (MQ=255)
                     cAGCTTCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGattatt                                                                                                                            <  1:2199445/150‑1 (MQ=255)
                         ttCGCGGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCAtt                                                                                                                        >  1:3062387/1‑150 (MQ=255)
                             cgGCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCa                                                                                                                    <  2:3317118/150‑1 (MQ=255)
                              ggCTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCaa                                                                                                                   >  2:1540157/1‑150 (MQ=255)
                                cTAACTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGt                                                                                                                 <  2:3799995/150‑1 (MQ=255)
                                   aCTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTa                                                                                                              <  2:2478122/150‑1 (MQ=255)
                                    cTGCAGCAGGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTaa                                                                                                             >  1:2262987/1‑150 (MQ=255)
                                        agcagGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCa                                                                                                         >  1:2842527/1‑150 (MQ=255)
                                        agcagGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCa                                                                                                         >  2:1443502/1‑150 (MQ=255)
                                           agGCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATaa                                                                                                      <  1:26342/150‑1 (MQ=255)
                                             gCTGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAAt                                                                                                    <  1:3623757/150‑1 (MQ=255)
                                               tGTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCt                                                                                                  <  2:580990/150‑1 (MQ=255)
                                                gTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCtt                                                                                                 <  1:1592796/150‑1 (MQ=255)
                                                gTTGAGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCtt                                                                                                 <  2:3723324/150‑1 (MQ=255)
                                                    aGTTCAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTc                                                                                             <  1:582075/150‑1 (MQ=255)
                                                        cAGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACaa                                                                                         >  2:2679625/1‑150 (MQ=255)
                                                         aGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACaaa                                                                                        >  2:465037/1‑150 (MQ=255)
                                                         aGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACaaa                                                                                        <  1:1553169/150‑1 (MQ=255)
                                                         aGCTGGCGTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACaaa                                                                                        >  2:411958/1‑150 (MQ=255)
                                                                gTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAg                                                                                 >  2:1654319/1‑150 (MQ=255)
                                                                gTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAg                                                                                 <  1:137105/150‑1 (MQ=255)
                                                                gTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAg                                                                                 >  1:3850283/1‑150 (MQ=255)
                                                                gTGAAATCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAg                                                                                 <  1:1921424/150‑1 (MQ=255)
                                                                   aaaTCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCa                                                                              >  2:2000701/1‑150 (MQ=255)
                                                                   aaaTCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCa                                                                              <  2:2087145/150‑1 (MQ=255)
                                                                     aTCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTATAGCATAAATCTTTGTCACAATGGTGGAGGCAAt                                                                            >  2:2389969/1‑150 (MQ=255)
                                                                     aTCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAAt                                                                            >  1:4484207/1‑150 (MQ=255)
                                                                     aTCGAGTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAAt                                                                            <  2:320506/150‑1 (MQ=255)
                                                                          gTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAg                                                                       <  1:4212383/150‑1 (MQ=255)
                                                                          gTAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAg                                                                       <  1:4212547/150‑1 (MQ=255)
                                                                           tAATTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgt                                                                      >  1:2251795/1‑150 (MQ=255)
                                                                            aaTTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtg                                                                     >  1:2815741/1‑150 (MQ=255)
                                                                            aaTTTCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtg                                                                     <  1:4713390/150‑1 (MQ=255)
                                                                               ttCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtggtg                                                                  <  1:4086695/150‑1 (MQ=255)
                                                                               ttCAGGAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtggtg                                                                  >  1:171475/1‑150 (MQ=255)
                                                                                   ggAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGtgtgtg                                                              >  1:2122930/1‑150 (MQ=255)
                                                                                   ggAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGtgtgtg                                                              <  1:1443502/150‑1 (MQ=255)
                                                                                   ggAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGtgtgtg                                                              >  1:3300068/1‑150 (MQ=255)
                                                                                    gAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGa                                                             >  1:3154456/1‑150 (MQ=255)
                                                                                    gAAATGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGa                                                             <  1:873669/150‑1 (MQ=255)
                                                                                      aaTGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACa                                                           >  1:3967537/1‑150 (MQ=255)
                                                                                       aTGCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACaa                                                          >  1:4128060/1‑150 (MQ=255)
                                                                                         gCTTATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATa                                                        <  1:237356/150‑1 (MQ=255)
                                                                                           ttATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAag                                                      <  2:3980255/150‑1 (MQ=255)
                                                                                           ttATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAag                                                      <  2:3160911/150‑1 (MQ=255)
                                                                                            tATGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAaga                                                     <  1:3477449/150‑1 (MQ=255)
                                                                                             aTGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAagag                                                    <  2:1800641/150‑1 (MQ=255)
                                                                                             aTGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAagag                                                    >  2:3659765/1‑150 (MQ=255)
                                                                                             aTGATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAagag                                                    >  1:2906814/1‑150 (MQ=255)
                                                                                               gATAAAACCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGtc                                                  >  2:1823740/1‑149 (MQ=255)
                                                                                                  aaaaCCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCg                                               <  1:753965/150‑1 (MQ=255)
                                                                                                    aaCCCGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGc                                             <  1:3831424/150‑1 (MQ=255)
                                                                                                       ccGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAgg                                          <  2:2907330/150‑1 (MQ=255)
                                                                                                        cGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa                                         >  1:1846532/1‑150 (MQ=255)
                                                                                                        cGGACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa                                         <  2:3955297/150‑1 (MQ=255)
                                                                                                          gACATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACa                                       >  2:347200/1‑150 (MQ=255)
                                                                                                            cATAGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACAtt                                     >  1:893007/1‑150 (MQ=255)
                                                                                                                   cccTCCTGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGa                              >  2:484483/1‑150 (MQ=255)
                                                                                                                    cctcctGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGaa                             <  1:4608432/150‑1 (MQ=255)
                                                                                                                        ctGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAg                         <  2:2208871/150‑1 (MQ=255)
                                                                                                                        ctGTGGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAg                         <  1:1961685/150‑1 (MQ=255)
                                                                                                                         tgtgGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGc                        >  1:1622281/1‑150 (MQ=255)
                                                                                                                         tgtgGCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGc                        >  1:1223800/1‑150 (MQ=255)
                                                                                                                            ggCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCAt                     >  1:3744780/1‑150 (MQ=255)
                                                                                                                             gCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATg                    >  1:3375458/1‑150 (MQ=255)
                                                                                                                             gCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATg                    >  1:3483400/1‑150 (MQ=255)
                                                                                                                             gCTAACGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATg                    <  1:148890/150‑1 (MQ=255)
                                                                                                                                 aCGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCaa                >  1:2636740/1‑150 (MQ=255)
                                                                                                                                 aCGCCTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCaa                <  2:462646/150‑1 (MQ=255)
                                                                                                                                  cGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCaaa               <  1:1613659/150‑1 (MQ=255)
                                                                                                                                    ccTCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCAAAcc             <  2:3061935/150‑1 (MQ=255)
                                                                                                                                      tCAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCAAACCCg           >  1:2696282/1‑150 (MQ=255)
                                                                                                                                       cAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCAAACCCgg          <  2:828819/150‑1 (MQ=255)
                                                                                                                                       cAATGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCAAACCCgg          >  1:3100712/1‑150 (MQ=255)
                                                                                                                                         aTGAATTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCAAACCCGGaa        >  1:1672347/1‑150 (MQ=255)
                                                                                                                                             aTTAAAGTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCAAACCCGGaacaac    >  2:1036943/1‑150 (MQ=255)
                                                                                                                                              ttAAAGTTAAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGCCATGGCAAACCCGGAACAACt   <  2:2610642/150‑1 (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: