Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,088,091 G→C 11.0% intergenic (+34/+16) yiiG → / ← frvR DUF3829 family lipoprotein/putative frv operon regulator; contains a PTS EIIA domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,088,0910GC11.0% 273.0 / 21.8 91intergenic (+34/+16)yiiG/frvRDUF3829 family lipoprotein/putative frv operon regulator; contains a PTS EIIA domain
Reads supporting (aligned to +/- strand):  ref base G (48/33);  new base C (4/6);  total (52/39)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

CTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCG  >  NC_000913/4087943‑4088239
                                                                                                                                                    |                                                                                                                                                    
cTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggc                                                                                                                                                     >  2:274419/1‑150 (MQ=255)
  gggATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGgcgc                                                                                                                                                   >  2:1688746/1‑150 (MQ=255)
      tAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATcc                                                                                                                                               <  1:1546563/150‑1 (MQ=255)
        aaGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCgc                                                                                                                                             >  1:158805/1‑149 (MQ=255)
         agagCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGac                                                                                                                                            <  1:1418054/150‑1 (MQ=255)
         agagCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGac                                                                                                                                            <  2:1130297/150‑1 (MQ=255)
         agagCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGac                                                                                                                                            <  1:23165/150‑1 (MQ=255)
          cagCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGaca                                                                                                                                           <  1:395581/149‑1 (MQ=255)
           cgCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGGTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGacac                                                                                                                                          <  1:1484030/149‑1 (MQ=255)
           agCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGacac                                                                                                                                          >  2:1374660/1‑150 (MQ=255)
             cAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGacactg                                                                                                                                        <  2:1604005/150‑3 (MQ=255)
               acatCAGGCCGCAAACTCACGCTGGATTGTCGAAGACTCTTTCACTATAGCATTACGCGAGTACAACGATATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTTATGTAAAATGCAGTGTCGGATGCGGCGCATCCGACACAGtt                                                                                                                                      <  2:143377/147‑1 (MQ=255)
                 ttCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCa                                                                                                                                    >  2:1370854/1‑150 (MQ=255)
                   cAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGTGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCTTCCGACACAGTTCAcc                                                                                                                                  >  2:250472/1‑150 (MQ=255)
                   cAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcc                                                                                                                                  >  2:797505/1‑150 (MQ=255)
                   cAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCCACAAAATTCAcc                                                                                                                                  >  2:1090673/1‑150 (MQ=255)
                        cGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtg                                                                                                                             >  1:1476843/1‑150 (MQ=255)
                           ataCTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCAGTACGCGAGTACAACTAAATGGTTGATGACTATAATAGCCTGCGATAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGa                                                                                                                          <  1:296665/148‑1 (MQ=255)
                            aaCTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCTGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                                                                         >  1:130603/1‑150 (MQ=255)
                                   gCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTc                                                                                                                  <  1:135177/150‑1 (MQ=255)
                                   gCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTc                                                                                                                  <  2:663083/150‑1 (MQ=255)
                                    cTGGATGGTCTAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTcc                                                                                                                 <  1:1604462/150‑1 (MQ=255)
                                       gATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATa                                                                                                              <  1:775760/150‑1 (MQ=255)
                                        attgTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATaa                                                                                                             <  2:435516/147‑1 (MQ=255)
                                        aTGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCTACACAGTTCACCCCTTGTTAAGTAACTCCATaa                                                                                                             >  1:1134647/1‑150 (MQ=255)
                                         tGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATaat                                                                                                            >  2:1027889/1‑149 (MQ=255)
                                         tGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATaaa                                                                                                            >  1:1554849/1‑150 (MQ=255)
                                           gTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATg                                                                                                          <  1:77128/150‑1 (MQ=255)
                                           gTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACACGTAACCCCATGTGAAGTAACTCCATAAATg                                                                                                          >  1:1055900/1‑150 (MQ=255)
                                             cGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                                                                                        <  1:1350373/150‑1 (MQ=255)
                                             cGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                                                                                        >  2:621401/1‑150 (MQ=255)
                                              gAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAAAGCCt                                                                                                       >  1:639899/1‑150 (MQ=255)
                                               aaGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGTTTGATGACTATAAAAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTc                                                                                                      <  1:1317662/150‑1 (MQ=255)
                                                  aCTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTc                                                                                                      >  2:1488411/1‑147 (MQ=255)
                                                  aCTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTc                                                                                                      <  1:1488411/147‑1 (MQ=255)
                                                   ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCatat                                                                                                  >  2:124686/1‑150 (MQ=255)
                                                      tttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCt                                                                                               <  1:678191/150‑1 (MQ=255)
                                                       ttCCCGATAGCATTACGCGAGTAAAACGAATTGGTTTATGACTATAATAGCCTGCGTTAAACGTTTTTGATTATAAATGCATTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                                                                              <  1:1090673/150‑1 (MQ=255)
                                                        tCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGCCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTTACTCCATAAATGCCTCATATCCTGt                                                                                             >  2:1514435/1‑150 (MQ=255)
                                                        tCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGCCACAGTTCACCCCTTGTGAAGTAACTCCATAAATGGCTCATATCCTGt                                                                                             >  1:445879/1‑150 (MQ=255)
                                                           cGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTaa                                                                                          >  2:1515227/1‑150 (MQ=255)
                                                                 cATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                                                                                    <  2:1314887/150‑1 (MQ=255)
                                                                      cgcgAGTACAACGAAATGTTTGATGTCTTTATAAGCCTGCGGTAATCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGTAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa                                                                               <  1:729547/150‑1 (MQ=255)
                                                                        cgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtc                                                                             <  1:293508/150‑1 (MQ=255)
                                                                          aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtctc                                                                           >  2:1640699/1‑150 (MQ=255)
                                                                            tACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg                                                                         >  2:169514/1‑150 (MQ=255)
                                                                              cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgcc                                                                       >  1:985298/1‑150 (MQ=255)
                                                                              aaaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgcc                                                                       <  2:1137383/149‑1 (MQ=255)
                                                                                aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCgg                                                                     >  2:938923/1‑150 (MQ=255)
                                                                                  gAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGAt                                                                   >  1:41306/1‑150 (MQ=255)
                                                                                     aTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCCGATGCGGCGCATCCCACACACTTCACCCCTTGTGAAGTAACTCCATAAATTCCTCATATCCAGTTAACCCGGCAATAATCTCCGACGGAAGCt                                                                >  2:1650144/1‑150 (MQ=255)
                                                                                        gTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac                                                             <  1:881277/150‑1 (MQ=255)
                                                                                        gTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTAACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCAGCCGGATGCTGac                                                             >  1:694652/1‑150 (MQ=255)
                                                                                          tgatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGacac                                                           >  2:224386/1‑150 (MQ=255)
                                                                                             tgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCCCCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAat                                                        >  2:1140505/1‑150 (MQ=255)
                                                                                                 tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGc                                                    <  2:158805/150‑1 (MQ=255)
                                                                                                 tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGc                                                    <  2:990396/150‑1 (MQ=255)
                                                                                                  ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCt                                                   >  1:12859/1‑150 (MQ=255)
                                                                                                  ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTAACCCCTTGTTAAGTAACTACATAAATGCCTCATATCCTGTTAACCCCGCAATAAACTCCGCCTGATGCTGAAACAaaattca                                                   >  1:153817/1‑147 (MQ=255)
                                                                                                     ataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGc                                                <  1:621401/150‑1 (MQ=255)
                                                                                                       aGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCa                                              <  2:1323828/150‑1 (MQ=255)
                                                                                                       aGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGAGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCTCCGGATGCTGACACAATATGCTTGCCa                                              >  2:1158883/1‑150 (MQ=255)
                                                                                                         ccTGCGTTATACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAATTAACTACATAAATGCCTCATATACTGTTAAACCGGCAATAATCTCCGCCGGATGCTTACAAAATATTCTTGCCata                                            >  2:562340/1‑148 (MQ=255)
                                                                                                         ccTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTACATAAATTCCTCATATCCGGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGa                                            >  2:569446/1‑150 (MQ=255)
                                                                                                            gCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAg                                         <  1:1640699/150‑1 (MQ=255)
                                                                                                               ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTa                                      >  2:9213/1‑150 (MQ=255)
                                                                                                                tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCCTCCGACACAGTTGACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCAGTTTACCCGGCAATAATCTCCGCCGGATGCTGGCACAATATGCTTGCCAGATAGCTaa                                     >  2:1585430/1‑150 (MQ=255)
                                                                                                                 aaaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTaaa                                    <  1:1383311/150‑1 (MQ=255)
                                                                                                                  aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAat                                   <  1:274419/150‑1 (MQ=255)
                                                                                                                       ttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATAtttt                              >  1:296619/1‑150 (MQ=255)
                                                                                                                       ttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATAtttt                              <  2:1085844/150‑1 (MQ=255)
                                                                                                                        tttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATAttttt                             >  1:10915/1‑150 (MQ=255)
                                                                                                                         ttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGAAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGCCACAATATGCTTGCCAGATAGCTAAATAttatta                            >  1:233856/1‑150 (MQ=255)
                                                                                                                           gatgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGc                          <  1:1056852/150‑1 (MQ=255)
                                                                                                                             tgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTc                        >  1:1017814/1‑150 (MQ=255)
                                                                                                                              gatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCa                       <  2:1528825/150‑1 (MQ=255)
                                                                                                                              gatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCa                       >  1:1022625/1‑150 (MQ=255)
                                                                                                                                tAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATg                     >  2:440388/1‑150 (MQ=255)
                                                                                                                                 aaaTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGa                    >  2:1157488/1‑150 (MQ=255)
                                                                                                                                  aaTGCAGTGTAGGATTCGGCGCATACGACACAGTTGACCCCATGTTAAGTAACTCCATAAATGCCTCATATACTGTTAACCCGGCAATAAAATACCCCGGATTCTGACACAATAATCTTGACATATATAAAAAAATTTTTATCACATTAc                   >  2:1380438/1‑150 (MQ=255)
                                                                                                                                        gtgtCGGATGCGGCGCATCCGACACAGTTCACCCGGTGTGAAGAAACTCCATAAAGGCCTCATATCCTGTTAGCCCGGCAATAATCTCCGCCGCATGCTGACACAATATGCTTGCCAGCTAGCTAAATATTTTTAGCTAATGACGCGCAt             >  2:685289/1‑150 (MQ=255)
                                                                                                                                        gtgtCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAACTGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCAt             >  1:286628/1‑150 (MQ=255)
                                                                                                                                        gtgtCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCGCCGCCGGATGCTGACACATGATGCTTGGCAGATAGCTAAATATTTTTAGCTCATGACGCGCAt             >  2:847538/1‑150 (MQ=255)
                                                                                                                                         tgtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCGAAATATTTTTAGCTCATGACGCGCATc            >  1:976463/1‑150 (MQ=255)
                                                                                                                                          gtCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcg           >  1:1523816/1‑150 (MQ=255)
                                                                                                                                          gtCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcg           >  2:868000/1‑150 (MQ=255)
                                                                                                                                          gtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGGCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcg           >  1:1085584/1‑150 (MQ=255)
                                                                                                                                          gtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcg           >  1:1631638/1‑150 (MQ=255)
                                                                                                                                          gtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcg           >  2:548381/1‑150 (MQ=255)
                                                                                                                                           tCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAAAAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcgg          <  2:1554849/150‑1 (MQ=255)
                                                                                                                                             ggATGGGGAGCATCCGACTAAGTTCACCCCTTTTGAAGTAACTCCATAAATGCCTCATAACGTGTTAACCCGTCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGAATGCTAAATATTTTAGCTCATGACGCGCATcggcgg       <  2:1016440/150‑1 (MQ=255)
                                                                                                                                                tGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcggcggcg     <  1:1370854/150‑1 (MQ=255)
                                                                                                                                                tGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATcggcggcg     <  2:1125100/150‑1 (MQ=255)
                                                                                                                                                   ggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGgcggcggcg  >  1:796046/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGGGATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCG  >  NC_000913/4087943‑4088239

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