Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,626,601 G→T P230P (CCG→CCT radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,626,6010GT96.9% 185.5 / NA 65P230P (CCG→CCTradADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (25/37);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGT  >  NC_000913/4626364‑4626846
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         
cccgtccgatatggAAACCAGCATCGTGCAGAGCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGAGGGACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAac                                                                                                                                                                                                                                             <  2:200376/239‑1 (MQ=255)
         gtCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATc                                                                                                                                                                                                                                    <  1:111345/250‑1 (MQ=255)
            ggAAACCAGCATCGCGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTTAAGGGCTGGAacagaacgggag                                                                                                                                                                                                                                   >  1:200376/1‑240 (MQ=255)
            ggAAACCAGCATCGAGCACATCTGCCTGTGTGCAGAAGAAGAGCAACCGAAGCTGAGGGTAATTGACGCGACCCAGTGGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTTAAACTGCGGCTTATCTTACGCGCTTCCCCAAAACGCGCGGTGTGGCGAGTGTCATGGTTGGACTCGTAACCAAAGATGGCTCGCTGGCTGGCCCTAAAGTGCTGGAACACTGTATCgcct                                                                                                                                                                                                                                <  2:425107/251‑4 (MQ=255)
                    gCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTg                                                                                                                                                                                                                        <  1:338391/251‑1 (MQ=255)
                      aTCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGGGATGCAAATGGCGGATGGACAGTCATAGCCTGGCAGCGTGGCGCAGGGGCGGGAAACGGAGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGAGGGGCGATTGTCATGGGTGGAAACGTAACCAAAGAGGGCTCGCTGGCGGGCCCTAAAGAGCTGGAAAACTGTATCGACAGTTCGGGGCt                                                                                                                                                                                                                      >  1:176099/1‑251 (MQ=255)
                              gATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTa                                                                                                                                                                                                                                      <  2:235181/227‑1 (MQ=255)
                              gATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTa                                                                                                                                                                                                                                      >  1:235181/1‑227 (MQ=255)
                                   gCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGCCTGTTCGGTGCTTTTGGATGGCg                                                                                                                                                                                                          >  1:76635/1‑249 (MQ=255)
                                      agATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGcc                                                                                                                                                                                                     <  1:40660/250‑1 (MQ=255)
                                        aTTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTAACGCATGTTCGGTGCTTTTGGATGGCGATGCCgg                                                                                                                                                                                                   >  1:347776/1‑250 (MQ=255)
                                                   agagCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTc                                                                                                                                                                                        <  1:148646/251‑1 (MQ=255)
                                                    gagCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCg                                                                                                                                                                                       <  1:447822/251‑1 (MQ=255)
                                                      ccaccgtaaGCTGATGGTAATTGACGCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTACCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCa                                                                                                                                                                                     <  1:127369/244‑1 (MQ=255)
                                                          ccGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGGGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGGTTTTGGATGGCGATGCCGAATCCCGGTTTCGCActtt                                                                                                                                                                                 >  1:423422/1‑248 (MQ=255)
                                                            gAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTTGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTGCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCActttac                                                                                                                                                                               <  1:104364/251‑6 (MQ=255)
                                                              aGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCCAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCCGGCGGGGCCCAAAATGCTGGAAAACTGTATCGCCTTTTTGGTGCTTTTTGGTGGCGATGCCGACTCCCCTTTTCGCActtttcgc                                                                                                                                                                             >  2:144475/1‑244 (MQ=255)
                                                                  gATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGcc                                                                                                                                                                         <  1:380788/251‑1 (MQ=255)
                                                                       tAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGGTTGTCATGGTTGGCCACGTAACCAAAGATGGCTCGCGGGCGGGCCCTAAAGTGCTGGAACACTGTATAGACTGTTTGGGGGGTTTGGGTGGGGATGCCGACTCCCCtttttgccctttaccccccccaca                                                                                                                                                                     >  2:9160/1‑230 (MQ=255)
                                                                            gACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAAcc                                                                                                                                                                <  2:305367/250‑1 (MQ=255)
                                                                                        gTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGt                                                                                                                                                   <  1:298075/251‑1 (MQ=255)
                                                                                        gTGATGCATATGGCGGATGTACAGACATCGCCTGGCAGCGGGGAGCAGGTGCGGGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACCCGTAACCAAAGATGGCCCGCTGGCGGGCCCTTAAGTGCTGGGACACTGGTACGCATGTTCGGTGCCTTTGGAAGGGCGATCAGCCTCCCGTTTTCCCAACTTACGCCGCCATAAAAACCCCttaggggaggg                                                                                                                                                   >  2:396544/1‑242 (MQ=255)
                                                                                           aTGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGaa                                                                                                                                                <  1:35723/251‑1 (MQ=255)
                                                                                             gCATAGGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATg                                                                                                                                              <  2:372209/251‑1 (MQ=255)
                                                                                               atatGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACCCGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGGATCGACTGTTCGGGGCTATTGAATGGCGTTGCCGACTCCCGGTTTCGCACTTTACGCAGCCATAAAAACCGCTTCAGGGCGGTGAATg                                                                                                                                              >  1:457861/1‑249 (MQ=255)
                                                                                                 atGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGGGTGGCGATTGTCATGGTTGGAAACGTAACCAAAGATGGCTCGCTGGGGGGCCCTAAAGAGGGGGAGCACTGTTTCGACTGGTCGGGGCGTTTGGATGGGGAGGCCGGCTCCCGTTTTTGCAATTTACGCCGCCCTAAAAACCCCTTCgggccgggggatgggat                                                                                                                                          >  2:292824/1‑236 (MQ=255)
                                                                                                     cGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCTGGAACCCTGGATCGACTGTTCGGTGCTTTTGGGGGGCGGTGCCGACCCCCGTTTTCGCACTTTACGCAGCCATAAAAACCggctcgggcgggggaatgggctgggg                                                                                                                                      >  2:45882/1‑226 (MQ=255)
                                                                                                            accgTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTAACTGACTCGCTTCTCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGTTCGCTGGCGGGCCCTAAAGTTCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATGAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTcg                                                                                                                               <  2:176099/248‑1 (MQ=255)
                                                                                                                       cTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGa                                                                                                                    >  1:51192/1‑250 (MQ=255)
                                                                                                                              gggggCGCTGGGGCGTGAAACGGTGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTTGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGc                                                                                                            <  2:16368/248‑1 (MQ=255)
                                                                                                                                tGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATAGTCATGGATGGACACGTAACCAAAGATGGGTCGCTGGCGGGCCCTATAGTGCTGGACCACTGTATCGACTGGTCGGGGCTTTTGGATGGCGATGCCTACTCCCCGTTTCTCACTTTTCGCAGCCATAAAAAACGCTTCGGCGCGGTGGATGAGACGTGCGACTTCGCCTAGGACCGAACAgggctg                                                                                                           >  2:293936/1‑248 (MQ=255)
                                                                                                                                    gcAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTg                                                                                                      >  1:137365/1‑251 (MQ=255)
                                                                                                                                       ggTGCATGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGATTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAg                                                                                                   <  1:9771/251‑1 (MQ=255)
                                                                                                                                           cGTGAAACGGCGGCTTATCGTACGCGCTTCGCCAAATCGCGCGGTGTGTCTATTTTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTGCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCTGCCATAACAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAg                                                                                               <  1:38551/251‑1 (MQ=255)
                                                                                                                                             tGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGACCACTGTATCGCCTGTTCGGTGCTTTTGGGGGGCGATGCCGACTCCCGTTTTCGCCCTTTACGCAGCCATAAAAACCGCTTCGGCGCGGGGGCTGAGCTGGGCGTCTTCGCGATGCCCGAACAGGGGCTGCGTggaggcaac                                                                                              >  2:359403/1‑242 (MQ=255)
                                                                                                                                             tGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCa                                                                                                <  2:197324/248‑1 (MQ=255)
                                                                                                                                             tGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCa                                                                                                >  1:197324/1‑248 (MQ=255)
                                                                                                                                                   ggcggcTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGTAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCtt                                                                                       <  1:175006/251‑1 (MQ=255)
                                                                                                                                                    gtggcttacctGACGCGCCTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTTGACACGTAACCAAAGATGGCGCGCTGGCGTGCCCTAAAGTGCTGGAACACTGGATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTTAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTc                                                                                      <  2:293833/241‑1 (MQ=255)
                                                                                                                                                    gcggcTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTc                                                                                      <  1:200037/251‑1 (MQ=255)
                                                                                                                                                      ggtttagctGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTGGCGGGCGGGCCCTAAAGTGCTGGAACACGGTATCGACTGTGCGGTGCTTTTTGATGGCGATGCCGACTCCCGTTGTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTGCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCgg                                                                                    <  2:162350/243‑1 (MQ=255)
                                                                                                                                                               tacgcgattCGCCAAAACTCGCGGTGGGGCGATTGTCATGGGTGGATACGGAACCAAAGATGGTCCGCTGGTGGGCCGTAAAGTGCTGGACCACTGTATCGACTGTACGTTGTTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGTCCGAACAGGGTCTGCGTGAATACAGCAATATTTCGGCAATTTTc                                                                            <  2:394924/243‑1 (MQ=255)
                                                                                                                                                               gacgcGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAAtt                                                                               <  1:241998/243‑1 (MQ=255)
                                                                                                                                                                  gcgcgtCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGGTGGGGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCGTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTaa                                                                        <  1:61136/246‑1 (MQ=255)
                                                                                                                                                                   cGCTTCGCAAAAACGCGGGGTGTGGCGATTGTCATGGTGGGGCGCGTACCCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACGCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAg                                                                       <  2:114656/251‑1 (MQ=255)
                                                                                                                                                                     cTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGCCTGCTCGGTGCCTTTGGATGGCGCTGCCGACTACCGTTTTCGCCATTTACGCACACATAAAAACCACTTCGGCGCGGGTAATGAGCAGGACTGCTTCGCGATAGCCCACCCGGGGCTGCGTGATGTAAGACACACTTAGGCAATTTTCTTTAGTc                                                                     >  2:405223/1‑251 (MQ=255)
                                                                                                                                                                       tCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTc                                                                     <  1:144475/249‑1 (MQ=255)
                                                                                                                                                                          ccAAAACGCGCGGTGTGGCGATTGTAATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCgcgga                                                                <  2:331747/251‑2 (MQ=255)
                                                                                                                                                                            aaaaCGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGCTGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCCCGATGGCCGAACAGGGGCTCCGTGAAGTCAGCAACCCTTCGGACATTTTTTTAAGTCgcgggg                                                               >  2:389398/1‑248 (MQ=255)
                                                                                                                                                                                 ggggggTGTGGGGATTGTCATGGTTGGACACGGAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAg                                                         <  2:9683/247‑1 (MQ=255)
                                                                                                                                                                                 gggcgGTGTGGCGATTGTCAGGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCTCTAAAGTGCTGGAACACTGTATCGACGGTTCGGTGCTTTTGGAGGGCGATGCCGACTCCCGTTTGCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAg                                                         <  2:469532/249‑1 (MQ=255)
                                                                                                                                                                                  cgcgGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAgg                                                        >  1:413023/1‑250 (MQ=255)
                                                                                                                                                                                  cgcgGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGAGGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGGACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGt                                                        <  1:32672/251‑1 (MQ=255)
                                                                                                                                                                                     ggTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGAc                                                     <  1:221702/251‑1 (MQ=255)
                                                                                                                                                                                              aTTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCTGGAACCCTGGATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCCGCCATAAAAAACGCTTCGGCGCGGTGAATGGGCTGGGCGTCTTCGCGAAGACCGAAAAGGGGCTGCGGGACGTCAGCAACCCTTCGGCAATTTTCTTAAGGCgcggagtgggagggccctccggccg                                            >  2:414251/1‑230 (MQ=255)
                                                                                                                                                                                                tGTCATGGTTGGACACGTAACCAACGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTGTCTTAATTCGCGGAGATGAAGTGACCTCCGGCAGc                                           <  1:483537/250‑1 (MQ=255)
                                                                                                                                                                                                tGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGa                                                      >  1:102891/1‑239 (MQ=255)
                                                                                                                                                                                                tGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGa                                                      <  2:102891/239‑1 (MQ=255)
                                                                                                                                                                                                      ggttggacACGTAACCAAAGATGGCTCGCTGGCGGGCCCTCAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGGGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGGGATGAAGTGCCCTCCGGCAGCTCCGTg                                     >  2:473478/8‑250 (MQ=255)
                                                                                                                                                                                                      ggttggacACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTAGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGCCTTCGCGATGACCGAACCGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCCTAAGTCGCGGAAGTGAAGTGACCTCCGGCAGCTccgtt                                     >  2:404505/8‑246 (MQ=255)
                                                                                                                                                                                                          ggacACGTAACCAAAGATGGCTCGGTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTGAATTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGTtggt                                <  1:130828/248‑1 (MQ=255)
                                                                                                                                                                                                                          aTGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACAATGTATCGCCTGTTCGGAGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCCTAAAAACCGCTTCGACGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGGTGGAGGTGTGGGAGGGAAcg                >  1:425064/1‑251 (MQ=255)
                                                                                                                                                                                                                            ggctcggtggcgggcCCTAAAGGGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGATTCACGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGGCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAAcgcg              <  2:302943/244‑1 (MQ=255)
                                                                                                                                                                                                                                        ggcCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCtggt  >  1:91866/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ggcCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTGTTCGCACTTTACGCAGCCATAAAACCCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCtggt  <  1:120539/251‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         
CAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGT  >  NC_000913/4626364‑4626846

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