Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,504,888 T→G 6.2% intergenic (+23/+17) insQ → / ← ydcO putative insertion element transposase InsQ/putative transport protein YdcO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,504,8880TG 6.2% 301.9 / 10.0 97intergenic (+23/+17)insQ/ydcOputative insertion element transposase InsQ/putative transport protein YdcO
Reads supporting (aligned to +/- strand):  ref base T (34/57);  new base G (2/4);  total (36/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

CGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCG  >  NC_000913/1504739‑1505035
                                                                                                                                                     |                                                                                                                                                   
cGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCtt                                                                                                                                                     >  1:974697/1‑150 (MQ=255)
 gCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTc                                                                                                                                                    <  1:829708/150‑1 (MQ=255)
  ccGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCa                                                                                                                                                   >  1:1527515/1‑150 (MQ=255)
    gATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGGCTCGCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAgg                                                                                                                                                 <  1:743180/150‑1 (MQ=255)
        tAAATGGCGCTCGCAACATTTTAGCGTCGGGGCACGCCGTTCTTGACTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGCATCCAAGTCCTTCAGGGCGt                                                                                                                                             >  2:192114/1‑150 (MQ=255)
        tAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCACGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGt                                                                                                                                             >  2:186146/1‑150 (MQ=255)
         aaaTGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTg                                                                                                                                            >  1:680359/1‑150 (MQ=255)
         aaaTGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGCCCTGAAGGACGTg                                                                                                                                            <  2:702709/150‑1 (MQ=255)
           aTGGCGCTCGCAACATTTTAGCGGCGGGGCACGACGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTgga                                                                                                                                          <  2:1534626/150‑1 (MQ=255)
           aTGCCGCTCGCTACATTTTCGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGCTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGCCCTGAAGGACGTgga                                                                                                                                          <  2:1099893/150‑1 (MQ=255)
            tGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGCCCTGAAGGACGTggat                                                                                                                                         <  1:1140434/150‑2 (MQ=255)
              gcgcTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGCCCTGAAGGACGTggatcc                                                                                                                                       <  2:305716/150‑4 (MQ=255)
               cgcTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGAt                                                                                                                                      >  1:1323876/1‑150 (MQ=255)
                  tCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTc                                                                                                                                   <  1:689780/150‑1 (MQ=255)
                  tCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTc                                                                                                                                   <  1:481602/150‑1 (MQ=255)
                   cGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATTTCa                                                                                                                                  <  2:406522/150‑1 (MQ=255)
                   cGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCa                                                                                                                                  <  1:1079179/150‑1 (MQ=255)
                   cGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCa                                                                                                                                  <  2:1271428/150‑1 (MQ=255)
                   cGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCa                                                                                                                                  <  1:208975/150‑1 (MQ=255)
                    gCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCaa                                                                                                                                 >  2:725620/1‑150 (MQ=255)
                    gCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCaa                                                                                                                                 <  2:413578/150‑1 (MQ=255)
                     cAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAat                                                                                                                                >  2:223562/1‑150 (MQ=255)
                     cAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAat                                                                                                                                <  1:1537412/150‑1 (MQ=255)
                        cATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATc                                                                                                                             >  2:657276/1‑150 (MQ=255)
                         aTTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCg                                                                                                                            >  2:1562380/1‑150 (MQ=255)
                          ttttAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCgg                                                                                                                           <  1:199581/150‑1 (MQ=255)
                           tttAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGt                                                                                                                          <  1:792336/150‑1 (MQ=255)
                           tttAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGt                                                                                                                          <  1:272797/150‑1 (MQ=255)
                           tttAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGt                                                                                                                          <  2:638569/150‑1 (MQ=255)
                            ttAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGtt                                                                                                                         >  2:927604/1‑150 (MQ=255)
                            ttAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGtt                                                                                                                         <  2:1527515/150‑1 (MQ=255)
                              aGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTc                                                                                                                       <  2:325168/150‑1 (MQ=255)
                              aGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTc                                                                                                                       <  2:596486/150‑1 (MQ=255)
                                cggcggGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCTACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTg                                                                                                                     <  1:260589/150‑1 (MQ=255)
                                cggcagGGCACGCCGTTATAGCCTTTGAAGAGAAGGTGCACACAGGCCGCCCGTTGAAGCAGGAACCCACAGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTg                                                                                                                     >  2:861027/1‑150 (MQ=255)
                                  gcggGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTc                                                                                                                   <  1:1158927/150‑1 (MQ=255)
                                   cggGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGTAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCa                                                                                                                  <  1:1267037/150‑1 (MQ=255)
                                          cGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTa                                                                                                           >  1:1025363/1‑150 (MQ=255)
                                               ttCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTc                                                                                                      <  1:109571/150‑1 (MQ=255)
                                               ttCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTc                                                                                                      <  2:702353/150‑1 (MQ=255)
                                                  ttGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAac                                                                                                   >  2:1474390/1‑150 (MQ=255)
                                                  ttGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCCTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAac                                                                                                   <  2:1351001/150‑1 (MQ=255)
                                                   tGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAaca                                                                                                  >  1:126361/1‑150 (MQ=255)
                                                       tgtgGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAAcaccac                                                                                              <  2:1531471/150‑1 (MQ=255)
                                                        gtgGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACg                                                                                             <  1:1284805/150‑1 (MQ=255)
                                                         tgGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGt                                                                                            >  1:1312012/1‑150 (MQ=255)
                                                             gagaTGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACa                                                                                        <  1:499481/150‑1 (MQ=255)
                                                                 tGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACg                                                                                    <  1:992287/150‑1 (MQ=255)
                                                                  ggTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGc                                                                                   >  1:660451/1‑150 (MQ=255)
                                                                     gCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTc                                                                                >  1:322832/1‑150 (MQ=255)
                                                                        gTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCTGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCgg                                                                             >  1:474739/1‑150 (MQ=255)
                                                                        gTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCgg                                                                             <  2:126361/150‑1 (MQ=255)
                                                                         tCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCggt                                                                            <  1:642302/150‑2 (MQ=255)
                                                                            ggCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAAt                                                                         >  1:277933/1‑150 (MQ=255)
                                                                               cgccCGTTGAAGCAGGAACCCAACGAAATGATTCAGGCGAAAGCCAGAACGTAGCAGTGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGATTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTaa                                                                      >  2:506558/1‑150 (MQ=255)
                                                                                gccCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAAt                                                                     >  2:1501523/1‑150 (MQ=255)
                                                                                 cccGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATc                                                                    >  2:370859/1‑150 (MQ=255)
                                                                                 cccGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATc                                                                    <  1:202463/150‑1 (MQ=255)
                                                                                  ccGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATcc                                                                   <  1:1129014/150‑1 (MQ=255)
                                                                                  ccGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATcc                                                                   <  2:1301160/150‑1 (MQ=255)
                                                                                    gTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATcccc                                                                 >  2:1443448/1‑150 (MQ=255)
                                                                                     ttGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCa                                                                <  1:949201/150‑1 (MQ=255)
                                                                                      tGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCaa                                                               >  2:570588/1‑150 (MQ=255)
                                                                                          gCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACg                                                           <  2:826740/150‑1 (MQ=255)
                                                                                            aGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCa                                                         >  2:85192/1‑150 (MQ=255)
                                                                                              gAACCCTCGGAAATGATTGAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGAGTGGAGGATGTCAATATCGGTTTCTGTCAGAGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGa                                                       <  2:722077/150‑1 (MQ=255)
                                                                                                aCCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAAc                                                     <  1:927604/150‑1 (MQ=255)
                                                                                                      cGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATcc                                               <  1:1205337/150‑1 (MQ=255)
                                                                                                       gAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATccc                                              >  2:443687/1‑150 (MQ=255)
                                                                                                         aaTGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCgc                                            >  2:1617831/1‑149 (MQ=255)
                                                                                                           tGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGAcc                                          <  1:191530/150‑1 (MQ=255)
                                                                                                           tGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGAcc                                          <  2:909864/150‑1 (MQ=255)
                                                                                                               tCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCg                                      <  2:1025363/150‑1 (MQ=255)
                                                                                                                 aGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTc                                    <  1:1562380/150‑1 (MQ=255)
                                                                                                                  ggCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAACCGCAGCACCAATCCCGACCAGCGTCa                                   >  2:1343148/1‑150 (MQ=255)
                                                                                                                    ccaCAGCCTGAACGAAGCAGGGATCCACGTCCTTCACGGCGTGGAGGATGTCAATATAGGTTTCTGTCAGCGATTAAAATCAACACCACGTAACAAACGCCTCCCGCAATTAATACCCAAAACGCAGAACCAATCCCGACCAGCGTCAAt                                 >  2:1010191/3‑150 (MQ=255)
                                                                                                                      aCAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTGTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATcc                               <  2:1106243/150‑1 (MQ=255)
                                                                                                                         gCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACt                            <  2:277933/150‑1 (MQ=255)
                                                                                                                         gCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACt                            <  1:85192/150‑1 (MQ=255)
                                                                                                                            tGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGc                         <  1:186146/150‑1 (MQ=255)
                                                                                                                            tGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGc                         <  1:192114/150‑1 (MQ=255)
                                                                                                                              aaCGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCg                       <  2:62276/150‑1 (MQ=255)
                                                                                                                                 gTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTa                    <  1:570588/150‑1 (MQ=255)
                                                                                                                                   aGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTAcc                  <  2:322832/150‑1 (MQ=255)
                                                                                                                                   aGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCACCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTAcc                  >  2:11474/1‑150 (MQ=255)
                                                                                                                                     cAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTa                    >  1:1125546/1‑146 (MQ=255)
                                                                                                                                     cAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTa                    <  2:1125546/146‑1 (MQ=255)
                                                                                                                                       atcctCCACGCCCTGAAGGACGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGaa              >  2:1536505/5‑150 (MQ=255)
                                                                                                                                         cctCCACGCCCTGAAGGACGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGCCCAGCGTCAATCCACTTGCCGTTACCAGAAAt            >  1:695087/3‑150 (MQ=255)
                                                                                                                                          aTCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATg           <  2:1323876/150‑1 (MQ=255)
                                                                                                                                            ccACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTATTCCCCAAAACGGAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGcc         <  2:974697/150‑1 (MQ=255)
                                                                                                                                            ccACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGcc         <  2:660451/150‑1 (MQ=255)
                                                                                                                                              aCGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATTccac       <  1:11474/150‑1 (MQ=255)
                                                                                                                                              aCGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGccac       <  2:1312012/150‑1 (MQ=255)
                                                                                                                                               cGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAATACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGccacc      >  1:1290337/1‑150 (MQ=255)
                                                                                                                                                gTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGccacca     <  1:1443448/150‑1 (MQ=255)
                                                                                                                                                 tCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCcaccac    >  1:550681/1‑150 (MQ=255)
                                                                                                                                                   cTTCAGGGAGTTGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTAAACACAACGTAACAAACCCCTCCTCCAATTAAACACCAAAACGCAGACCCAATCACGAACACCGTCAATCCACTTGACGTTACCAGAAATGCCaccaccc  >  2:667710/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCG  >  NC_000913/1504739‑1505035

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