Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,678,903 T→C 5.4% D159D (GAT→GAC)  ydgH → DUF1471 family periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,678,9030TC 5.4% 240.7 / 3.1 76D159D (GAT→GAC) ydgHDUF1471 family periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base C (2/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

TGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGA  >  NC_000913/1678755‑1679048
                                                                                                                                                    |                                                                                                                                                 
tGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGAt                                                                                                                                                   >  7:69305/1‑149 (MQ=255)
    gAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCa                                                                                                                                               >  2:391309/1‑149 (MQ=255)
     aaaGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCaa                                                                                                                                              >  1:243867/1‑149 (MQ=255)
      aaGATCAGGCTGTGCTGATTGAACCGTTTGACACGGGCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAAc                                                                                                                                             <  6:384578/149‑1 (MQ=255)
        gATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCa                                                                                                                                           <  8:69305/149‑1 (MQ=255)
            aGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGggc                                                                                                                                       >  7:238234/1‑149 (MQ=255)
               cTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGggcggc                                                                                                                                    <  1:137611/149‑1 (MQ=255)
                  ttCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCa                                                                                                                                   <  2:243867/147‑1 (MQ=255)
                       aTTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCa                                                                                                                               <  3:218618/146‑1 (MQ=255)
                       aTTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGc                                                                                                                             <  8:241213/148‑1 (MQ=255)
                          gAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTa                                                                                                                          <  6:236168/148‑1 (MQ=255)
                          gAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTAt                                                                                                                         >  5:397520/1‑149 (MQ=255)
                           aaCCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTAt                                                                                                                         >  4:374577/1‑148 (MQ=255)
                           aaCCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTAt                                                                                                                         >  2:356350/1‑148 (MQ=255)
                            aCCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTa                                                                                                                       <  7:134886/149‑1 (MQ=255)
                             ccGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTAc                                                                                                                      >  8:328962/1‑149 (MQ=255)
                              cGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACt                                                                                                                     >  1:326670/1‑149 (MQ=255)
                              cGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACt                                                                                                                     <  3:272309/149‑1 (MQ=255)
                                tttGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGc                                                                                                                   <  1:391309/149‑1 (MQ=255)
                                   gACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCAtt                                                                                                                >  1:213057/1‑149 (MQ=255)
                                    acacGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCAtt                                                                                                                >  1:39938/1‑148 (MQ=255)
                                    acacGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCAtt                                                                                                                >  7:400177/1‑148 (MQ=255)
                                    acacGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACAGCAtt                                                                                                                >  8:242364/1‑148 (MQ=255)
                                    acacGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATAGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTc                                                                                                               <  8:238234/149‑1 (MQ=255)
                                          tCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCtat                                                                                                         <  3:374577/149‑1 (MQ=255)
                                            aCCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATaa                                                                                                       <  2:326670/149‑1 (MQ=255)
                                                tCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATaaaa                                                                                                     >  1:49947/1‑147 (MQ=255)
                                                   aGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGTGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGAt                                                                                                >  1:392110/1‑149 (MQ=255)
                                                    gggCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATg                                                                                               >  1:363554/1‑149 (MQ=255)
                                                    gggCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATg                                                                                               <  1:356350/149‑1 (MQ=255)
                                                     ggCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGc                                                                                              >  7:133816/1‑149 (MQ=255)
                                                      gCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCt                                                                                             <  7:328962/149‑1 (MQ=255)
                                                        ttCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAgg                                                                                                                                         >  7:168333/1‑103 (MQ=255)
                                                        ttCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAgg                                                                                                                                         <  8:168333/103‑1 (MQ=255)
                                                             tagTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTTCATTCATCTATAAAAAAGATGCTAAGaaa                                                                                       <  8:183284/146‑1 (MQ=255)
                                                             tCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTAt                                                                                                                         >  8:389613/1‑114 (MQ=255)
                                                             tCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTAt                                                                                                                         <  7:389613/114‑1 (MQ=255)
                                                                tAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTa                                                                                   <  2:363554/149‑1 (MQ=255)
                                                                  gccagccaGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATc                                                                                 <  4:327042/149‑1 (MQ=255)
                                                                   ccagccagAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGTCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGATAAGAACCGTATCg                                                                                >  2:321354/1‑149 (MQ=255)
                                                                     agccagAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTc                                                                              >  1:255600/1‑149 (MQ=255)
                                                                        cagAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGACGCCAACCAGGGCGGCAACCAGc                                                                                                                             <  4:89582/99‑1 (MQ=255)
                                                                        cagAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGACGCCAACCAGGGCGGCAACCAGc                                                                                                                             >  3:89582/1‑99 (MQ=255)
                                                                                aaTGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGa                                                                   >  4:383281/1‑149 (MQ=255)
                                                                                aaTGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGACGCCAACCAGGGCGGCAACCAGCGTa                                                                                                                          <  5:108211/94‑1 (MQ=255)
                                                                                aaTGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGACGCCAACCAGGGCGGCAACCAGCGTa                                                                                                                          >  6:108211/1‑94 (MQ=255)
                                                                                 atgatgCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATaa                                                                                                       >  1:180104/1‑112 (MQ=255)
                                                                                 atgatgCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATaa                                                                                                       <  2:180104/112‑1 (MQ=255)
                                                                                         aTCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAg                                                                                                  >  3:103287/1‑109 (MQ=255)
                                                                                         aTCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAg                                                                                                  <  4:103287/109‑1 (MQ=255)
                                                                                            aCCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGc                                                       <  3:383281/149‑1 (MQ=255)
                                                                                            aCCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGc                                                       <  1:63867/149‑1 (MQ=255)
                                                                                            aCCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGc                                                       <  2:392110/149‑1 (MQ=255)
                                                                                               aaaGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGAtgt                                                                >  4:22569/1‑137 (MQ=255)
                                                                                               aaaGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGAtgt                                                                <  3:22569/137‑1 (MQ=255)
                                                                                                 aGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTAACGTCCAGAGCCCGGATGTGATCCCGGCAGAtt                                                  >  2:51284/1‑149 (MQ=255)
                                                                                                  gCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACAGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTc                                                 <  7:242364/149‑1 (MQ=255)
                                                                                                     gCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGAt                                                                  <  4:26636/129‑1 (MQ=255)
                                                                                                     gCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGAt                                                                  >  3:26636/1‑129 (MQ=255)
                                                                                                           gCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTCCGAAGCAgg                                        >  5:117233/1‑149 (MQ=255)
                                                                                                            cGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGa                                       >  2:373918/1‑149 (MQ=255)
                                                                                                                  gTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTg                                   >  8:64798/1‑147 (MQ=255)
                                                                                                                  gTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGACAGATTCCGAAGCAGGACGTGCa                                 <  8:133816/149‑1 (MQ=255)
                                                                                                                    gCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGc                               >  2:357608/1‑149 (MQ=255)
                                                                                                                    gCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAAACAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGc                               >  5:46147/1‑149 (MQ=255)
                                                                                                                        aCTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTg                           >  6:328734/1‑149 (MQ=255)
                                                                                                                         ctctTTCTACATCGTTCTTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTgg                          <  1:357608/149‑1 (MQ=255)
                                                                                                                         ctctTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTg                           <  7:268898/148‑1 (MQ=255)
                                                                                                                            tttCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTg                                   >  4:159603/1‑137 (MQ=255)
                                                                                                                            tttCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTg                                   <  3:159603/137‑1 (MQ=255)
                                                                                                                            tttCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGc                         >  5:219073/1‑147 (MQ=255)
                                                                                                                              tCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGc                         <  6:219073/145‑1 (MQ=255)
                                                                                                                                         cgtCAACTCGAAGCCAACCAGGGCGGCACCCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCGCTGGCTGCCGGTGGCGAAGc          <  6:397520/149‑1 (MQ=255)
                                                                                                                                             aaaTCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCcgc       >  4:335235/1‑148 (MQ=255)
                                                                                                                                              aaTCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAAc                                                                                      <  5:111420/68‑1 (MQ=255)
                                                                                                                                              aaTCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAAc                                                                                      >  6:111420/1‑68 (MQ=255)
                                                                                                                                                  gATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCgaaga  >  5:349649/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGCCGAAAGATCAGGCTGTTCTGATTGAACCGTTTGACACGGTCACCGTCCAGGGCTTCTATCGTAGCCAGCCAGAAGTCAATGATGCCATCACCAAAGCGGCAAAAGCGAAAGGTGCCTACTCTTTCTACATCGTTCGTCAAATCGATGCCAACCAGGGCGGCAACCAGCGTATTACTGCATTCATCTATAAAAAAGATGCTAAGAAACGTATCGTCCAGAGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGA  >  NC_000913/1678755‑1679048

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