Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,873,919 C→T 6.4% 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 6.4% 650.7 / 24.7 204G25D (GGT→GAT) cysCadenosine 5'‑phosphosulfate kinase
Reads supporting (aligned to +/- strand):  ref base C (94/95);  new base T (5/8);  total (100/104)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-01

Only 100 of 207 total aligned reads displayed.

GCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGATCGATAAAAATCAGCCCACCCGTCACCGGATTTTGT  >  NC_000913/2873648‑2874190
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                               
gtgTTCGGCGGGGGTGGGCGGTATAATGGCGGTATGCACCACCCGTCCGTCTTCAACCATCAAATTCGCCACGTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATCATCCGGGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCg                                                                                                                                                                                                                                                                               <  4:246288/272‑1 (MQ=255)
           ggTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTcg                                                                                                                                                                                                                                                                   <  2:280673/275‑1 (MQ=255)
                                   gCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGca                                                                                                                                                                                                                                                         <  2:334332/261‑1 (MQ=255)
                                              ccGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAg                                                                                                                                                                                                                                                                       >  4:12266/1‑236 (MQ=255)
                                                                                                          gCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTTTAGCCCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCCTGCAGCGCCATTATTTAACCCCCAGCAAACCGCCCCCCCCCCAGtcc                                                                                                                                                                    >  4:32083/1‑273 (MQ=255)
                                                                                                                         agaTCGCTGCAGACACCGTGGCGAACATTCACGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCGCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGTAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGcgacg                                                                                                                                                      <  2:318439/274‑1 (MQ=255)
                                                                                                                                tccATAACCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGGGTTTATGTAACGCCTCCTCCAGCGCCCCCGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAACACTACACCACGATAACCctgg                                                                                                                                                                                                                                                                            >  6:636525/3‑149 (MQ=255)
                                                                                                                                tGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtgg                                                                                                                                                                                                                                                                            <  5:806134/149‑1 (MQ=255)
                                                                                                                                 gCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtggt                                                                                                                                                                                                                                                                           <  8:1349904/149‑1 (MQ=255)
                                                                                                                                   aTAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGgtgt                                                                                                                                                                                                                                                                         <  5:431281/149‑1 (MQ=255)
                                                                                                                                    tAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTa                                                                                                                                                                                                                                                                        <  7:1220064/149‑1 (MQ=255)
                                                                                                                                    tAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTa                                                                                                                                                                                                                                                                        <  6:1218900/149‑1 (MQ=255)
                                                                                                                                        ccGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTc                                                                                                                                                                                                                                                                    >  7:1047960/1‑149 (MQ=255)
                                                                                                                                            ggCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTCTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACCAACCAGAGCATTCAATTCCa                                                                                                                                  >  2:185880/1‑275 (MQ=255)
                                                                                                                                              cGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttg                                                                                                                                                                                                                                                              <  5:1007114/149‑1 (MQ=255)
                                                                                                                                               gAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgt                                                                                                                                                                                                                                                             <  6:814831/149‑1 (MQ=255)
                                                                                                                                                aaCATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgtt                                                                                                                                                                                                                                                            >  6:138983/1‑149 (MQ=255)
                                                                                                                                                aaCATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgtt                                                                                                                                                                                                                                                            <  6:100981/149‑1 (MQ=255)
                                                                                                                                                aaCATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgtt                                                                                                                                                                                                                                                            <  7:457191/149‑1 (MQ=255)
                                                                                                                                                     tGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCTCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGcaca                                                                                                                                                                                                                                                       <  5:361505/149‑1 (MQ=255)
                                                                                                                                                     tGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTcg                                                                                                                                                                                                                                                                   <  4:373415/137‑1 (MQ=255)
                                                                                                                                                     tGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGcaca                                                                                                                                                                                                                                                       <  6:1104480/149‑1 (MQ=255)
                                                                                                                                                       tCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtgg                                                                                                                                                                                                                                                                            >  7:782327/1‑126 (MQ=255)
                                                                                                                                                       tCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGt                                                                                                                                                                                                                                                     <  5:790542/149‑1 (MQ=255)
                                                                                                                                                       tCGCCATCCAGCAGATACGTACTGACGCCCAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCa                                                                                                                       <  2:38671/275‑1 (MQ=255)
                                                                                                                                                          ccATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGAc                                                                                                                                                                                                                                                  >  6:51760/1‑149 (MQ=255)
                                                                                                                                                               cagcagATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGAt                                                                                                                                                                                                                                             <  7:1087495/149‑1 (MQ=255)
                                                                                                                                                               cagcagATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGAt                                                                                                                                                                                                                                             <  8:699075/149‑1 (MQ=255)
                                                                                                                                                               cagcagATACGTACTGACGCCGAATTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtggt                                                                                                                                                                                                                                                                           <  7:950945/119‑1 (MQ=255)
                                                                                                                                                                     aTACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCg                                                                                                                                                                                                                                                <  6:1473531/140‑1 (MQ=255)
                                                                                                                                                                     aTACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTNt                                                                                                                                                                                                                                       <  7:340928/149‑1 (MQ=255)
                                                                                                                                                                     aTACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAgacg                                                                                                                                                                                                                               <  1:53701/157‑1 (MQ=255)
                                                                                                                                                                         gTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCa                                                                                                                                                                                                                                   >  7:1415237/1‑149 (MQ=255)
                                                                                                                                                                            cTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAgac                                                                                                                                                                                                                                >  5:231810/1‑149 (MQ=255)
                                                                                                                                                                            cTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAgac                                                                                                                                                                                                                                <  5:685905/149‑1 (MQ=255)
                                                                                                                                                                             tGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtgg                                                                                                                                                                                                                                                                            <  7:975909/104‑1 (MQ=255)
                                                                                                                                                                               aCGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTg                                                                                                                                                                     >  2:87024/1‑205 (MQ=255)
                                                                                                                                                                               aCGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGt                                                                                                >  2:206330/1‑274 (MQ=255)
                                                                                                                                                                                  ccGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGca                                                                                                                                                                                                                                                         <  6:872685/118‑1 (MQ=255)
                                                                                                                                                                                  ccGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGtt                                                                                                                                                                                                                          >  7:161574/1‑149 (MQ=255)
                                                                                                                                                                                       tttATGTAACGCCTCCTCCATCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACTGTGACCGGATGGCTATGCCAGACGACGTTTTCGt                                                                                                                                                                                                                     >  8:797731/1‑149 (MQ=255)
                                                                                                                                                                                       tttATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTc                                                                                                                                                                                                                                                                    <  8:387259/102‑1 (MQ=255)
                                                                                                                                                                                          aTGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTa                                                                                                                                                                                                                                                                        >  5:173929/1‑95 (MQ=255)
                                                                                                                                                                                           tGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATg                                                                                                                                                                                                                 <  5:91434/149‑1 (MQ=255)
                                                                                                                                                                                           tGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACt                                                                                   <  2:170478/275‑1 (MQ=255)
                                                                                                                                                                                           actctccagtactccAGCTCCCCGGCCAACATTGATTTACCGGACCCGGAGGGGTCGGTAAACCACAGCACTACACCCCTATGACCGTGGTGTAGCTCGCGTTGCTGCACAGTGACCGGATTGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCGAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACt                                                                                   <  1:31813/264‑1 (MQ=255)
                                                                                                                                                                                             tAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCa                                                                                                                                                                                                               >  5:1296970/1‑149 (MQ=255)
                                                                                                                                                                                              aaCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTAgt                                                                                                                                                                                                                                                                      <  8:13053/93‑2 (MQ=39)
                                                                                                                                                                                              aaCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTAgt                                                                                                                                                                                                                                                                      <  8:1160992/93‑2 (MQ=38)
                                                                                                                                                                                              aaCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTAgt                                                                                                                                                                                                                                                                      >  7:1160992/1‑92 (MQ=39)
                                                                                                                                                                                               aCGCCTCCTCCGGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAgc                                                                                                                                                                                                             >  7:477177/1‑149 (MQ=255)
                                                                                                                                                                                               aCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCa                                                                                                                                                                                                                                   <  5:255353/127‑1 (MQ=255)
                                                                                                                                                                                               aCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAgc                                                                                                                                                                                                             >  7:1319217/1‑149 (MQ=255)
                                                                                                                                                                                               aCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTCGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAgc                                                                                                                                                                                                             >  8:1101686/1‑149 (MQ=255)
                                                                                                                                                                                                 gCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAgcgc                                                                                                                                                                                                           >  8:394751/1‑149 (MQ=255)
                                                                                                                                                                                                   ctcctcCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCa                                                                                                                                                                                                         >  6:159539/1‑149 (MQ=255)
                                                                                                                                                                                                    tcctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCat                                                                                                                                                                                                        >  8:652476/1‑149 (MQ=255)
                                                                                                                                                                                                    tcctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGCTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCat                                                                                                                                                                                                        >  6:910041/1‑149 (MQ=255)
                                                                                                                                                                                                      ctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCatta                                                                                                                                                                                                      >  5:1233222/1‑149 (MQ=255)
                                                                                                                                                                                                      ctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGCCATGCAGCGCCatta                                                                                                                                                                                                      >  6:530468/1‑149 (MQ=255)
                                                                                                                                                                                                       tccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCattat                                                                                                                                                                                                     >  6:687794/1‑149 (MQ=255)
                                                                                                                                                                                                       tccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCattat                                                                                                                                                                                                     <  8:199689/149‑1 (MQ=255)
                                                                                                                                                                                                           gcgcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAttatttat                                                                                                                                                                                                 <  6:231810/149‑1 (MQ=255)
                                                                                                                                                                                                            cgcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATc                                                                                                                                                                                                <  5:741887/149‑1 (MQ=255)
                                                                                                                                                                                                               cccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATcccc                                                                                                                                                                                             >  5:503824/1‑149 (MQ=255)
                                                                                                                                                                                                                ccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATccccc                                                                                                                                                                                            <  5:1482476/149‑1 (MQ=255)
                                                                                                                                                                                                                ccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATccccc                                                                                                                                                                                            <  5:881621/149‑1 (MQ=255)
                                                                                                                                                                                                                  ggCGACCGTTGATTTACCGGACCCGGGGAGGCCGGGAAACAACAGGACAAAACCACGTTGACCGTGGGGTAGCTGGCGTTGTGGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGGGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCAcc                                                            <  4:188293/275‑1 (MQ=255)
                                                                                                                                                                                                                  ggCGACCGTTGAGTTACCGGACCCGGAGAGGCCGGTAAACAACAGCACTCCACCCAGATGACCGTTGTGTTGCTCGCGTTGTTGCACAGTGACCGGCTGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCAcc                                                            <  4:206330/275‑1 (MQ=255)
                                                                                                                                                                                                                        cGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATc                                                                                                                                                                                    >  5:609348/1‑149 (MQ=255)
                                                                                                                                                                                                                         gTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcg                                                                                                                                                                                   >  7:387950/1‑149 (MQ=255)
                                                                                                                                                                                                                         gTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcg                                                                                                                                                                                   >  7:879309/1‑149 (MQ=255)
                                                                                                                                                                                                                          ttGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATACCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcgc                                                                                                                                                                                  >  5:439207/1‑149 (MQ=255)
                                                                                                                                                                                                                             aTTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcgcgcg                                                                                                                                                                               >  5:1035824/1‑149 (MQ=255)
                                                                                                                                                                                                                                 aCCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGcc                                                                                                                                                                           >  5:1362598/1‑149 (MQ=255)
                                                                                                                                                                                                                                   cGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATcccc                                                                                                                                                                                             >  7:764361/1‑129 (MQ=255)
                                                                                                                                                                                                                                   cGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATcccc                                                                                                                                                                                             <  8:764361/129‑1 (MQ=255)
                                                                                                                                                                                                                                     gACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAg                                                                                                                                                                       >  5:868587/1‑149 (MQ=255)
                                                                                                                                                                                                                                        ccGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGc                                                                                                                                                                    <  5:1401763/149‑1 (MQ=255)
                                                                                                                                                                                                                                         cGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCGCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGcgat                                     <  2:277155/275‑1 (MQ=255)
                                                                                                                                                                                                                                             gagGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGaaa                                                                                                                                                               <  5:599431/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   ggTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGcg                                                                                                                                                         >  8:939889/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    gTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGcga                                                                                                                                                        >  5:827528/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    gTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGcga                                                                                                                                                        <  7:394751/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     tAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGcgac                                                                                                                                                       >  8:287226/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        aCCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGaa                                                                                                                                                    >  6:439606/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         ccACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAAc                                                                                                                                                   <  8:1115231/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         ccACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAAc                                                                                                                                                   <  6:407330/149‑1 (MQ=255)
                                                                                                                                                                                                                                                          cacaGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCatta                                                                                                                                                                                                      >  8:156138/1‑97 (MQ=255)
                                                                                                                                                                                                                                                           acaGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCa                                                                                                                                                 >  7:796417/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               cACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGc                                                                                                                                             >  5:1332725/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                aCTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCa                                                                                                                                            <  8:1184702/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                aCTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCa                                                                                                                                            >  8:716254/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 cTACACCACGGTCATCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCAt                                                                                                                                           >  6:289717/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  tACACCACGGTCATCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCAtt                                                                                                                                          >  5:1452248/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    caccacGGTCATCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCa                                                                                                                                        >  7:12792/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     accacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCaa                                                                                                                                       <  6:756042/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ccacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTg                                                                                                                                                            >  7:1333185/1‑127 (MQ=255)
                                                                                                                                                                                                                                                                      ccacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAt                                                                                                                                      >  8:465047/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      ccacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAt                                                                                                                                      <  6:366894/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                               
GCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGATCGATAAAAATCAGCCCACCCGTCACCGGATTTTGT  >  NC_000913/2873648‑2874190

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