Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,937,983 T→C 6.4% I112T (ATC→ACC)  pykA → pyruvate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,937,9830TC 6.4% 372.1 / 10.3 109I112T (ATC→ACC) pykApyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base T (53/49);  new base C (3/4);  total (56/53)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 110 total aligned reads displayed.

TGGCTATTCTGGGTGACCTCCAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGT  >  NC_000913/1937842‑1938131
                                                                                                                                             |                                                                                                                                                    
tgGCTATTCTGGGTGACCTCCAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGATACCGACttt                                                                                                                                               <  3:1442250/149‑3 (MQ=255)
 ggCTATTCTGGGTGACCTCCAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGATACCGACttt                                                                                                                                               >  2:734115/1‑146 (MQ=255)
 ggCTATTCTGGGTGACCTCCAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGGTTTCCGCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGATACCGACttt                                                                                                                                               <  1:734115/148‑3 (MQ=255)
                  tCCAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGgcctgc                                                                                                                                     <  1:1005319/141‑1 (MQ=255)
                       gggCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGAcgtcgt                                                                                                                        <  1:63618/149‑1 (MQ=255)
                         gCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGc                                                                                                                      >  4:804068/1‑149 (MQ=255)
                          cccAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGa                                                                                                                              >  3:362085/1‑140 (MQ=255)
                          cccAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGa                                                                                                                              <  4:362085/140‑1 (MQ=255)
                           ccAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGgc                                                                                                                                         <  4:859713/128‑1 (MQ=255)
                           ccAAAATCCGTATATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGgc                                                                                                                                         >  3:859713/1‑128 (MQ=255)
                                 tCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGAc                                                                                                              >  3:1294709/1‑149 (MQ=255)
                                     tgtATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGAcg                                                                                                                            <  3:825478/131‑1 (MQ=255)
                                     tgtATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGAcg                                                                                                                            >  4:825478/1‑131 (MQ=255)
                                     tgtATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATcc                                                                                                          >  4:1325185/1‑149 (MQ=255)
                                     tgtATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATcc                                                                                                          >  4:682966/1‑149 (MQ=255)
                                      gtATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCct                                                                                                         <  3:682966/149‑1 (MQ=255)
                                      gtATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCct                                                                                                         <  3:1325185/149‑1 (MQ=255)
                                            aCCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCGACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATc                                                                                                           <  4:833963/141‑1 (MQ=255)
                                            aCCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGgcctgcc                                                                                                                                    >  4:1020070/1‑116 (MQ=255)
                                            aCCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGgcctgcc                                                                                                                                    >  4:572571/1‑116 (MQ=255)
                                            aCCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGgcctgcc                                                                                                                                    <  3:572571/116‑1 (MQ=255)
                                            aCCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGgcctgcc                                                                                                                                    <  3:1020070/116‑1 (MQ=255)
                                            aCCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATc                                                                                                           >  3:833963/1‑141 (MQ=255)
                                            aCCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCctgctgct                                                                                                   >  2:1255572/1‑149 (MQ=255)
                                               tttAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGa                                                                                                >  2:166645/1‑149 (MQ=255)
                                               tttAAAGAAAGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGa                                                                                                <  1:1069910/149‑1 (MQ=255)
                                                   aagaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGAcg                                                                                                                            <  1:315867/117‑1 (MQ=255)
                                                   aagaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGAcg                                                                                                                            <  1:68326/117‑1 (MQ=255)
                                                   aagaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGAcg                                                                                                                            >  2:315867/1‑117 (MQ=255)
                                                   aagaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGAcg                                                                                                                            >  2:68326/1‑117 (MQ=255)
                                                    agaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATc                                                                                                                                                     >  2:1432857/1‑91 (MQ=255)
                                                    agaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATc                                                                                                                                                     <  1:1432857/91‑1 (MQ=255)
                                                           aaaGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGt                                                                                                                 <  3:752577/120‑1 (MQ=255)
                                                           aaaGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGt                                                                                                                 >  4:752577/1‑120 (MQ=255)
                                                           aaaGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTg                                                                                                                >  4:1382744/1‑121 (MQ=255)
                                                           aaaGTTTTCCTCAATATTGGGGATAAATTCCTGCTCCACGCTAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTg                                                                                                                <  3:1382744/121‑1 (MQ=255)
                                                              gTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCa                                                                                 <  3:804068/149‑1 (MQ=255)
                                                              gTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCa                                                                                 <  4:499097/149‑1 (MQ=255)
                                                                    cTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGa                                                                                                                                                   >  4:1106329/1‑77 (MQ=255)
                                                                    cTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGa                                                                                                                                                   <  3:1106329/77‑1 (MQ=255)
                                                                     tCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTaaaa                                                                          <  4:1294709/149‑1 (MQ=255)
                                                                     tCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTaaaa                                                                          >  3:1057612/1‑149 (MQ=255)
                                                                      cAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTaaaa                                                                          <  3:1067160/148‑1 (MQ=255)
                                                                      cAATATTGGGGATAAATTCCTGCTCGACGCCAACCTCGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTaaaa                                                                          >  4:1067160/1‑148 (MQ=255)
                                                                         tatTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGa                                                                                                                                                   <  4:1142529/72‑1 (MQ=255)
                                                                         tatTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGa                                                                                                                                                   >  3:1142529/1‑72 (MQ=255)
                                                                             ggggATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGt                                                                                         >  4:1459598/1‑126 (MQ=255)
                                                                             ggggATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGt                                                                                         <  3:1459598/126‑1 (MQ=255)
                                                                              gggATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAATGGCctgcctgc                                                                                                                                 <  3:611658/85‑1 (MQ=255)
                                                                              gggATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAATGGCctgcctgc                                                                                                                                 >  4:611658/1‑85 (MQ=255)
                                                                              gggATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGaa                                                                 >  2:928353/1‑149 (MQ=255)
                                                                               ggATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCctgcctgc                                                                                                                                 <  1:272962/84‑1 (MQ=255)
                                                                                  tAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTc                                                                                        >  3:566420/1‑122 (MQ=255)
                                                                                    aaTTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAg                                                           <  4:1057612/149‑1 (MQ=255)
                                                                                     aTTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAgg                                                          >  4:117559/1‑149 (MQ=255)
                                                                                         cTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTc                                                             >  3:113183/1‑142 (MQ=255)
                                                                                         cTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTc                                                             <  4:113183/142‑1 (MQ=255)
                                                                                         cTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCAt                                                      >  3:968591/1‑149 (MQ=255)
                                                                                          tGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTgg                                                                                                                  <  3:128932/88‑1 (MQ=255)
                                                                                          tGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTgg                                                                                                                  >  4:128932/1‑88 (MQ=255)
                                                                                             tCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGaaa                                                  >  4:637199/1‑149 (MQ=255)
                                                                                             tCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGaaa                                                  <  3:117559/149‑1 (MQ=255)
                                                                                                  gCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGAcgt                                                                                                                           >  2:361758/1‑71 (MQ=255)
                                                                                                  gCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGtt                                             <  2:1232929/149‑1 (MQ=255)
                                                                                                            gTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTc                                   >  2:961741/1‑149 (MQ=255)
                                                                                                                  gggAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTc                             >  2:233157/3‑149 (MQ=255)
                                                                                                                   tGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGaaa                                                  <  2:427322/127‑1 (MQ=255)
                                                                                                                   tGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGaaa                                                  >  1:427322/1‑127 (MQ=255)
                                                                                                                        gCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTc                                                                                        <  4:1135293/84‑1 (MQ=255)
                                                                                                                        gCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTc                                                                                        >  3:1135293/1‑84 (MQ=255)
                                                                                                                        gCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTaaaa                                                                          <  1:1293925/98‑1 (MQ=255)
                                                                                                                        gCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTaaaa                                                                          >  2:1293925/1‑98 (MQ=255)
                                                                                                                          gACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGaaa                                                  >  4:122192/1‑120 (MQ=255)
                                                                                                                          gACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGaa                                      >  4:1137270/1‑132 (MQ=255)
                                                                                                                          gACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGaa                                      <  3:1137270/132‑1 (MQ=255)
                                                                                                                          gACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTcc                     >  2:405189/1‑149 (MQ=255)
                                                                                                                          gACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCAGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGaaa                                                  <  3:122192/120‑1 (MQ=255)
                                                                                                                           aCAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGcc                                                                                                                     <  1:954033/52‑1 (MQ=255)
                                                                                                                               aGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTc                                                                                        >  1:1366154/1‑77 (MQ=255)
                                                                                                                               aGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTc                                                                                        <  2:1366154/77‑1 (MQ=255)
                                                                                                                               aGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAggg                                                         >  1:550820/1‑108 (MQ=255)
                                                                                                                               aGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAggg                                                         <  2:550820/108‑1 (MQ=255)
                                                                                                                               aGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCctct                >  3:739461/1‑149 (MQ=255)
                                                                                                                                tttgtaGTCGATACCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCggtg                         >  1:1197876/6‑139 (MQ=255)
                                                                                                                                tttgtaGTCGATACCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCggtg                         >  1:804475/6‑139 (MQ=255)
                                                                                                                                tttgtaGTCGATACCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCggtg                         <  2:1197876/134‑1 (MQ=255)
                                                                                                                                tttgtaGTCGATACCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCggtg                         <  2:804475/134‑1 (MQ=255)
                                                                                                                                gAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTgg                                                                   >  4:1379008/1‑97 (MQ=255)
                                                                                                                                gAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTgg                                                                   <  3:1379008/97‑1 (MQ=255)
                                                                                                                                gAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCtctc               <  1:928353/149‑1 (MQ=255)
                                                                                                                                gAAAAAGGCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCtctc               >  2:783476/1‑149 (MQ=255)
                                                                                                                                 aaaaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAg                                                           <  4:497306/104‑1 (MQ=255)
                                                                                                                                 aaaaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAg                                                           >  3:497306/1‑104 (MQ=255)
                                                                                                                                  aaaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGaa                                                                 <  2:1282187/97‑1 (MQ=255)
                                                                                                                                  aaaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGaa                                                                 >  1:1282187/1‑97 (MQ=255)
                                                                                                                                   aaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCcaa            >  1:219959/1‑149 (MQ=255)
                                                                                                                                       tCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAAt        <  4:739461/149‑1 (MQ=255)
                                                                                                                                          gTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATaaa     <  1:783476/149‑1 (MQ=255)
                                                                                                                                          gTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATaaa     <  2:219959/149‑1 (MQ=255)
                                                                                                                                             tCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGt  <  1:961741/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
TGGCTATTCTGGGTGACCTCCAGGGGCCCAAAATCCGTGTATCCACCTTTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGT  >  NC_000913/1937842‑1938131

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