Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,822,679 A→C 6.2% intergenic (‑40/+29) recX ← / ← recA regulatory protein for RecA/DNA recombination and repair protein; ssDNA‑dependent ATPase; synaptase; ssDNA and dsDNA binding protein; ATP‑dependent homologous DNA strand exchanger; recombinase A; LexA autocleavage cofactor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,822,6790AC 6.2% 672.3 / 32.6 195intergenic (‑40/+29)recX/recAregulatory protein for RecA/DNA recombination and repair protein; ssDNA‑dependent ATPase; synaptase; ssDNA and dsDNA binding protein; ATP‑dependent homologous DNA strand exchanger; recombinase A; LexA autocleavage cofactor
Reads supporting (aligned to +/- strand):  ref base A (85/98);  new base C (5/7);  total (90/105)
Fisher's exact test for biased strand distribution p-value = 7.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

Only 100 of 196 total aligned reads displayed.

CGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTT  >  NC_000913/2822531‑2822825
                                                                                                                                                    |                                                                                                                                                  
cgAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGa                                                                                                                                                    <  3:487971/149‑1 (MQ=255)
 gAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGAt                                                                                                                                                   >  3:546982/1‑149 (MQ=255)
   gTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGc                                                                                                                                                 >  3:101465/1‑149 (MQ=255)
        cGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGAccc                                                                                                                                            <  3:1470386/149‑1 (MQ=255)
          tCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGTCGCATCTGCGGCCCtt                                                                                                                                          <  3:125122/149‑1 (MQ=255)
           cGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGTCGCATCTGCGGCCCttt                                                                                                                                         <  2:1645615/149‑2 (MQ=255)
            gCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtgt                                                                                                                                        >  2:724557/1‑149 (MQ=255)
            gCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtgt                                                                                                                                        <  1:891185/149‑1 (MQ=255)
            gCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtgt                                                                                                                                        >  1:687972/1‑149 (MQ=255)
             cAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtgtg                                                                                                                                       <  2:508381/149‑1 (MQ=255)
              aGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTgtgt                                                                                                                                      >  4:91988/1‑149 (MQ=255)
              aGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTgtgt                                                                                                                                      <  4:1088884/149‑1 (MQ=255)
               gTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGTCGCATCTGCGGCCCTtttgct                                                                                                                                     <  1:380192/149‑3 (MQ=255)
               gTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTa                                                                                                                                     >  1:475859/1‑149 (MQ=255)
               gTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTa                                                                                                                                     >  2:92794/1‑149 (MQ=255)
                      cTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACa                                                                                                                              >  3:99918/1‑149 (MQ=255)
                       tCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGTATATGCCGCGCGACGGGATGTTAATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCCCAGACGCGTCCCTGGTGTATCAACCaa                                                                                                                             <  4:882561/149‑1 (MQ=255)
                       tCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACaa                                                                                                                             >  4:669865/1‑149 (MQ=255)
                       tCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACaa                                                                                                                             <  1:1695119/149‑1 (MQ=255)
                       tCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACaa                                                                                                                             >  1:896597/1‑149 (MQ=255)
                       tCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACaa                                                                                                                             >  2:466590/1‑149 (MQ=255)
                         aCTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATg                                                                                                                                                  <  2:836859/126‑1 (MQ=255)
                         aCTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATg                                                                                                                                                  >  1:836859/1‑126 (MQ=255)
                         aCTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGa                                                                                                                           >  4:1315211/1‑149 (MQ=255)
                         aCTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGa                                                                                                                           >  1:628661/1‑149 (MQ=255)
                         aCTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGa                                                                                                                           >  3:1556995/1‑149 (MQ=255)
                         aCTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGa                                                                                                                           <  4:1679498/149‑1 (MQ=255)
                               aTCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTa                                                                                                                     >  1:1410411/1‑149 (MQ=255)
                               aTCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTa                                                                                                                     >  4:706252/1‑149 (MQ=255)
                                  gcgcACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTaaaa                                                                                                                  <  2:428604/149‑1 (MQ=255)
                                   cgcACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCCGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTaaaaa                                                                                                                 <  4:1156941/149‑1 (MQ=255)
                                      aCCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGTCGCATCTGCGGCCCttt                                                                                                                                         <  2:1143513/122‑2 (MQ=255)
                                        cGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTc                                                                                                            <  3:1175042/149‑1 (MQ=255)
                                        cGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTc                                                                                                            <  4:1193035/149‑1 (MQ=255)
                                        cGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTc                                                                                                            <  2:275326/149‑1 (MQ=255)
                                          ccAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGt                                                                                                          <  4:764909/149‑1 (MQ=255)
                                             gAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtgt                                                                                                                                        <  2:88996/116‑1 (MQ=255)
                                                 gCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTc                                                                                                   <  1:610772/149‑1 (MQ=255)
                                                  cGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTCGATTCTGTCATGGCATATCCTTACAACATAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCt                                                                                                  <  3:1458786/149‑1 (MQ=255)
                                                   gTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTAt                                                                                                                                    >  2:1087504/1‑114 (MQ=255)
                                                   gTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTg                                                                                                 <  2:1058895/149‑1 (MQ=255)
                                                            gATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCtt                                                                                        >  1:548552/1‑149 (MQ=255)
                                                            gATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTCCTGCTACGCCtt                                                                                        >  1:1280975/1‑149 (MQ=255)
                                                             aTCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTc                                                                                       >  1:564373/1‑149 (MQ=255)
                                                               ccAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGc                                                                                     <  1:1613842/149‑1 (MQ=255)
                                                                 aaCAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTa                                                                                   <  2:896597/149‑1 (MQ=255)
                                                                 aaCAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTa                                                                                   >  3:1085481/1‑149 (MQ=255)
                                                                 aaCAGGCGAGCATATGCCGGGCGACGGGATGTGGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGAAGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTa                                                                                   >  4:1478190/1‑149 (MQ=255)
                                                                   cAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTatc                                                                                 <  2:592008/149‑1 (MQ=255)
                                                                    aGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTatca                                                                                <  4:1229550/149‑1 (MQ=255)
                                                                      gCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTatcatc                                                                              <  2:16078/149‑1 (MQ=255)
                                                                        gAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTa                                                                            <  1:429251/149‑1 (MQ=255)
                                                                          gCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACa                                                                          <  2:82928/149‑1 (MQ=255)
                                                                          gCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACa                                                                          <  4:1388913/149‑1 (MQ=255)
                                                                          gCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACa                                                                          <  4:1471048/149‑1 (MQ=255)
                                                                           cATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACag                                                                         >  4:627044/1‑149 (MQ=255)
                                                                            atatGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACaga                                                                        <  1:1661405/149‑1 (MQ=255)
                                                                             tatGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACagag                                                                       <  2:1410411/149‑1 (MQ=255)
                                                                              atGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAgaga                                                                      >  1:1536204/1‑149 (MQ=255)
                                                                                     gCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCgg                                                               <  4:546982/149‑1 (MQ=255)
                                                                                     gCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCgg                                                               <  1:924410/149‑1 (MQ=255)
                                                                                           ggATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTatca                                                                                >  4:1614318/1‑126 (MQ=255)
                                                                                           ggATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGa                                                         <  3:1373610/149‑1 (MQ=255)
                                                                                             aTGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAGAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGa                                                         <  3:719299/147‑1 (MQ=255)
                                                                                             aTGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGt                                                       <  2:1708774/149‑1 (MQ=255)
                                                                                               gTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTc                                                     >  2:1210800/1‑149 (MQ=255)
                                                                                               gTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTc                                                     >  1:1003753/1‑149 (MQ=255)
                                                                                                ttGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCg                                                    >  4:360792/1‑149 (MQ=255)
                                                                                                ttGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCg                                                    >  4:151811/1‑149 (MQ=255)
                                                                                                ttGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCg                                                    >  2:759512/1‑149 (MQ=255)
                                                                                                 tGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCgg                                                   <  3:9830/149‑1 (MQ=255)
                                                                                                 tGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCgg                                                   <  2:628661/149‑1 (MQ=255)
                                                                                                   aTTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGt                                                 <  4:1082175/149‑1 (MQ=255)
                                                                                                    ttCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGtt                                                >  3:1431010/1‑149 (MQ=255)
                                                                                                        gTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtgt                                                                                                                                        <  3:1199999/57‑1 (MQ=255)
                                                                                                           aTGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTcagc                                         >  3:935353/1‑149 (MQ=255)
                                                                                                              gCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCagcagc                                      <  3:266866/149‑1 (MQ=255)
                                                                                                                   tCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTc                                 <  3:1315211/149‑1 (MQ=255)
                                                                                                                      ttACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACg                              <  3:438528/149‑1 (MQ=255)
                                                                                                                      ttACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACg                              <  3:163886/149‑1 (MQ=255)
                                                                                                                      ttACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACg                              >  1:1132766/1‑149 (MQ=255)
                                                                                                                      ttACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACg                              >  3:657161/1‑149 (MQ=255)
                                                                                                                      ttACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACg                              >  4:438528/1‑149 (MQ=255)
                                                                                                                            cTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACttt                        <  2:564373/149‑1 (MQ=255)
                                                                                                                             ttAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTttc                       >  1:748305/1‑149 (MQ=255)
                                                                                                                             ttAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGAATAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTttc                       >  1:993922/1‑149 (MQ=255)
                                                                                                                                aaaaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTttcttc                    >  1:513885/1‑149 (MQ=255)
                                                                                                                                aaaaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTttcttc                    >  3:860089/1‑149 (MQ=255)
                                                                                                                                 aaaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTtcttct                   >  1:1710654/1‑149 (MQ=255)
                                                                                                                                 aaaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTtcttct                   <  2:687972/149‑1 (MQ=255)
                                                                                                                                  aaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTtctc                  >  2:1026484/1‑149 (MQ=255)
                                                                                                                                  aaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTtctc                  >  3:425132/1‑149 (MQ=255)
                                                                                                                                  aaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTtctc                  <  3:74261/149‑1 (MQ=255)
                                                                                                                                   aaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCg                 >  2:991058/1‑149 (MQ=255)
                                                                                                                                    aGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGa                <  4:245619/149‑1 (MQ=255)
                                                                                                                                      cAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGAtc              <  3:1034162/149‑1 (MQ=255)
                                                                                                                                      cAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGAtc              <  2:1021029/149‑1 (MQ=255)
                                                                                                                                       aaaaGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGAtct             >  1:1004837/1‑149 (MQ=255)
                                                                                                                                          aGGGTCGCATCTGCGGCCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGaaa                                                                    >  4:587193/1‑91 (MQ=255)
                                                                                                                                          aGGGTCGCATCTGCGGCCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGaaa                                                                    <  3:587193/91‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTT  >  NC_000913/2822531‑2822825

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