Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,264 G→A Q29Q (CAG→CAA creD → putative inner membrane protein CreD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,638,2640GA97.9% 145.5 / NA 49Q29Q (CAG→CAAcreDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (20/27);  minor base C (1/0);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGG  >  NC_000913/4638021‑4638506
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
cAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTTGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGAtt                                                                                                                                                                                                                                               <  1:329964/250‑1 (MQ=255)
     aGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGaa                                                                                                                                                                                                                                         <  1:314724/251‑1 (MQ=255)
            ccGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCCCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGccg                                                                                                                                                                                                                                  <  2:183997/251‑3 (MQ=255)
                  gtttAACGGCGAAGTCTCGCGGCGCAACGTGCAGGAATGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTTCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATGCGGCAAGGGATTGTCGAACGTGCCTATTAcc                                                                                                                                                                                                                            <  2:56344/251‑1 (MQ=255)
                  gtttAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTAc                                                                                                                                                                                                                             <  1:508561/250‑1 (MQ=255)
                         ggCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCg                                                                                                                                                                                                                      >  1:103262/1‑250 (MQ=255)
                             aaGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGCCTTCACCGTCACGTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTATTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGAtgtg                                                                                                                                                                                                                 <  1:586035/251‑1 (MQ=255)
                             aaGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGAtgtg                                                                                                                                                                                                                 <  2:186844/251‑1 (MQ=255)
                                          cAACGTGCATGAAGGTGTCTTGCTGTCCTCGCTGTGACTTCACCTTCACTTCACATAGCTGCAAATTCTTACCACCTAGCCTTCGTATTCTGTTGCCATTGCAAAGGAGAAGACTATGTTGAAAGCCCCCCTGTTCTGGAAAATGACCACCTTGGTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTc                                                                                                                                                                                                    <  2:587465/251‑1 (MQ=255)
                                                             gTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGcc                                                                                                                                                                                 >  2:315728/1‑251 (MQ=255)
                                                              tGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCg                                                                                                                                                                                >  1:294568/1‑251 (MQ=255)
                                                                        gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCg                                                                                                                                                                      <  1:587711/251‑1 (MQ=255)
                                                                         cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGCTTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGCAGATGCGATTCGCCAAAGTACCAGCGGGCCGAAAAAACTCGt                                                                                                                                                                     >  2:631852/1‑251 (MQ=255)
                                                                                  cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGTCCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCt                                                                                                                                                            <  1:558172/251‑1 (MQ=255)
                                                                                                  aGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCtgtg                                                                                                                                            <  1:598509/251‑1 (MQ=255)
                                                                                                      tCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTAc                                                                                                                                                                                          <  2:220687/201‑1 (MQ=255)
                                                                                                      tCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTAc                                                                                                                                                                                          >  1:220687/1‑201 (MQ=255)
                                                                                                        aaaTTCTTCCCACATATCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATTTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAATCTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAg                                                                                                                                      <  2:418704/251‑1 (MQ=255)
                                                                                                            tcttcCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCtt                                                                                                                                   <  1:504562/250‑1 (MQ=255)
                                                                                                                     aTAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGc                                                                                                                         >  1:164429/1‑251 (MQ=255)
                                                                                                                     aTAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTGGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGTAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGc                                                                                                                         <  1:84071/251‑1 (MQ=255)
                                                                                                                                  ccTGCTGACATTGCAAAGGAGAAGTCTATGTTGAAATCCTCCCTGTTCTGGTACATGTCCCCCTTGTTTGGTGCAGTATTGCTGTTGTTTATTCCGATAAAGCTGACTCGGCAAGGGATTGTCGAACGTGGCGATTATCTTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATTCGGTGTAGGAAGAGGATa                                                                                                             <  2:524426/250‑1 (MQ=255)
                                                                                                                                        gcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAAcc                                                                                                        <  1:315728/249‑1 (MQ=255)
                                                                                                                                           atagCAAAGGAGAATATTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTTGTGCAGTATTGCTGTTGGTGATTCCGATAATGTTGATTCGGCAAGTGATTGGCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAATCCGTGGa                                                                                                   <  2:123975/248‑1 (MQ=255)
                                                                                                                                                   ggAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAAc                                                                                           >  1:420196/1‑251 (MQ=255)
                                                                                                                                                   ggAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAAc                                                                                           <  1:45422/251‑1 (MQ=255)
                                                                                                                                                     agaagaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAACACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAACCCGTGGAg                                                                                                  >  2:58300/1‑242 (MQ=255)
                                                                                                                                                     agaagaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAg                                                                                                  <  1:58300/242‑1 (MQ=255)
                                                                                                                                                       aagaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGa                                                                                        >  1:232798/1‑249 (MQ=255)
                                                                                                                                                        agaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGt                                                                                     <  1:220653/251‑1 (MQ=255)
                                                                                                                                                                     aTCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCAttttt                                                                        <  1:527712/251‑1 (MQ=255)
                                                                                                                                                                          cccgtttcgggaaaatgaccccattGTTTGGTGTAGTATTGATGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGATCGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGTGGGCCGCAAAAACTCGTTGGTCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATTCATTTTTGGTTa                                                                   <  2:535481/227‑1 (MQ=255)
                                                                                                                                                                             tGTTCTGGAAAATGAACACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCt                                                                >  1:352183/1‑251 (MQ=255)
                                                                                                                                                                              gTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTg                                                               >  1:536789/1‑251 (MQ=255)
                                                                                                                                                                                  tGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGcc                                                                                                                                                                                 >  1:445684/1‑134 (MQ=255)
                                                                                                                                                                                  tGGAAAATGACCACCGTGTGTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATGCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCATCGGGcc                                                                                                                                                                                 <  2:445684/134‑1 (MQ=255)
                                                                                                                                                                                    gAAAATGACCACCTTGTTTGGAGCAGAATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCAAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGGGGAGCGGAAACGAGGTTTTATCCATTATTGGTTACCTggggcat                                                         >  2:555033/1‑245 (MQ=255)
                                                                                                                                                                                     aaaaTGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGACGTTTTATCCATTTTTGGTTACCTGAGTCAtt                                                        >  1:673996/1‑251 (MQ=255)
                                                                                                                                                                                                        gTGCAGTTTTGCTGTCGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGGGCCGATGACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACGCGTTGGGACGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGGTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAg                                     <  2:181873/251‑1 (MQ=255)
                                                                                                                                                                                                        gTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCACGTGATAGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGAGTCGCCAAAGTACCATCGGGCCGCAAAAACTCGATGTGCCACTTCTCGCTATTACTGTGACCGAGCCTTATACGGTGCAGGAAGAGGATAAAACCGAGGCGCGGACCGTAAGTTTTAACCATTCTTGGTTACCTGAGTCATTGATGGTTGCTGGCCATCa                                      >  2:191940/1‑250 (MQ=255)
                                                                                                                                                                                                                   cTGTTGTTGATTCCGATAATGATGATTCGGCAAGTGATTGTCGAACGTGCCGCTTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTAAGTCATTGACGGTTGATAGAAATCAGAACGTGga                             >  2:453262/1‑248 (MQ=255)
                                                                                                                                                                                                                         ttgATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTTCCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGaa                                                                                             <  2:34448/179‑1 (MQ=255)
                                                                                                                                                                                                                         ttgATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGaa                                                                                             >  1:34448/1‑179 (MQ=255)
                                                                                                                                                                                                                             ttCCGATAATGCTGATTCGGCCAGTGATTGTCGAACGTGCCGATCACCGTAGCGATGTGGAAGATGCGCTTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGATTATCGATATTCCTGTGACCGCGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAACTTTTATCCCTTTTTGGTTACCTGCGACATTGGTGGTTGATGGCATTCAGTACGTGGAAGAACGCAAGAt                 >  2:38523/1‑250 (MQ=255)
                                                                                                                                                                                                                              tCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAACGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGCGCTTTATACGGTGCAGGAAGAGGATAAAAGCGTGGAGCGGAAACGAAGTTTTATCCATTTGTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGACGAACGCaggaa                 >  2:563015/1‑245 (MQ=255)
                                                                                                                                                                                                                                    aaTGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGACGAACGCAAGATAGGGAttt         >  1:102153/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aaTGCTGATTCGGAAAGTGATTGTCGAACGTGCCGAGTACCGTAGAGATGGGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGAttt         <  2:547763/251‑1 (MQ=255)
                                                                                                                                                                                                                                      tGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTtat       >  1:201852/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aTTCGGCAAGTGATTGTCGAACGTGCAGATTACCGTAGCGATGTGGAAGATGCGATTCTCCAAAGTACCAGGGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGGGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGCTTTATTCCgg  <  2:562475/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
CAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGG  >  NC_000913/4638021‑4638506

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