Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 994,950 T→A 7.9% intergenic (+26/+99) PP_0859 → / ← PP_0860 ketoglutaramate omega‑amidase/sulfite reductase flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947994,9500TA 7.9% 335.7 / 14.1 113intergenic (+26/+99)PP_0859/PP_0860ketoglutaramate omega‑amidase/sulfite reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (61/43);  new base A (5/4);  total (66/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 115 total aligned reads displayed.

GAGCCTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTACAGC  >  NC_002947/994801‑995099
                                                                                                                                                     |                                                                                                                                                       
gagCCTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCt                                                                                                                                                         >  1:727034/1‑150 (MQ=255)
   ccTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGc                                                                                                                                                      <  2:2028547/150‑1 (MQ=255)
       cTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCcgc                                                                                                                                                  <  1:2049890/150‑1 (MQ=255)
        tCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCcgcg                                                                                                                                                 >  1:1462835/1‑150 (MQ=255)
              ccGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAggg                                                                                                                                           >  1:1433530/1‑150 (MQ=255)
                     ggCGGACGGGGTGTTCACTGCAGGACTGAGCGCGGCGGAACTGGCGGCGCACCGGGCCAAGTTACAGTAGAGGCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGccc                                                                                                                                    <  1:1311237/150‑1 (MQ=255)
                            ggggTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTgccagc                                                                                                                             >  2:1563645/1‑150 (MQ=255)
                              ggTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAg                                                                                                                                             >  2:1328683/1‑132 (MQ=255)
                              ggTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAg                                                                                                                                             <  1:1328683/132‑1 (MQ=255)
                              ggTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTgccagcc                                                                                                                            <  2:555152/149‑1 (MQ=255)
                              ggTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTgccagcc                                                                                                                            >  1:555152/1‑149 (MQ=255)
                                 gtTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAAtt                                                                                                                        <  1:826172/150‑1 (MQ=255)
                                   tCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGGCCGGCCCTTTCGCGGGCAAGCCCGCGAAGAAGCCGGGCCTgctgtcagt                                                                                                                          >  1:729673/1‑139 (MQ=21)
                                   tCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGGCCGGCCCTTTCGCGGGCAAGCCCGCGAAGAAGCCGGGCCTgctgtcagt                                                                                                                          <  2:729673/146‑8 (MQ=21)
                                   tCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTc                                                                                                                      <  1:464357/150‑1 (MQ=255)
                                    cACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTcc                                                                                                                     >  1:396088/1‑150 (MQ=255)
                                             aCTGGGCGCGGCGGAACTGGCGGCGTACCGGCCCCAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCCGCTTCTTCGCGGGCTTCCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGt                                                                                                            <  1:94754/150‑1 (MQ=255)
                                              cTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACGGGcgggt                                                                                                           >  1:1150185/1‑146 (MQ=255)
                                                 aGCGCGGCGGAACTGGCTGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CAc                                                                                                      >  1:27787/1‑150 (MQ=255)
                                                 aGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CAc                                                                                                      >  2:2214025/1‑150 (MQ=255)
                                                   cgcgGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCg                                                                                                    >  2:2188852/1‑150 (MQ=255)
                                                   cgcgGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCg                                                                                                    <  1:1563645/150‑1 (MQ=255)
                                                      ggcggAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTg                                                                                                 >  2:725081/1‑150 (MQ=255)
                                                          gAACTGGCGGCGTACCGGGCCAAGTTGCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTgccgc                                                                                             >  2:30298/1‑150 (MQ=255)
                                                              tGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGc                                                                                                                                                      <  2:487807/91‑1 (MQ=255)
                                                              tGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGc                                                                                                                                                      >  1:487807/1‑91 (MQ=255)
                                                              tGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAACAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTgccgccgcc                                                                                         >  2:1613540/1‑150 (MQ=255)
                                                                gcggcgTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTgccagccc                                                                                                                           >  1:1967629/1‑115 (MQ=255)
                                                                gcggcgTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTgccagccc                                                                                                                           <  2:1967629/116‑2 (MQ=255)
                                                                     gTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTAt                                                                                  >  2:807081/1‑150 (MQ=255)
                                                                      tACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTgcca                                                                                                                               <  2:574782/106‑1 (MQ=255)
                                                                      tACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTgcca                                                                                                                               >  1:574782/1‑106 (MQ=255)
                                                                      tACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTAtt                                                                                 <  1:923905/150‑1 (MQ=255)
                                                                        ccGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTcg                                                                               >  1:312320/1‑150 (MQ=255)
                                                                        ccGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTcg                                                                               >  1:2046423/1‑150 (MQ=255)
                                                                         cGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTcgc                                                                              <  2:2046423/150‑1 (MQ=255)
                                                                             ccAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGt                                                                          <  1:807081/150‑1 (MQ=255)
                                                                              cAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccgg                                                                                                                                       >  2:1729546/1‑90 (MQ=255)
                                                                              cAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccgg                                                                                                                                       <  1:1729546/90‑1 (MQ=255)
                                                                                 gTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAAc                                                                      <  2:2066467/150‑1 (MQ=255)
                                                                                   tCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAAccc                                                                    >  2:1372121/1‑150 (MQ=255)
                                                                                     ccGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGTGCCGGCCCTGCAAGCCCATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGc                                                                  >  2:1717445/1‑150 (MQ=255)
                                                                                     ccGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGACCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCggg                                                                           >  1:1760303/1‑141 (MQ=255)
                                                                                     ccGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGACCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCggg                                                                           <  2:1760303/141‑1 (MQ=255)
                                                                                     ccGGCCAATCTGGATGCCGATACATTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGc                                                                  >  2:2271718/1‑150 (MQ=255)
                                                                                      cGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCt                                                                 >  2:793390/1‑150 (MQ=255)
                                                                                          cAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTa                                                             <  1:2163621/150‑1 (MQ=255)
                                                                                          cAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTa                                                             >  1:244314/1‑150 (MQ=255)
                                                                                             tCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTacac                                                          >  2:1373804/1‑150 (MQ=255)
                                                                                             tCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTacac                                                          >  2:923081/1‑150 (MQ=255)
                                                                                             tCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGCCTTACCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTacac                                                          >  2:1373746/1‑150 (MQ=255)
                                                                                             tCTGGATGCCGATACCTTCGAGCTGCACTGACAACAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTacac                                                          <  1:1613540/150‑1 (MQ=255)
                                                                                                ggATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAcc                                                       >  2:1485315/1‑150 (MQ=255)
                                                                                                ggATGCCGATACCTTCGAGCTGCACGGACAGCAGGCCCGGCTGCTTCGCGGGCTTGCCCGCGACAGGGCCGGCCCTGCCAGCAAATTACCCAGGCAGGTGCCATCTCTGATCTGCTACCGGAGTCGCGAGTAAACCCGCTCCTACACAcc                                                       >  2:1093581/1‑150 (MQ=255)
                                                                                                    gCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGAt                                                    <  1:1252864/149‑1 (MQ=255)
                                                                                                    gCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGAt                                                    >  2:1252864/1‑149 (MQ=255)
                                                                                                    gCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATc                                                   <  1:1485315/150‑1 (MQ=255)
                                                                                                    gCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATc                                                   <  2:1462835/150‑1 (MQ=255)
                                                                                                    gCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATc                                                   <  1:967421/150‑1 (MQ=255)
                                                                                                     ccGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccggcc                                                                                                                                     <  1:2336260/69‑1 (MQ=255)
                                                                                                     ccGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccggcc                                                                                                                                     >  2:2336260/1‑69 (MQ=255)
                                                                                                      cGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCCCCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCaa                                                 >  2:72160/1‑150 (MQ=255)
                                                                                                      cGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCaa                                                 >  2:407797/1‑150 (MQ=255)
                                                                                                      cGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCaa                                                 <  1:725081/150‑1 (MQ=255)
                                                                                                          aCCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTaca                                                           <  1:2108918/136‑1 (MQ=255)
                                                                                                            cTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATAcaca                                           >  1:2150558/1‑150 (MQ=255)
                                                                                                             ttCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CAc                                                                                                      <  2:483095/90‑1 (MQ=255)
                                                                                                             ttCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CAc                                                                                                      >  1:483095/1‑90 (MQ=255)
                                                                                                             ttCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTc                                                                <  1:562799/128‑1 (MQ=255)
                                                                                                             ttCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTc                                                                >  2:562799/1‑128 (MQ=255)
                                                                                                             ttCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACAt                                          >  1:228261/1‑150 (MQ=255)
                                                                                                             ttCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACAt                                          <  1:923081/150‑1 (MQ=255)
                                                                                                              tCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATCCGCGGGTAAACCCGCTCCTACACTCCGATCAATACACAtc                                         <  2:336422/150‑1 (MQ=255)
                                                                                                                  gCTGCACTCACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGc                                     <  1:899136/150‑1 (MQ=255)
                                                                                                                   cTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGcc                                    >  1:2198305/1‑150 (MQ=255)
                                                                                                                    tGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCg                                   >  2:917034/1‑150 (MQ=255)
                                                                                                                     gaaCTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGa                                  >  2:1843981/3‑150 (MQ=255)
                                                                                                                     gCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGAt                                                    >  2:2246635/1‑132 (MQ=255)
                                                                                                                     gCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGAt                                                    <  1:2246635/132‑1 (MQ=255)
                                                                                                                      cACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGat                                 >  2:784875/1‑150 (MQ=255)
                                                                                                                         tGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccggcc                                                                                                                                     <  1:344016/49‑1 (MQ=255)
                                                                                                                         tGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccggcc                                                                                                                                     >  2:344016/1‑49 (MQ=255)
                                                                                                                         tGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATc                              <  1:1372121/150‑1 (MQ=255)
                                                                                                                            cagcagGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGcc                           >  1:493557/1‑150 (MQ=255)
                                                                                                                                    ccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCg                   >  1:809720/1‑150 (MQ=255)
                                                                                                                                     cGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGa                  >  2:1240426/1‑150 (MQ=255)
                                                                                                                                        cttcttCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGa                  >  2:354990/1‑147 (MQ=255)
                                                                                                                                        cttcttCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGa                  <  1:354990/147‑1 (MQ=255)
                                                                                                                                         ttcttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCaa              >  2:1004448/1‑150 (MQ=255)
                                                                                                                                         ttcttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAAAAAGATCATCCACACCGCTACCTACACTTCTCGCCCAGATCGCCCATCCTCGATCaa              >  2:1192676/1‑150 (MQ=21)
                                                                                                                                           tttcGCGGGCAAGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGAt                                                    <  2:2024580/109‑1 (MQ=255)
                                                                                                                                           tttcGCGGGCAAGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATacac                                            >  1:2048748/2‑118 (MQ=255)
                                                                                                                                           tttcGCGGGCAAGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATacac                                            <  2:2048748/117‑1 (MQ=255)
                                                                                                                                           tttcGCGGGCAAGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAAct            >  2:1379044/2‑150 (MQ=255)
                                                                                                                                           tttcGCGGGCAAGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAAct            >  2:1895426/2‑150 (MQ=255)
                                                                                                                                            ttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActg           >  2:101466/1‑150 (MQ=255)
                                                                                                                                            ttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActg           >  1:2008266/1‑150 (MQ=255)
                                                                                                                                              cgcgGGGTTGCCGGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActgct         <  2:727034/150‑1 (MQ=255)
                                                                                                                                                    cTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTAcag   <  2:789488/150‑1 (MQ=255)
                                                                                                                                                     ttGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTAcagc  >  2:638573/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
GAGCCTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTT‑‑CACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTGTACAGC  >  NC_002947/994801‑995099

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