Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,873,919 C→T 5.2% G25D (GGT→GAT)  cysC ← adenosine 5'‑phosphosulfate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,873,9190CT 5.2% 628.0 / 16.0 193G25D (GGT→GAT) cysCadenosine 5'‑phosphosulfate kinase
Reads supporting (aligned to +/- strand):  ref base C (91/91);  new base T (4/6);  total (96/97)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-02

Only 100 of 197 total aligned reads displayed.

GGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGATCGATAAAAATCAGCCCACCCGTCACCGGAT  >  NC_000913/2873647‑2874185
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                          
ggcgggcGGGGCGGTGGTGCGAGATACAAGGGGGAAGCACCGCCGGGCCGGCTTCAACCATCAAATTCGCCACTTCACCGGCGCGACGGATATTCCCTTTACGATCGGCATCGCTAAAACCGGGATCGCTGCATAGGCCGTGGCGAACATTGTCGCCATCCAGCAGATACGGACTGACGCCGAGGTTATGTAACGCCGCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGAcc                                                                                                                                                                                                                                                                           <  4:574511/268‑1 (MQ=255)
            ggTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTcg                                                                                                                                                                                                                                                              <  2:239309/275‑1 (MQ=255)
                        tAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCCAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTCTGTAACGCCTCCTCACGCGCCCCGGCGCCCGTTGATTTACCGGACCCGGAGAGCCCGCTACACCCCAGAACCACACCACGATCCACGCGTTGTACCCCCCgttgtttccc                                                                                                                                                                                                                                                   >  4:304387/1‑270 (MQ=255)
                                                    ttCAACCCTCAAATTCGCCACTTCCCCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGt                                                                                                                                                                                                                      <  2:298884/275‑1 (MQ=255)
                                                                                              ctctTTACGATCGGCATCGCTAAAACCGAGATCGCTTCATAATCCGTGGCGAACATTGTCGCCATCCATCAGATACGTACTGACGCCGAGGTTATGTAACGCCTCCTCCAGCGCCCCGGCGTCCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGGCCGTGGTGTAGCTcg                                                                                                                                                                                                                                                              <  1:190229/193‑1 (MQ=255)
                                                                                                                           gaTCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAg                                                                                                                                                                                                                                                                  >  2:200243/1‑160 (MQ=255)
                                                                                                                               gCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTcg                                                                                                                                                                                                                                                              <  4:200243/160‑1 (MQ=255)
                                                                                                                                cTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTACCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCCTTATTTATCCCCCAGCAAATCGCGCGCGCCCCACTTCCGCAAGTGCCGACGAACCaga                                                                                                                                          >  2:281527/1‑275 (MQ=255)
                                                                                                                                 tGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCg                                                                                                                                                                                                                                                                          >  4:83465/1‑146 (MQ=255)
                                                                                                                                    aTAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGgtgt                                                                                                                                                                                                                                                                    <  8:691879/149‑1 (MQ=255)
                                                                                                                                     tAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTa                                                                                                                                                                                                                                                                   <  7:440487/149‑1 (MQ=255)
                                                                                                                                     tAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTa                                                                                                                                                                                                                                                                   <  5:332156/149‑1 (MQ=255)
                                                                                                                                       aTCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTAgt                                                                                                                                                                                                                                                                 <  7:113769/149‑2 (MQ=255)
                                                                                                                                        tCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTTCCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCCGTGCGGAAAGTGGCGACGAACCAGAGCATTCa                                                                                                                                   >  2:460479/1‑274 (MQ=255)
                                                                                                                                         ccGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTc                                                                                                                                                                                                                                                               >  7:938401/1‑149 (MQ=255)
                                                                                                                                             ggCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCaa                                                                                                                                  >  4:274220/1‑270 (MQ=255)
                                                                                                                                               cGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttg                                                                                                                                                                                                                                                         <  5:225042/149‑1 (MQ=255)
                                                                                                                                               cGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTGTATGTAACTCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCgg                                                                                                                                                                                                                                          <  1:54591/164‑1 (MQ=255)
                                                                                                                                                gAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGgtgt                                                                                                                                                                                                                                                                    <  5:1062197/137‑1 (MQ=255)
                                                                                                                                                gAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGgtgt                                                                                                                                                                                                                                                                    >  6:1062197/1‑137 (MQ=255)
                                                                                                                                                gAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgt                                                                                                                                                                                                                                                        >  8:1045323/1‑149 (MQ=255)
                                                                                                                                                gAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgt                                                                                                                                                                                                                                                        <  6:507803/149‑1 (MQ=255)
                                                                                                                                                 aaCATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgtt                                                                                                                                                                                                                                                       <  6:343331/149‑1 (MQ=255)
                                                                                                                                                 aaCATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgtt                                                                                                                                                                                                                                                       >  5:803905/1‑149 (MQ=255)
                                                                                                                                                 aaCATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgtt                                                                                                                                                                                                                                                       <  7:1045323/149‑1 (MQ=255)
                                                                                                                                                    atgttCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgtt                                                                                                                                                                                                                                                       >  6:1071210/5‑146 (MQ=255)
                                                                                                                                                    aTTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGca                                                                                                                                                                                                                                                    >  7:98410/1‑149 (MQ=255)
                                                                                                                                                      tGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGtt                                                                                                                                                                                                                     >  2:707606/1‑178 (MQ=255)
                                                                                                                                                        tCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGt                                                                                                                                                                                                                                                <  7:291523/149‑1 (MQ=255)
                                                                                                                                                          gCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGa                                                                                                                                                                                                                                              >  5:86433/1‑149 (MQ=255)
                                                                                                                                                          gCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCg                                                                                                                                                                                                                 <  4:707606/178‑1 (MQ=255)
                                                                                                                                                              tCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAAt                                                                                                            >  2:427745/1‑275 (MQ=255)
                                                                                                                                                                cagcagATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGAt                                                                                                                                                                                                                                        <  7:54620/149‑1 (MQ=255)
                                                                                                                                                                         cGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGCTGGCTATGCCAGACGACGTTTTCGTCAt                                                                                                                                                                                                             >  4:512650/1‑167 (MQ=255)
                                                                                                                                                                            aCTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAga                                                                                                                                                                                                                            <  8:768904/149‑1 (MQ=255)
                                                                                                                                                                             cTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAg                                                                                                                                                                                                         <  2:512650/167‑1 (MQ=255)
                                                                                                                                                                              tGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAgacg                                                                                                                                                                                                                          >  5:873749/1‑149 (MQ=255)
                                                                                                                                                                                  gCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCCCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGacgacgc                                                                                                                                                                                                                      >  6:23600/1‑148 (MQ=255)
                                                                                                                                                                                   acGATTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTcgcg                                                                                                                                                                                                                                                            <  3:140824/109‑1 (MQ=255)
                                                                                                                                                                                    cGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAGCCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCaga                                                                                                                                          <  4:433671/223‑1 (MQ=255)
                                                                                                                                                                                    cGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGttt                                                                                                                                                                                                                    <  6:75568/149‑1 (MQ=255)
                                                                                                                                                                                    cGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGttt                                                                                                                                                                                                                    >  8:483626/1‑149 (MQ=255)
                                                                                                                                                                                      aGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTg                                                                                     >  2:437860/1‑274 (MQ=255)
                                                                                                                                                                                        tttATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGt                                                                                                                                                                                                                >  5:44149/1‑149 (MQ=255)
                                                                                                                                                                                           aTGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCAt                                                                                                                                                                                                             <  6:965052/149‑1 (MQ=255)
                                                                                                                                                                                            tGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATg                                                                                                                                                                                                            <  8:965747/149‑1 (MQ=255)
                                                                                                                                                                                                 cGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGa                                                                                                                                                            >  1:130128/1‑192 (MQ=255)
                                                                                                                                                                                                  gCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAgcgc                                                                                                                                                                                                      >  7:210442/1‑149 (MQ=255)
                                                                                                                                                                                                  gCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTg                                                                                                                                                       >  4:658555/1‑196 (MQ=255)
                                                                                                                                                                                                     tcctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCat                                                                                                                                                                                                   >  5:270715/1‑149 (MQ=255)
                                                                                                                                                                                                     tcctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGt                                                                                                                                                        <  3:130128/192‑1 (MQ=255)
                                                                                                                                                                                                      cctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGcga                                                                                                                                                   <  2:658555/196‑1 (MQ=255)
                                                                                                                                                                                                       ctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGAt                                                                                                                                                                                                                                        >  5:1012103/1‑110 (MQ=255)
                                                                                                                                                                                                       ctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCatta                                                                                                                                                                                                 >  5:823885/1‑149 (MQ=255)
                                                                                                                                                                                                           aGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAttattta                                                                                                                                                                                             <  8:407278/149‑1 (MQ=255)
                                                                                                                                                                                                              gcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATcc                                                                                                                                                                                          <  6:873749/149‑1 (MQ=255)
                                                                                                                                                                                                                cccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATcccc                                                                                                                                                                                        <  7:483626/149‑1 (MQ=255)
                                                                                                                                                                                                                  cGGCGACCGTTGATTTACCGGACCGGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCa                                                                                                                                                                                      <  6:803905/149‑1 (MQ=255)
                                                                                                                                                                                                                  cGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAgc                                                                                                                           >  1:135710/1‑207 (MQ=255)
                                                                                                                                                                                                                          gTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCatta                                                                                                                                                                                                 <  8:751217/130‑1 (MQ=255)
                                                                                                                                                                                                                          gTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCatta                                                                                                                                                                                                 >  7:751217/1‑130 (MQ=255)
                                                                                                                                                                                                                          gTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGAGGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTATTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcg                                                                                                                                                                              >  8:592264/1‑149 (MQ=255)
                                                                                                                                                                                                                           ttGATTTTCCGGACCCGGGGAGCCCGGTAAACCACAGCACTACACCCCGATGACCGGGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTca                                               <  4:460479/275‑1 (MQ=255)
                                                                                                                                                                                                                           ttGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcgc                                                                                                                                                                             >  7:645955/1‑149 (MQ=255)
                                                                                                                                                                                                                            tGATTTACCGGACCCGGGGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTca                                               <  2:11757/274‑1 (MQ=255)
                                                                                                                                                                                                                                  aCCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGcc                                                                                                                                                                      >  8:242095/1‑149 (MQ=255)
                                                                                                                                                                                                                                   ccGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTACCCCCCAGCAAATCGCGCGCGccc                                                                                                                                                                     <  6:270715/149‑1 (MQ=255)
                                                                                                                                                                                                                                     ggACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGAt                                                                                                      >  2:310230/1‑210 (MQ=255)
                                                                                                                                                                                                                                         ccGGAGAGGCCGGTAAACCACAGCACTACACCACGATCACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGc                                                                                                                                                               >  6:101275/1‑149 (MQ=255)
                                                                                                                                                                                                                                          cGGAGAGTCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCg                                                                                                                                                              >  6:165089/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 gCCGGTAAACCACAGCACTACACCACGATGACCgtggt                                                                                                                                                                                                                                                                      >  8:558107/1‑38 (MQ=255)
                                                                                                                                                                                                                                                  ccGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTgg                                                                                                                                                      <  8:938401/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   cGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGc                                                                                                                                                     >  5:557528/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    ggTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGcg                                                                                                                                                    >  7:389657/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         aCCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAgc                                                                                                                                                                                                        <  8:983266/92‑1 (MQ=255)
                                                                                                                                                                                                                                                          ccACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAAc                                                                                                                                              <  5:119133/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           cacaGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAAcc                                                                                                                                             >  6:276589/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              aGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGATCGATAAAAATCAGCCCAcc             <  2:63020/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                cACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGc                                                                                                                                        <  7:558263/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                cACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTcc                                                                                                                              >  4:729065/1‑159 (MQ=255)
                                                                                                                                                                                                                                                                  cTACACCACGGTCATCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAgc                                                                                                                                                                                                        >  5:457054/1‑83 (MQ=38)
                                                                                                                                                                                                                                                                  cTACACCACGGTCATCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCAt                                                                                                                                      >  5:995549/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  cTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCAt                                                                                                                                      >  7:872348/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   tACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAga                                                                                                                                                                                                                            <  8:463865/62‑1 (MQ=255)
                                                                                                                                                                                                                                                                      accacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCa                                                                                                                                       <  5:652861/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                      accacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCa                                                                                                                                       >  6:652861/1‑144 (MQ=255)
                                                                                                                                                                                                                                                                      accacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCaa                                                                                                                                  >  8:112580/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        cacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAtt                                                                                                                                <  8:633077/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                        cacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGt                                                 >  2:160014/1‑228 (MQ=255)
                                                                                                                                                                                                                                                                         cccaTGACCGTGGTGTAGCTTGCGTTGTTGCAAAGTAACCGGATGGCAATGCCAGACGCCGTTTTCGTCATGCAGCGCCATTATTTATCCCCCCGCAAATCGCGCGCGCCCCACTGCGGAAAGTTGCGACTAACCAGAGCATTCAATTc                                                                                                                               <  8:210442/146‑1 (MQ=255)
                                                                                                                                                                                                                                                                         acGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCaga                                                                                                                                          <  2:742511/138‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGGGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCt                                                                                  >  2:33155/1‑193 (MQ=255)
                                                                                                                                                                                                                                                                          cGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTcc                                                                                                                              <  6:450628/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCa                                                                                                                                                                                      <  6:86615/92‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCa                                                                                                                             <  6:1020548/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCa                                                                                                                             <  5:622573/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aTGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTcaca                                             <  4:160014/228‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aTGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGATCGATAAAAATCAGCCCACCCGTCa        <  2:546035/265‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGGGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACt                                                                              <  4:33155/193‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGtt                                                                                                                          >  6:931782/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                          
GGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGATCGATAAAAATCAGCCCACCCGTCACCGGAT  >  NC_000913/2873647‑2874185

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