Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,706,790 C→A 10.7% T148K (ACG→AAG)  rsxB → SoxR iron‑sulfur cluster reduction factor component, putative iron‑sulfur protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,706,7900CA10.7% 178.5 / 15.5 67T148K (ACG→AAG) rsxBSoxR iron‑sulfur cluster reduction factor component, putative iron‑sulfur protein
Reads supporting (aligned to +/- strand):  ref base C (22/38);  new base A (4/3);  total (26/41)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.48e-01

AAAATCAACCGTTGCGCCCCAGGTGGCGAAGCTGTGATG‑CTAAAAATTGCCGAGTTGCTTAATGTCGAGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCA  >  NC_000913/1706555‑1707025
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                            
aaaaCCACCCGTTGCGCCCAAGGGGGGGAAGGTGGGATG‑TTAAAAATGCCCGAGTTGGTTAAGGGCGGGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACggg                                                                                                                                                                                                                                           <  2:1003536/239‑1 (MQ=255)
    tCACCCGTTGCCCCCCAGGTGGCGAAGCTGTGATG‑CTAAAAATTGCCGAGTTGCTTAATGTCGAGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGc                                                                                                                                                                                                                                       <  1:994761/239‑1 (MQ=255)
    ccccccGTTTCCCCCCGGGTGGGGAAGCT‑TGGTGTGTAAAAATTGCCGGGTTGCTTAATGTGGAGCCGCCGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGc                                                                                                                                                                                                                                       <  2:724693/235‑1 (MQ=255)
    aattttgtataagaggcagTGGGGAAGCTTGGATT‑GTAAAAATTGCCGATTTGCTTAATGTCGAGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGc                                                                                                                                                                                                                                       <  2:1113314/221‑1 (MQ=255)
         cgtgtcgcCCCAGGTGGCGAAGCTGTGATG‑CTAAAAATGTCCGAGTTGCTTAATGTCGAGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAAttt                                                                                                                                                                                                                                  <  1:1313446/234‑1 (MQ=255)
               gcCCCAGGTGGCGAAGCTGTGATG‑CTAAAAATTGCCGAGTTGCTTAATGTCGAGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATgtgt                                                                                                                                                                                                                            <  1:381476/239‑1 (MQ=255)
                       tGGCGAAGCTGTGATG‑CTAAAAATTGCCGAGTTGCTTAATGTCGAGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATcc                                                                                                                                                                                                                      <  1:1172238/237‑1 (MQ=255)
                                  tgATG‑CTAAAAATTGCCGAGTTGCTTAATGTCGAGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTTGCGCTACCCGTGCCATGCCTACGGTAATGAGGGGTCTCTGTACGGGCTGCAATTTTTttgtctgtctctttttaccatc                                                                                                                                                                                                          >  1:1113314/1‑220 (MQ=255)
                                      g‑CTAAAAATTGCCGAGTTGCTTAATGGCGAGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAAt                                                                                                                                                                                                                                    <  1:492845/208‑1 (MQ=255)
                                                   cggtgtGCTTATGGGCGGGCCGGAGCCGGTGGGTGGGGAGGCGGAAGGGATAACGCCTGCGGGGATGGGGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAAc                                                                                                                                                                                         <  2:117544/234‑1 (MQ=255)
                                                       aggCTTAATGGCGGGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGt                                                                                                                                                                                     <  1:1199852/237‑1 (MQ=255)
                                                          tttAAGGTCGGGCCGCAGCCGCTGGATGGGGAAGGGCAAGGGAAAACGCCGGGGGGGATGGGGGCGGTTATTGAAGAAAATAACTGTATTGGCTGCAATAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGc                                                                                                                                                                                  <  2:1053139/238‑1 (MQ=255)
                                                                    agccgcagccgcTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGGATTCAGGCGTGTCCGGGAGACCCCATCGGTGGCGCTCCCCATGCCATTCATACGGTAAATGGGGATCTCTTTCCGGGCTGCAATTTATGTGTTTGTCACTGCCAGACCACCCGAATCTCGTTGCAACCGGGCGCAAAAaca                                                                                                                                                                           >  2:1199852/1‑236 (MQ=255)
                                                                     gccgcagacgcTGGATGGCGAAGCGAAAGAGATAACGGCTGCGCGGATGGTGGTGGTTATTGATGAAAATAACTGTTTTGGCTGCAATAAATGTAGTCCGGCGTGTCCGGTAGACGCCCTCGTTGGCGATACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTAAGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAATCACCt                                                                                                                                                                        <  2:400386/238‑1 (MQ=255)
                                                                             cgcTGGATGGCGAAACGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAAGGTATTCAGGCGGGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTcc                                                                                                                                                                  <  2:1297133/236‑1 (MQ=255)
                                                                                  gATGGGGAATCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAAt                                                                                                                                                           <  2:929960/238‑1 (MQ=255)
                                                                                          aGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTg                                                                                                                                                   <  1:894092/238‑1 (MQ=255)
                                                                                            cgcAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGGAATGAGTGATCTCTTTACGGGCTGCCATTTATTGGTTGTTCCTGTCCCCACGCACTCAATCTCGTTGCAACCGGGCGCCGAAACACCTGACCCCCGGAAATGGGactggcaa                                                                                                                                                >  2:893394/1‑232 (MQ=255)
                                                                                             gcAAGAGATAACGCCTGGGGGGTTGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTAGTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTTCCATGCATACGGTAATGAGTGATCTCTGTAAGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGaa                                                                                                                                                 <  1:627525/237‑1 (MQ=255)
                                                                                                     tAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTcc                                                                                                                                        <  1:246943/238‑1 (MQ=255)
                                                                                                        ccccTGCGCGGTTTGTTGGTTTTTTTGTTGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTTCCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTg                                                                                                                                    <  2:584606/237‑1 (MQ=255)
                                                                                                             gtgtgtgtggtggCGGTTTTTGATGAAAATAAGTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTGGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTAt                                                                                                                               <  2:44873/232‑1 (MQ=255)
                                                                                                                        tgggttttttGATGAAAATAACTTTTTTGGCTGGCCTAAATGTTTTGAGGCGTGTTCGGTAGACGCCCTCGGTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTg                                                                                                                     <  2:745665/230‑1 (MQ=255)
                                                                                                                             ttttttagtaaaagtaCTGTTTTGGCTGCCCTAAATGTATTCAGGCGTGTCCGGTAGGCGCCATGGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAaca                                                                                                                <  2:856228/223‑1 (MQ=255)
                                                                                                                             ttactggttaaaaTGACTGTCATGGCCGCACTAAATTTATTCAGGCGTGTCCGGTAAGCGCCATCGGTGGCGCTACCCGTGGCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAaca                                                                                                                <  1:270953/229‑1 (MQ=255)
                                                                                                                              ttttGATGAAAATAACTGTGTTGGCTGCACTAAATGGATTCAGGGGTGTCCGGTTGACGCCATCGGTGGCGCTACCCGTGCCATTGATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACAcc                                                                                                              <  2:461951/237‑1 (MQ=255)
                                                                                                                                 tgatgaAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTccc                                                                                                                        >  1:1005842/1‑226 (MQ=255)
                                                                                                                                 tgatgaAAAAAACTGTATTGGCTGCCCTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTccc                                                                                                                        <  2:1005842/226‑1 (MQ=255)
                                                                                                                                  tatgaAAATTATTGTTTTTGTTGGCCTAAATGTATCCCGGTGGTGCCGGTAGGCGCCATCGTTGGCCGTACCCTTGCCATGCATACGGTAATGAGGTATCTCTGGACGGGCTTCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑ttc                                                                                                         <  2:405185/238‑3 (MQ=255)
                                                                                                                                     gaAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGAACCCTGGAAATGGGATCTGAACACCATTCCCGGGCGGTACATTCCCGTGGAacacacagctta                                                                                                       >  2:865718/1‑231 (MQ=255)
                                                                                                                                             cTGTACTGGCTTCCCTCAATGTATTCAGGGTTGTCCGGTAGACGCCATCGTTGGCGCTATCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCGCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAGATGGGATCTGAACACCATTCCCGTCCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATtc                                                                                              <  1:497421/239‑1 (MQ=255)
                                                                                                                                                  ttGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGc                                                                                         >  2:947143/1‑239 (MQ=255)
                                                                                                                                                     gCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGCCTCCTGGCAATTGGATCTGAACACCATTCCCGTGCGGATCATTCCCGTGGCACCACA‑TGCTTAATTTTTTTACAGCAtt                                                                                      >  2:123356/1‑239 (MQ=255)
                                                                                                                                                       tGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCTTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAAAACCCTTCCCGTGCGTTCCATTCCCCGGGAACACCA‑TCCTTAAATTATTTCCTGCCTTcc                                                                                    >  2:570517/1‑238 (MQ=255)
                                                                                                                                                          aCTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTc                                                                                                                                         >  2:288992/1‑184 (MQ=255)
                                                                                                                                                          aCTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTc                                                                                                                                         <  1:288992/184‑1 (MQ=255)
                                                                                                                                                          aCTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTAGCATTCCCGTGGAACACCA‑TTCTTAAGTTATTCTCTGCATTCAGaa                                                                                 >  1:784338/1‑239 (MQ=255)
                                                                                                                                                          aCTAAAAGTTTTCAGGCCTGGCCGGTAGACGCCATCGTTTGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGa                                                                                  <  1:1305046/238‑1 (MQ=255)
                                                                                                                                                               aTGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGCACATGGATTCTGAACACCATTCCCGAGCGTATCTTTCCCGTGGAACACAC‑ACCTTAAGTTATTTTCTCCTTTCaaaaaaaaa                                                                            >  2:1128855/1‑238 (MQ=255)
                                                                                                                                                                      cAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGACCACCATTCCCGGGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTATGAATTCCGAAaaaataaaa                                                                        >  2:986705/1‑236 (MQ=255)
                                                                                                                                                                      cAGGCGTGGTCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCCTACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACCCCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCt                                                                     <  1:49687/239‑1 (MQ=255)
                                                                                                                                                                            tgtCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTAAGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGAtt                                                               >  1:123922/1‑239 (MQ=255)
                                                                                                                                                                            tgtCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCGGTAAGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTCCATCTCGTTGCAACCGGTCGCAGAACCACCTGACTCCTGGAACTGGGATCTGACCCCCATTCCCGTGCGTATCTTTCCCGTGGAACACAC‑TCCCTAAGTCATACTCTGCATTCAGAaaaaaaaaaaaaggggat                                                                >  2:608956/1‑230 (MQ=255)
                                                                                                                                                                              tCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGAACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAAAACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGCTTTc                                                             >  1:486360/1‑239 (MQ=255)
                                                                                                                                                                                   tAGACGACATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTGGCAGAAACACCTGACTCCTGAAACTGGGAACTGAAAACCATTCCCGTGCGTATCATTCCCGTGGAAAACCA‑TGCTTAAGTTATTCTTTGCATTACGAAAAAATAAAATCAGGGATTTAAAcgg                                                        >  2:677456/1‑239 (MQ=255)
                                                                                                                                                                                       cGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAAcggcgg                                                     <  1:865718/238‑1 (MQ=255)
                                                                                                                                                                                       cGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTAAGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACggcggc                                                    >  1:1208202/1‑239 (MQ=255)
                                                                                                                                                                                        gCCGTCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAAATTAGGTGTTGAACCGTGCCCAACGCACTGCATCGCGTTGCAACCGCTAGCAGAAACACCTCACTCCTGGAATGGGGTTCTGAACACACTTCCCGTGCCAATCATTACCGTGAAAAACCC‑ACTTTAAGTTATTCTCTGCCGTCGGAAAAAAAAAAATCtagaagttaaactcgggca                                                   >  2:1246586/1‑221 (MQ=255)
                                                                                                                                                                                        gCCATCGTTGGAGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGGGCCCGACGCACGGCATCTCGTTGCAACCGGTCGAAGAAACACCTGACTCCTGGAAATGGGAACTGGACACCATTACCTTGCGTGTCCTTCCCGTGGAA‑ACCACTGCTTGAGTTATTATCTGCATTCAGAAAAAATAACATCTGGGATTTCacaggaggc                                                    >  2:1028373/1‑230 (MQ=255)
                                                                                                                                                                                                tGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCAtccatcca                                           <  1:947143/239‑1 (MQ=255)
                                                                                                                                                                                                 ggCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCAc                                          <  1:662064/239‑1 (MQ=255)
                                                                                                                                                                                                    ctccaccggtcccttcATACGGGAATTGGTGATTTCTGGAAGGGCCTCAATTTATGGGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCg                                        <  2:1208202/223‑1 (MQ=255)
                                                                                                                                                                                                         ccGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAATCGCCTGGCTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCAATCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATg                                  <  1:422077/239‑1 (MQ=255)
                                                                                                                                                                                                                aTGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAAcc                            <  1:123356/238‑1 (MQ=255)
                                                                                                                                                                                                                    atttgtgaaTGAGGTATCCCTGTACGGGCTGCAATTTATGTGTTGACCCGTGCCCGACGCACTGAATTTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTCTTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAg                         <  1:1246586/230‑1 (MQ=255)
                                                                                                                                                                                                                      aCGGTAATGAGTGAGCTCTGTACGGGCTGCAATTTATGTGTTCATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACCCCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCa                     <  1:1028373/239‑1 (MQ=255)
                                                                                                                                                                                                                       cGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCaa                    <  1:570517/239‑1 (MQ=255)
                                                                                                                                                                                                                       cGGTAATGAGTGATCTCTGGACGGGCTTCAATTTAAGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCaa                    <  1:1128855/239‑1 (MQ=255)
                                                                                                                                                                                                                       cGGTAATGAGGGATCTCTGTACGGGCTTCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCaa                    <  1:986705/239‑1 (MQ=255)
                                                                                                                                                                                                                             tGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAAcggc                                                       >  1:845455/1‑198 (MQ=255)
                                                                                                                                                                                                                             tGAGGGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAAcggc                                                       <  2:845455/198‑1 (MQ=255)
                                                                                                                                                                                                                             tGAATGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATGCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAACAAATAAAATCTGGGATTTCAACGGCGGACTCCCTCCACCGGAGATGAAAACCCAGGCCAACggga               >  1:842833/1‑236 (MQ=255)
                                                                                                                                                                                                                               aGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGTAGATGAAAACCAAGTCCAACGGTaca             >  1:605377/1‑238 (MQ=255)
                                                                                                                                                                                                                                    tctcTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATTTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGCATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAATTAAACTATGGGATTCCAACAAGGGCACCCTCACCCCGGGGATTAAAACCCAACCCAAAggacaaccccc        >  2:1155189/1‑229 (MQ=255)
                                                                                                                                                                                                                                        tGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGGACACCCCTgccc    >  1:260953/1‑236 (MQ=255)
                                                                                                                                                                                                                                        tGTAAGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAACATAAAATCTGGGATTTCAAAGGGGGCATTCATCAACCGGAAATGAAAACCCAGTCCAAAGGGAAACCCCTgggcg   >  2:1363051/1‑238 (MQ=255)
                                                                                                                                                                                                                                         gTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTAAACGGCGGCATCCATCCACCGGAAGTGAAAACCCAATCCAACGGAACAACCCTGCCCCa  >  2:197785/1‑239 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                            
AAAATCAACCGTTGCGCCCCAGGTGGCGAAGCTGTGATG‑CTAAAAATTGCCGAGTTGCTTAATGTCGAGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCA‑TGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCA  >  NC_000913/1706555‑1707025

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