Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,706,790 C→A 9.9% 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,7900CA 9.9% 168.0 / 13.1 62T148K (ACG→AAG) rsxBSoxR iron‑sulfur cluster reduction factor component; putative iron‑sulfur protein
Reads supporting (aligned to +/- strand):  ref base C (20/36);  new base A (4/2);  total (24/38)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-01

AAAATCAACCGTTGCGCCCCAGGTGGCGAAGCTGTGATGCTAAAAATTGCCGAGTTGCTTAATGTCGAGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCA  >  NC_000913/1706555‑1707025
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                            
aaaaCCACCCGTTGCGCCCAAGGGGGGGAAGGTGGGATGTTAAAAATGCCCGAGTTGGTTAAGGGCGGGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACggg                                                                                                                                                                                                                                           <  2:1003536/239‑1 (MQ=255)
    tCACCCGTTGCCCCCCAGGTGGCGAAGCTGTGATGCTAAAAATTGCCGAGTTGCTTAATGTCGAGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGc                                                                                                                                                                                                                                       <  1:994761/239‑1 (MQ=255)
    ccccccGTTTCCCCCCGGGTGGGGAAGCTTGGTGTGTAAAAATTGCCGGGTTGCTTAATGTGGAGCCGCCGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGc                                                                                                                                                                                                                                       <  2:724693/235‑1 (MQ=255)
         cgtgtcgcCCCAGGTGGCGAAGCTGTGATGCTAAAAATGTCCGAGTTGCTTAATGTCGAGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAAttt                                                                                                                                                                                                                                  <  1:1313446/234‑1 (MQ=255)
               gcCCCAGGTGGCGAAGCTGTGATGCTAAAAATTGCCGAGTTGCTTAATGTCGAGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATgtgt                                                                                                                                                                                                                            <  1:381476/239‑1 (MQ=255)
                       tGGCGAAGCTGTGATGCTAAAAATTGCCGAGTTGCTTAATGTCGAGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATcc                                                                                                                                                                                                                      <  1:1172238/237‑1 (MQ=255)
                                      gCTAAAAATTGCCGAGTTGCTTAATGGCGAGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAAt                                                                                                                                                                                                                                    <  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)
                                                                                                                             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)
                                                                                                                                  tatgaAAATTATTGTTTTTGTTGGCCTAAATGTATCCCGGTGGTGCCGGTAGGCGCCATCGTTGGCCGTACCCTTGCCATGCATACGGTAATGAGGTATCTCTGGACGGGCTTCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCAttc                                                                                                         <  2:405185/238‑3 (MQ=255)
                                                                                                                                     gaAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGAACCCTGGAAATGGGATCTGAACACCATTCCCGGGCGGTACATTCCCGTGGAacacacagctta                                                                                                       >  2:865718/1‑231 (MQ=255)
                                                                                                                                             cTGTACTGGCTTCCCTCAATGTATTCAGGGTTGTCCGGTAGACGCCATCGTTGGCGCTATCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCGCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAGATGGGATCTGAACACCATTCCCGTCCGTATCATTCCCGTGGAACACCATGCTTAAGTTATtc                                                                                              <  1:497421/239‑1 (MQ=255)
                                                                                                                                                  ttGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGc                                                                                         >  2:947143/1‑239 (MQ=255)
                                                                                                                                                     gCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGCCTCCTGGCAATTGGATCTGAACACCATTCCCGTGCGGATCATTCCCGTGGCACCACATGCTTAATTTTTTTACAGCAtt                                                                                      >  2:123356/1‑239 (MQ=255)
                                                                                                                                                       tGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCTTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAAAACCCTTCCCGTGCGTTCCATTCCCCGGGAACACCATCCTTAAATTATTTCCTGCCTTcc                                                                                    >  2:570517/1‑238 (MQ=255)
                                                                                                                                                          aCTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTc                                                                                                                                         >  2:288992/1‑184 (MQ=255)
                                                                                                                                                          aCTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTc                                                                                                                                         <  1:288992/184‑1 (MQ=255)
                                                                                                                                                          aCTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTAGCATTCCCGTGGAACACCATTCTTAAGTTATTCTCTGCATTCAGaa                                                                                 >  1:784338/1‑239 (MQ=255)
                                                                                                                                                          aCTAAAAGTTTTCAGGCCTGGCCGGTAGACGCCATCGTTTGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGa                                                                                  <  1:1305046/238‑1 (MQ=255)
                                                                                                                                                               aTGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGCACATGGATTCTGAACACCATTCCCGAGCGTATCTTTCCCGTGGAACACACACCTTAAGTTATTTTCTCCTTTCaaaaaaaaa                                                                            >  2:1128855/1‑238 (MQ=255)
                                                                                                                                                                      cAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGACCACCATTCCCGGGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTATGAATTCCGAAaaaataaaa                                                                        >  2:986705/1‑236 (MQ=255)
                                                                                                                                                                      cAGGCGTGGTCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCCTACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACCCCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCt                                                                     <  1:49687/239‑1 (MQ=255)
                                                                                                                                                                            tgtCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTAAGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGAtt                                                               >  1:123922/1‑239 (MQ=255)
                                                                                                                                                                            tgtCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCGGTAAGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTCCATCTCGTTGCAACCGGTCGCAGAACCACCTGACTCCTGGAACTGGGATCTGACCCCCATTCCCGTGCGTATCTTTCCCGTGGAACACACTCCCTAAGTCATACTCTGCATTCAGAaaaaaaaaaaaaggggat                                                                >  2:608956/1‑230 (MQ=255)
                                                                                                                                                                              tCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGAACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAAAACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGCTTTc                                                             >  1:486360/1‑239 (MQ=255)
                                                                                                                                                                                   tAGACGACATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTGGCAGAAACACCTGACTCCTGAAACTGGGAACTGAAAACCATTCCCGTGCGTATCATTCCCGTGGAAAACCATGCTTAAGTTATTCTTTGCATTACGAAAAAATAAAATCAGGGATTTAAAcgg                                                        >  2:677456/1‑239 (MQ=255)
                                                                                                                                                                                       cGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAAcggcgg                                                     <  1:865718/238‑1 (MQ=255)
                                                                                                                                                                                       cGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTAAGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACggcggc                                                    >  1:1208202/1‑239 (MQ=255)
                                                                                                                                                                                        gCCATCGTTGGAGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGGGCCCGACGCACGGCATCTCGTTGCAACCGGTCGAAGAAACACCTGACTCCTGGAAATGGGAACTGGACACCATTACCTTGCGTGTCCTTCCCGTGGAAACCACTGCTTGAGTTATTATCTGCATTCAGAAAAAATAACATCTGGGATTTCacaggaggc                                                    >  2:1028373/1‑230 (MQ=255)
                                                                                                                                                                                                tGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCAtccatcca                                           <  1:947143/239‑1 (MQ=255)
                                                                                                                                                                                                 ggCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCAc                                          <  1:662064/239‑1 (MQ=255)
                                                                                                                                                                                                         ccGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAATCGCCTGGCTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCAATCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATg                                  <  1:422077/239‑1 (MQ=255)
                                                                                                                                                                                                                aTGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAAcc                            <  1:123356/238‑1 (MQ=255)
                                                                                                                                                                                                                    atttgtgaaTGAGGTATCCCTGTACGGGCTGCAATTTATGTGTTGACCCGTGCCCGACGCACTGAATTTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTCTTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAg                         <  1:1246586/230‑1 (MQ=255)
                                                                                                                                                                                                                      aCGGTAATGAGTGAGCTCTGTACGGGCTGCAATTTATGTGTTCATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACCCCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCa                     <  1:1028373/239‑1 (MQ=255)
                                                                                                                                                                                                                       cGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCaa                    <  1:570517/239‑1 (MQ=255)
                                                                                                                                                                                                                       cGGTAATGAGTGATCTCTGGACGGGCTTCAATTTAAGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCaa                    <  1:1128855/239‑1 (MQ=255)
                                                                                                                                                                                                                       cGGTAATGAGGGATCTCTGTACGGGCTTCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCaa                    <  1:986705/239‑1 (MQ=255)
                                                                                                                                                                                                                             tGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAAcggc                                                       >  1:845455/1‑198 (MQ=255)
                                                                                                                                                                                                                             tGAGGGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAAcggc                                                       <  2:845455/198‑1 (MQ=255)
                                                                                                                                                                                                                             tGAATGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATGCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAACAAATAAAATCTGGGATTTCAACGGCGGACTCCCTCCACCGGAGATGAAAACCCAGGCCAACggga               >  1:842833/1‑236 (MQ=255)
                                                                                                                                                                                                                               aGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGTAGATGAAAACCAAGTCCAACGGTaca             >  1:605377/1‑238 (MQ=255)
                                                                                                                                                                                                                                    tctcTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATTTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGCATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAATTAAACTATGGGATTCCAACAAGGGCACCCTCACCCCGGGGATTAAAACCCAACCCAAAggacaaccccc        >  2:1155189/1‑229 (MQ=255)
                                                                                                                                                                                                                                        tGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGGACACCCCTgccc    >  1:260953/1‑236 (MQ=255)
                                                                                                                                                                                                                                        tGTAAGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAACATAAAATCTGGGATTTCAAAGGGGGCATTCATCAACCGGAAATGAAAACCCAGTCCAAAGGGAAACCCCTgggc  >  2:1363051/1‑238 (MQ=255)
                                                                                                                                                                                                                                         gTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTAAACGGCGGCATCCATCCACCGGAAGTGAAAACCCAATCCAACGGAACAACCCTGCCCCa  >  2:197785/1‑239 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                            
AAAATCAACCGTTGCGCCCCAGGTGGCGAAGCTGTGATGCTAAAAATTGCCGAGTTGCTTAATGTCGAGCCGCAGCCGCTGGATGGCGAAGCGCAAGAGATAACGCCTGCGCGGATGGTGGCGGTTATTGATGAAAATAACTGTATTGGCTGCACTAAATGTATTCAGGCGTGTCCGGTAGACGCCATCGTTGGCGCTACCCGTGCCATGCATACGGTAATGAGTGATCTCTGTACGGGCTGCAATTTATGTGTTGATCCGTGCCCGACGCACTGCATCTCGTTGCAACCGGTCGCAGAAACACCTGACTCCTGGAAATGGGATCTGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCA  >  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: