Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,671,759 C→A 72.7% R179R (CGG→AGG)  hcaB → 2,3‑dihydroxy‑2,3‑dihydrophenylpropionate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,671,7590CA72.7% 219.5 / 118.0 146R179R (CGG→AGG) hcaB2,3‑dihydroxy‑2,3‑dihydrophenylpropionate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (22/18);  new base A (53/53);  total (75/71)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.19e-01

Only 100 of 146 total aligned reads displayed.

CTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGC  >  NC_000913/2671662‑2671855
                                                                                                 |                                                                                                
       ccTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATg                                                                                         >  2:1092500/1‑100 (MQ=255)
       ccTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATg                                                                                         <  1:4010997/100‑1 (MQ=255)
       ccTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATg                                                                                         <  2:1136323/100‑1 (MQ=255)
       ccTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATg                                                                                         <  1:2353900/100‑1 (MQ=255)
       ccTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATg                                                                                         <  1:1880531/100‑1 (MQ=255)
        cTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATgg                                                                                        >  1:1065898/1‑100 (MQ=255)
        cTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATgg                                                                                        >  1:4008162/1‑100 (MQ=255)
        cTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATgg                                                                                        <  1:3176804/100‑1 (MQ=255)
        cTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGCGAATgg                                                                                        >  1:3577696/1‑100 (MQ=255)
          ggcggcGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCg                                                                                      >  1:3866165/1‑100 (MQ=255)
          ggcggcGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCg                                                                                      <  1:3246576/100‑1 (MQ=255)
          ggcggcGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCg                                                                                      <  2:4123749/100‑1 (MQ=255)
            cggcggTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTc                                                                                    <  1:2278244/100‑1 (MQ=255)
            cggcggTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTc                                                                                    >  1:3079222/1‑100 (MQ=255)
             ggcggTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCg                                                                                   >  2:3723826/1‑100 (MQ=255)
              gcggTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCgg                                                                                  <  2:1303949/100‑1 (MQ=255)
              gcggTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCgg                                                                                  >  1:1553829/1‑100 (MQ=255)
              gcggTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCgg                                                                                  <  2:769618/100‑1 (MQ=255)
                gatGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGcc                                                                                >  2:189861/3‑100 (MQ=255)
                  tGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCg                                                                              <  2:188981/100‑1 (MQ=255)
                    gCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAGTGGCGTCGGCCCgtg                                                                            >  2:768166/1‑100 (MQ=255)
                     cccGCTGTACACCGCCAGTAAACATGCCGCAACCCGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAACTGAGGGTGAATCGc                                                                                       >  2:2513474/1‑88 (MQ=255)
                      ccGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGtgtg                                                                          <  1:414535/100‑1 (MQ=255)
                      ccGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGtgtg                                                                          >  2:680022/1‑100 (MQ=255)
                            tACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGGTTAGGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATgg                                                                    <  2:43891/100‑1 (MQ=255)
                            tACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATgg                                                                    <  1:1055653/100‑1 (MQ=255)
                              cacCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGcc                                                                  >  1:4175683/1‑100 (MQ=255)
                              cacCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGcc                                                                  <  1:3539047/100‑1 (MQ=255)
                               acCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGcc                                                                  <  1:3477113/99‑1 (MQ=255)
                                ccGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAg                                                                >  1:3348711/1‑100 (MQ=255)
                                  gCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCATCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCg                                                              <  1:780634/100‑1 (MQ=255)
                                  gCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAg                                                                >  2:1690688/1‑98 (MQ=255)
                                  gCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCg                                                              >  2:2024047/1‑100 (MQ=255)
                                  gCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCg                                                              >  2:135639/1‑100 (MQ=255)
                                   ccAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGa                                                             <  1:1950101/100‑1 (MQ=255)
                                   ccAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGACTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGa                                                             >  2:2563015/1‑100 (MQ=255)
                                    cAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGAc                                                            <  2:3686268/100‑1 (MQ=255)
                                    cAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGAc                                                            >  1:3657484/1‑100 (MQ=255)
                                     aGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGAcc                                                           <  2:368739/100‑1 (MQ=255)
                                     aGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGAcc                                                           <  1:2226498/100‑1 (MQ=255)
                                        aaaCATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGAACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTgc                                                        >  2:3156664/1‑100 (MQ=255)
                                         aaCATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTgcg                                                       <  1:1263295/100‑1 (MQ=255)
                                         aaCATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGa                                                             >  2:2269817/1‑94 (MQ=255)
                                           cATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGGGGTATGGCCAGCGACCTGcgcg                                                     >  1:1409796/1‑100 (MQ=255)
                                             tGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGc                                                   <  2:1399931/100‑1 (MQ=255)
                                              gCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGcc                                                  >  2:2204461/1‑100 (MQ=255)
                                              gCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGcc                                                  <  1:797260/100‑1 (MQ=255)
                                                cGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCca                                                >  2:711026/1‑100 (MQ=255)
                                                  cAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCcaca                                              >  2:4110334/1‑100 (MQ=255)
                                                  cAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCcaca                                              <  2:3066202/100‑1 (MQ=255)
                                                  cAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCAACCTGCGCGGCCcaca                                              <  2:2328333/100‑1 (MQ=255)
                                                    aCCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAgg                                            <  2:560859/100‑1 (MQ=255)
                                                     ccGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGgc                                           <  1:3109444/100‑1 (MQ=255)
                                                       ggACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGgcgc                                         >  1:204850/1‑100 (MQ=255)
                                                         aCTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTc                                       >  1:927922/1‑100 (MQ=255)
                                                          cTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCg                                      >  1:2689964/1‑100 (MQ=255)
                                                           ttATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCgg                                     <  1:1830524/100‑1 (MQ=255)
                                                            tATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCggg                                    <  2:3929739/100‑1 (MQ=255)
                                                              ttCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCa                                  >  1:2523410/1‑100 (MQ=255)
                                                                cGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCaaa                                >  1:2366524/1‑100 (MQ=255)
                                                                cGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCaaa                                >  1:904355/1‑100 (MQ=255)
                                                                cGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCaaa                                >  1:1121465/1‑100 (MQ=255)
                                                                cGCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCaaa                                >  2:607898/1‑100 (MQ=255)
                                                                 nccAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAg                               >  2:456361/2‑100 (MQ=255)
                                                                 gCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAg                               <  1:1422559/100‑1 (MQ=255)
                                                                 gCCAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAg                               <  1:1389813/100‑1 (MQ=255)
                                                                   cAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTg                             <  2:1981193/100‑1 (MQ=255)
                                                                   cAACTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTg                             >  2:3912270/1‑100 (MQ=255)
                                                                     aCTGGCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCGAAGTGaa                           <  2:1131681/100‑1 (MQ=255)
                                                                      cTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGaaa                          >  2:3347709/1‑100 (MQ=255)
                                                                      cTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGCGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGaaa                          >  2:247599/1‑100 (MQ=255)
                                                                        ggCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAAcc                        >  1:2339002/1‑100 (MQ=255)
                                                                        ggCTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAAcc                        <  1:820807/100‑1 (MQ=255)
                                                                          cTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTc                      <  1:2308548/100‑1 (MQ=255)
                                                                          cTTATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCt                       >  2:1749393/1‑99 (MQ=255)
                                                                          attATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTc                      <  2:347787/99‑1 (MQ=255)
                                                                            tATGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGa                    >  1:3259221/1‑100 (MQ=255)
                                                                             aTGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGAt                   <  2:2331499/100‑1 (MQ=255)
                                                                              tGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATa                  <  2:2849249/100‑1 (MQ=255)
                                                                              tGAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATa                  <  2:3528030/100‑1 (MQ=255)
                                                                               gAACTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATaa                 >  2:3931083/1‑100 (MQ=255)
                                                                                 aCTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATg               <  1:3214760/100‑1 (MQ=255)
                                                                                  cTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGc              <  1:1503295/100‑1 (MQ=255)
                                                                                  cTGGCACCGAAAGTGAGGGTGCATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGc              <  2:1503121/100‑1 (MQ=255)
                                                                                  cTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGc              >  1:2026938/1‑100 (MQ=255)
                                                                                  cTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGc              >  1:2552208/1‑100 (MQ=255)
                                                                                  cTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGc              <  1:3109511/100‑1 (MQ=255)
                                                                                  cTGGCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGc              >  2:1691214/1‑100 (MQ=255)
                                                                                    ggCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAg            >  2:1238467/1‑100 (MQ=255)
                                                                                    ggCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAg            >  2:388087/1‑100 (MQ=255)
                                                                                    ggCACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAg            >  2:3235997/1‑100 (MQ=255)
                                                                                      cACCGAAAGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGtc          >  2:1470804/1‑100 (MQ=255)
                                                                                       aCCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGtct         >  2:2453052/1‑100 (MQ=255)
                                                                                       aCCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGtct         >  1:3714315/1‑100 (MQ=255)
                                                                                       aCCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGtct         >  1:3714073/1‑100 (MQ=255)
                                                                                        ccGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGtctc        >  2:3122195/1‑100 (MQ=255)
                                                                                         cGAAATTGAGGGTGAATTGCGTCGGCCCGTGTGGTATGTCCAGTGACCTGTGCGGCCCACAGTCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTctct       <  2:3529439/100‑1 (MQ=255)
                                                                                             aGTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACg   >  1:3311735/1‑100 (MQ=255)
                                                                                             aGTGAGGGTGAAAGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACg   >  1:2742481/1‑100 (MQ=255)
                                                                                              gTGAGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGc  >  1:961057/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATGAACTGGCACCGAAAGTGCGGGTGAATGGCGTCGGCCCGTGTGGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGC  >  NC_000913/2671662‑2671855

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