Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,354 C→T L59L (CTC→CTT 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,3540CT98.3% 189.3 / ‑4.4 59L59L (CTC→CTTcreDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (31/27);  minor base G (0/1);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

CACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTT  >  NC_000913/4638114‑4638590
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                            
cacaTAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTc                                                                                                                                                                                                                                    >  1:82772/1‑251 (MQ=255)
    tAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAACCTCGTTGGGCCGCTTATCGCTATTCCtgt                                                                                                                                                                                                                                >  2:107097/1‑251 (MQ=255)
       cTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGAc                                                                                                                                                                                                                             >  2:126207/1‑251 (MQ=255)
         tCAAATTCTTCCCACATAGCCTTCTTATCCTGCGGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGAcc                                                                                                                                                                                                                            <  2:200301/250‑1 (MQ=255)
          cAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGa                                                                                                                                                                                                                          <  1:53258/251‑1 (MQ=255)
           aaaTTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCg                                                                                                                                                                                                                           >  1:448762/1‑249 (MQ=255)
            aaTTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGc                                                                                                                                                                                                                        <  2:189622/251‑1 (MQ=255)
              ttcttcCCACATAGCCTTCGTATCCTGCTGACATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTACGCTAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTCTCGCTATTCCTGTGACCGAGCtt                                                                                                                                                                                                                      >  1:538537/1‑251 (MQ=255)
                   cccACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTtata                                                                                                                                                                                                                  <  1:220783/250‑1 (MQ=255)
                        aTAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGTGCCGCTTATCGCTAt                                                                                                                                                                                                                                      <  2:519944/225‑1 (MQ=255)
                        aTAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTAt                                                                                                                                                                                                                                      >  1:519944/1‑225 (MQ=255)
                          atccTTCGTATCCTGCTGCCATTGCTAAGGAGAAGACTAGGTTGTACTCCCCCTTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATGGCTGTTGTTGATTCCGATAATGCTTATTCGGCAAGTGATTGGCGACCGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCa                                                                                                                                                                                                           <  1:508/247‑1 (MQ=255)
                           gccTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAgg                                                                                                                                                                                                         <  1:88778/249‑1 (MQ=255)
                                 gTATCCTGCTGCCATTGCAAAGGCGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAgg                                                                                                                                                                                                   >  2:468220/1‑251 (MQ=255)
                                 gTATACTGCTGCCAGTGCAAAGGAGAAGACTATGGTGAAATCCCCCCTGTTCTGGAAAATGACCAACTTGTTTGGTGCAGTAGTGCTGTTGGTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGGAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTGATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAgg                                                                                                                                                                                                   <  1:533357/251‑1 (MQ=255)
                                   aTCCTGCTGCCATTGCACAGGCGAAGACTATGTTGAAATCCCCACTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCAGATAATGCTGCTTCGGCAAGTGATTGTCGAACGTGCCGGTTACCGTAGCGAAGTGGAAGATGCGATTCGCCAAAGGACCAGCGGGCCGCCAAAAGTCGTTGGGCCGCTTATCGCGATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGAt                                                                                                                                                                                                 >  1:205328/1‑251 (MQ=255)
                                   aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGGGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGa                                                                                                                                                                                                  <  2:246951/250‑1 (MQ=255)
                                     ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATTTGGAAGATGCGATGCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTAGTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATaa                                                                                                                                                                                               <  1:380060/251‑1 (MQ=255)
                                                 gCAAAGGAGAAGACTATGTTGAAAGCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCg                                                                                                                                                                                   <  1:176506/251‑1 (MQ=255)
                                                        agaagaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAgg                                                                                                                                                                                                   <  1:82821/228‑1 (MQ=255)
                                                        agaagaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAgg                                                                                                                                                                                                   >  2:82821/1‑228 (MQ=255)
                                                        agaagaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACgg                                                                                                                                                                            >  2:499610/1‑250 (MQ=255)
                                                                  gTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATcc                                                                                                                                                                  >  1:450811/1‑251 (MQ=255)
                                                                     gAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCAtt                                                                                                                                                               <  1:499610/251‑1 (MQ=255)
                                                                     gAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTGGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATCACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCGGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCAt                                                                                                                                                                <  1:509522/250‑1 (MQ=255)
                                                                          ccccccTGTTCTGGAAAATGCCCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCAGCTTCGGCCAGTGATTGTCGCACGTGCCGATTACCGCAGCGATGGGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAACACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGTAGTTTTTTCCATTTTTgg                                                                                                                                                          >  1:229981/1‑251 (MQ=255)
                                                                               cTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTAcc                                                                                                                                                     >  2:142449/1‑251 (MQ=255)
                                                                                  ttCTGGAAACTGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCGAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGCGGAAGATGCTATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTg                                                                                                                                                   >  2:511508/1‑250 (MQ=255)
                                                                                           atgaccaccttGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAttgat                                                                                                                                         <  2:13155/241‑1 (MQ=255)
                                                                                              accaccttGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAttgatggt                                                                                                                                      >  1:179128/8‑251 (MQ=255)
                                                                                               ccaccttGTTTGGGGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGGAGCGATGTGGAAGAGGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAttgatggtt                                                                                                                                     <  2:516058/245‑1 (MQ=255)
                                                                                                caccttGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAttgatggttg                                                                                                                                    >  1:170958/6‑251 (MQ=255)
                                                                                                     tGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGCTTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGc                                                                                                                               >  1:409751/1‑251 (MQ=255)
                                                                                                      gTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGc                                                                                                                               <  2:218688/250‑1 (MQ=255)
                                                                                                        ttGGTGCAGTAGTGCTGTTGTTGCGTCCGATAATGCTGAATCGGCAAGTGATTGGCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGGACCAGCGGGCCGCAAAAACGCGTTGGGCCGCTTATCGCTATTCGTGTGACCGAGCTTTATACGGGTCAGGAATAGGAGAAATCCGTGGAGCGGAAACGAATTTTTATCCATTTTTGGTTACCTGTGTCAttgaagttggacggcaa                                                                                                                             <  2:289193/250‑14 (MQ=255)
                                                                                                           gTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCagtc                                                                                                                           >  1:113067/1‑246 (MQ=255)
                                                                                                                gTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGt                                                                                                                    >  1:345456/1‑251 (MQ=255)
                                                                                                                  tttGCTGTAGTTGATTCCGATAATGCTGAGTCGGCAAGTGATTGTCGAACGTGCCGATGACCGTAGCGATGTTGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAGCTCGTTTGGCCGCTTATCGCTATTCCTGTGACCGTGCTGTAGACGGTGCAGGAAGAGGATAAAACCGTGGAGCGTAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTg                                                                                                                   <  2:458514/249‑1 (MQ=255)
                                                                                                                    tGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGga                                                                                                                 <  1:173608/250‑1 (MQ=255)
                                                                                                                              gtatcctataatgatgtgtCGGCAAGTTATTGACGAACGTGCCGAGTACCGTAGCGGTGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGTAGGAAGAGGATAAAGCCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGGCATTGATGGTTGATGGCAAGCAGAACGTGTAAGAACGCAAGa                                                                                                      <  2:123054/233‑1 (MQ=255)
                                                                                                                              gATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAg                                                                                                       <  1:525646/250‑1 (MQ=255)
                                                                                                                                 tCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGAt                                                                                                     >  1:535000/1‑249 (MQ=255)
                                                                                                                                  ccGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGGGGAAGAACGCAAGATAgg                                                                                                  <  1:239448/251‑1 (MQ=255)
                                                                                                                                    gATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAggg                                                                                                 >  1:18549/1‑250 (MQ=255)
                                                                                                                                        aTGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTta                                                                                            <  1:152267/251‑1 (MQ=255)
                                                                                                                                                       ctGATTGTCGCACGTGCCGATTACCGAAGCGATGAGGAAGAAGCGATTCGTCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATTGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTg                                                                             >  2:180379/2‑251 (MQ=255)
                                                                                                                                                         gATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATCCGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGc                                                                           >  2:263800/1‑251 (MQ=255)
                                                                                                                                                                 aaCGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCCGGAAGAGGATAAAACCGTGGAGCGGAAACGCAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGAt                                                                   >  2:507680/1‑251 (MQ=255)
                                                                                                                                                                      gccgATTACCGTAGCGATGTGGAAGATGCGAGTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGGCTGGCACAGTGATTTAAc                                                              <  2:179128/248‑1 (MQ=255)
                                                                                                                                                                         gATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGGACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACgct                                                           <  2:525228/251‑3 (MQ=255)
                                                                                                                                                                                 tAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCCAGATAGGGATTTATACCGGTCAGGTCTGGCACCGTGATTTAACGCTAAAAGCCg                                                   >  1:582033/1‑251 (MQ=255)
                                                                                                                                                                                    cGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGAtt                                                >  2:308262/1‑251 (MQ=255)
                                                                                                                                                                                      aTGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGCCTGGCACAGTGATTTAACGCTAAAAGCCGATTTc                                              >  1:135426/1‑251 (MQ=255)
                                                                                                                                                                                           gAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGt                                         <  1:315329/251‑1 (MQ=255)
                                                                                                                                                                                                   aaTTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGCGCTTTATACGGAGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCAGATTTCGATGTTTCGCGTc                                 >  2:174587/2‑251 (MQ=255)
                                                                                                                                                                                                     ttCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCtt                               >  1:587707/1‑251 (MQ=255)
                                                                                                                                                                                                             aGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGAGGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATTCCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGagct                       <  1:76957/251‑4 (MQ=255)
                                                                                                                                                                                                                               aaaCTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCat     >  1:196988/1‑250 (MQ=255)
                                                                                                                                                                                                                                  cTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCtt  >  1:177790/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                            
CACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTT  >  NC_000913/4638114‑4638590

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