Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 28,116 T→A 6.2% L275Q (CTG→CAG)  rihC → ribonucleoside hydrolase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091328,1160TA 6.2% 226.8 / 4.4 65L275Q (CTG→CAG) rihCribonucleoside hydrolase 3
Reads supporting (aligned to +/- strand):  ref base T (29/32);  new base A (2/2);  total (31/34)
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.51e-01

GTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACA  >  NC_000913/27831‑28389
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                 
gTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGACTTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTggg                                                                                                                                                                                                                                                                                >  1:35566/1‑289 (MQ=255)
           gTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAgccagcca                                                                                                                                                                                                                                                                     <  2:416659/289‑1 (MQ=255)
                      tttGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGGAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTg                                                                                                                                                                                                                                                          >  1:314588/1‑289 (MQ=255)
                      tttGCATGTCACCAGTCAGGAAATATTAACTCCGGACTGTCTCTCTACCCTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGt                                                                                                                                                                                                                                                           <  1:112628/288‑1 (MQ=255)
                         ggATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCa                                                                                                                                                                                                                                                       <  2:294247/289‑1 (MQ=255)
                               cACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGAt                                                                                                                                                                                                                                                 <  2:35566/289‑1 (MQ=255)
                                            atatTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAg                                                                                                                                                                                                                                    >  2:163430/1‑289 (MQ=255)
                                                       tGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCAt                                                                                                                                                                                                                                                      <  2:220729/260‑1 (MQ=255)
                                                       tGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCAt                                                                                                                                                                                                                                                      >  1:220729/1‑260 (MQ=255)
                                                       tGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCcagcag                                                                                                                                                                                                                         >  2:390307/1‑289 (MQ=255)
                                                          cTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAgtg                                                                                                                                                                                                                       <  2:274481/288‑1 (MQ=255)
                                                              ctctctACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAgtgggtg                                                                                                                                                                                                                   >  2:15947/1‑288 (MQ=255)
                                                                 tctACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTg                                                                                                                                                                                                               <  1:106852/289‑1 (MQ=255)
                                                                      acTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTg                                                                                                                                                                                                          <  1:15947/289‑1 (MQ=255)
                                                                           cGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGaa                                                                                                                                                                                                                                      >  1:146623/1‑256 (MQ=255)
                                                                           cGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGaa                                                                                                                                                                                                                                      <  2:146623/256‑1 (MQ=255)
                                                                           cGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTgg                                                                                                                                                                                                      >  2:68448/1‑288 (MQ=255)
                                                                           cGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTgg                                                                                                                                                                                                      >  2:376327/1‑288 (MQ=255)
                                                                           cGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCAGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTgg                                                                                                                                                                                                      >  2:458278/1‑288 (MQ=255)
                                                                           cGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCAGGGCAAGCCAGCCAATGGACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTg                                                                                                                                                                                                       >  2:458274/1‑287 (MQ=255)
                                                                                           cGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAAc                                                                                                                                                                                       <  2:159198/289‑1 (MQ=255)
                                                                                                aaaTGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTcac                                                                                                                                                                                >  2:163315/1‑289 (MQ=255)
                                                                                                aaaTGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTcac                                                                                                                                                                                <  2:390905/289‑1 (MQ=255)
                                                                                                                  ttAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGt                                                                                                                                                              <  1:163315/289‑1 (MQ=255)
                                                                                                                  ttAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGt                                                                                                                                                              <  1:163430/289‑1 (MQ=255)
                                                                                                                     gCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGCGGTCGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGtcat                                                                                                                                                           >  2:115781/1‑289 (MQ=255)
                                                                                                                         cTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGtcattca                                                                                                                                                        >  1:423131/1‑288 (MQ=255)
                                                                                                                           aCCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCg                                                                                                                                                     >  1:143776/1‑289 (MQ=255)
                                                                                                                                   ggCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAgccagcca                                                                                                                                                                                                                                                                     >  1:273285/1‑169 (MQ=255)
                                                                                                                                   ggCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAgccagcca                                                                                                                                                                                                                                                                     <  2:273285/169‑1 (MQ=255)
                                                                                                                                         aTGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGTTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCttt                                                                                                                                         <  2:303053/287‑1 (MQ=255)
                                                                                                                                         aTGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCttt                                                                                                                                         <  2:143776/287‑1 (MQ=255)
                                                                                                                                              aggCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTAAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGTGTGGCTGCGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCACTCATCGACTCATGCCTTTCACTGat                                                                                                                                  <  1:261790/287‑1 (MQ=255)
                                                                                                                                                  ggCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATAtc                                                                                                                               >  1:469527/1‑288 (MQ=255)
                                                                                                                                                  ggCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATAtc                                                                                                                               <  2:423131/288‑1 (MQ=255)
                                                                                                                                                    cTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAag                                                                                                                                                                                                                                                                            >  2:378009/1‑145 (MQ=255)
                                                                                                                                                    cTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAag                                                                                                                                                                                                                                                                            <  1:378009/145‑1 (MQ=255)
                                                                                                                                                     ttGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATAtccctc                                                                                                                           <  2:314588/289‑1 (MQ=255)
                                                                                                                                                                              tCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTa                                                                                                  >  1:302639/1‑289 (MQ=255)
                                                                                                                                                                              tCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTa                                                                                                  >  2:385397/1‑289 (MQ=255)
                                                                                                                                                                                  ctggctggTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCa                                                                                               <  1:390307/288‑1 (MQ=255)
                                                                                                                                                                                           gcgcCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCa                                                                                                                                                                                                                                                       >  2:363689/1‑127 (MQ=255)
                                                                                                                                                                                           gcgcCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCa                                                                                                                                                                                                                                                       <  1:363689/127‑1 (MQ=255)
                                                                                                                                                                                            cgcCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAAc                                                                                                                                                                                       <  2:73624/190‑1 (MQ=255)
                                                                                                                                                                                            cgcCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAAc                                                                                                                                                                                       >  1:73624/1‑190 (MQ=255)
                                                                                                                                                                                                         tCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGTATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCg                                                                       <  2:505583/289‑1 (MQ=255)
                                                                                                                                                                                                         tCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCg                                                                       <  2:302639/289‑1 (MQ=255)
                                                                                                                                                                                                         tCA‑‑CTCTCAAACCCTGTTTGATGGCAGTGTAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATGGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCg                                                                       <  1:68448/289‑1 (MQ=255)
                                                                                                                                                                                                               tttcAAACCCGGTTGTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAAGCACGACGGTGGGTGATATCGACGGTGGCCTGGGCAAGCCAGCCAAAGTACAGGTGGCATTGGATCTGGGTGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTGATTTCTAATTATCAGCGTTTTTGGCTAGCGGCGTAGCGATgc                                                                   <  1:115781/287‑1 (MQ=255)
                                                                                                                                                                                                               tctcAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATg                                                                    >  1:435270/1‑288 (MQ=255)
                                                                                                                                                                                                                    aaCCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTa                                                                                                             >  2:120715/1‑242 (MQ=255)
                                                                                                                                                                                                                    aaCCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTa                                                                                                             <  1:120715/242‑1 (MQ=255)
                                                                                                                                                                                                                               ggggCAGTGGAAACGCAGGGCGAATTTACCTCAGGGACGACGGTGGTTGATATCGACGGTTGCCAGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGGGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaa                                                   <  1:458278/287‑1 (MQ=255)
                                                                                                                                                                                                                                tgGAAGTGGGAACTCAGGGCGAATTTACCTCAGGCACGCCGTTGGTTGATATCGACGGTTGCCAGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaa                                                   <  1:458274/288‑1 (MQ=255)
                                                                                                                                                                                                                                                gggCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTAc                                                                                                                                                                                                                                                               >  2:317673/1‑68 (MQ=255)
                                                                                                                                                                                                                                                gggCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTAc                                                                                                                                                                                                                                                               <  1:317673/68‑1 (MQ=255)
                                                                                                                                                                                                                                                gggCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTa                                  >  1:233893/1‑289 (MQ=255)
                                                                                                                                                                                                                                                           ccTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCGTTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                       <  2:412487/289‑1 (MQ=255)
                                                                                                                                                                                                                                                           ccTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                       >  1:376472/1‑289 (MQ=255)
                                                                                                                                                                                                                                                               aGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATa                     >  2:190703/1‑287 (MQ=255)
                                                                                                                                                                                                                                                                         gtggtTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATa                     >  1:286016/1‑277 (MQ=255)
                                                                                                                                                                                                                                                                         gtggtTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATa                     <  2:286016/277‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 tatCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACa  >  1:503089/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                            gCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                       <  1:4523/256‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            gCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                       >  2:4523/1‑256 (MQ=255)
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                 
GTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACA  >  NC_000913/27831‑28389

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