Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,836,905 G→C 5.6% S278T (AGC→ACC)  ynjB → putative ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,836,9050GC 5.6% 364.1 / 8.6 103S278T (AGC→ACC) ynjBputative ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (52/45);  new base C (2/4);  total (54/49)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

AATATCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAACTG  >  NC_000913/1836758‑1837053
                                                                                                                                                   |                                                                                                                                                    
aaTATCTTGATGTGCTGCATCCGTATTCGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAgc                                                                                                                                                     >  2:27526/1‑149 (MQ=255)
   atgttGATGTGCTGCATACGTATATGTTGCGCGAAGGCAAGGAGTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCGGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAgcgg                                                                                                                                                   >  4:561778/4‑148 (MQ=255)
     cTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAgcggcg                                                                                                                                                 >  3:656088/1‑148 (MQ=255)
     cTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGa                                                                                                                                                >  1:714233/1‑149 (MQ=255)
     cTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGa                                                                                                                                                <  1:130229/149‑1 (MQ=255)
         aTGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCGCTGGCGAtttt                                                                                                                                            <  1:202707/149‑2 (MQ=255)
         aTGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCGCTGGCGAttt                                                                                                                                             <  3:65972/148‑1 (MQ=255)
          tgtgCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGc                                                                                                                                           >  2:996074/1‑149 (MQ=255)
            tgctgcATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCt                                                                                                                                         <  1:28995/149‑1 (MQ=255)
            tgctgcATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCt                                                                                                                                         <  2:122902/149‑1 (MQ=255)
             ggtgcATCCGTATTTGTGGCACGCAGGATAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTACGCCTGTCGCTGACCTTTATCCCCGCACATGCGCAGCTAAAAATCGCCAGCGGCGATTTGCCt                                                                                                                                         >  4:657699/3‑148 (MQ=255)
                  aTCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGc                                                                                                                                   >  1:404001/1‑149 (MQ=255)
                   tCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGc                                                                                                                                   <  4:894120/148‑1 (MQ=255)
                    ccGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCa                                                                                                                                  >  3:129579/1‑148 (MQ=255)
                      gTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGtt                                                                                                                               >  2:607742/1‑149 (MQ=255)
                      gTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGtt                                                                                                                               >  2:208611/1‑149 (MQ=255)
                       tATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGtt                                                                                                                               >  3:1010964/1‑148 (MQ=255)
                         tgtgtgGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTAc                                                                                                                             <  4:429970/146‑1 (MQ=255)
                            gtgGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGtt                                                                                                                         >  2:524500/1‑149 (MQ=255)
                            gtgGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGtt                                                                                                                         >  1:589158/1‑149 (MQ=255)
                              ggCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGttt                                                                                                                        <  4:1008546/148‑1 (MQ=255)
                              ggCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGttt                                                                                                                        <  4:701001/148‑1 (MQ=255)
                                cgcgAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTg                                                                                                                      <  3:524825/148‑1 (MQ=255)
                                      ggaaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGAt                                                                                                                                               >  2:231286/1‑117 (MQ=255)
                                      ggaaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGAt                                                                                                                                               <  1:231286/117‑1 (MQ=255)
                                       gaaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTcgc                                                                                                              >  2:720478/1‑149 (MQ=255)
                                       gaaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTcgc                                                                                                              >  2:909555/1‑149 (MQ=255)
                                       gaaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTcgc                                                                                                              >  1:257774/1‑149 (MQ=255)
                                        aaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTTGCTTTcgc                                                                                                              >  3:911778/1‑148 (MQ=255)
                                        aaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTcgc                                                                                                              >  4:367346/1‑148 (MQ=255)
                                        aaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTcgc                                                                                                              >  3:197151/1‑148 (MQ=255)
                                         aaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTcgcg                                                                                                             <  4:1028183/148‑1 (MQ=255)
                                         aaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGa                                                                                                            <  2:462695/149‑1 (MQ=255)
                                          aggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGa                                                                                                            >  4:81770/1‑148 (MQ=255)
                                          aggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAg                                                                                                           <  2:757791/149‑1 (MQ=255)
                                           ggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAgg                                                                                                          >  2:205187/1‑149 (MQ=255)
                                              tttCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGga                                                                                                       >  1:685375/1‑149 (MQ=255)
                                                  ccGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGgatgat                                                                                                   <  2:737958/149‑1 (MQ=255)
                                                  ccGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGgatgat                                                                                                   >  2:887216/1‑149 (MQ=255)
                                                  ccGCCTTCACCCGCGCGGATGGATGCTCTGCCGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAggg                                                                                                         >  2:601050/1‑143 (MQ=255)
                                                   cGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGgatgat                                                                                                   <  4:206579/148‑1 (MQ=255)
                                                   cGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGAtt                                                                                                  <  2:872826/149‑1 (MQ=255)
                                                   cGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGAtt                                                                                                  <  2:103153/149‑1 (MQ=255)
                                                   cGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGAtt                                                                                                  <  2:326932/149‑1 (MQ=255)
                                                    gCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGAtt                                                                                                  <  4:196577/148‑1 (MQ=255)
                                                                gcgGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAATCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCAt                                                                                      <  4:898614/148‑1 (MQ=255)
                                                                       gttgCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCGGTCACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCa                                                                               >  4:939635/3‑148 (MQ=255)
                                                                       gatgCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCa                                                                               <  4:295240/148‑1 (MQ=255)
                                                                       gatgCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCAc                                                                              <  1:27526/149‑1 (MQ=255)
                                                                          gCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCAttt                                                                                    >  3:16266/1‑140 (MQ=255)
                                                                          gCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCAttt                                                                                    <  4:16266/140‑1 (MQ=255)
                                                                            tctGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCAtt                                                                          >  3:399235/1‑148 (MQ=255)
                                                                              tgctgAAAGCCGGCACATTGCGCCTTTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCt                                                                       <  2:817921/149‑1 (MQ=255)
                                                                                ctgAAAGCCGGCACATTGCGCCTGTCGGTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTg                                                                      >  3:907463/1‑148 (MQ=255)
                                                                                  gAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGcc                                                                    >  4:1039206/1‑148 (MQ=255)
                                                                                  gAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGcc                                                                    >  4:139743/1‑148 (MQ=255)
                                                                                   aaaGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCaa                                                                  >  2:378000/1‑149 (MQ=255)
                                                                                   aaaGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCa                                                                   <  4:555699/148‑1 (MQ=255)
                                                                                   aaaGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATCTCGTCACCATTCCTGCCa                                                                   >  3:137161/1‑148 (MQ=255)
                                                                                      gCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAAcg                                                                <  4:925400/148‑1 (MQ=255)
                                                                                      gCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATCGGCAACGTGCATTTCGTCACCATTCCTGCCAAcg                                                                >  4:14368/1‑148 (MQ=255)
                                                                                          gCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGTCAACGCGAAt                                                           <  1:107200/149‑1 (MQ=255)
                                                                                          gCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAAt                                                           <  1:996074/149‑1 (MQ=255)
                                                                                            acaTTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGc                                                         <  1:607742/149‑1 (MQ=255)
                                                                                              aTTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGcc                                                        <  3:561778/148‑1 (MQ=255)
                                                                                               ttGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAAcgc                                                               <  1:436259/140‑1 (MQ=255)
                                                                                               ttGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAAcgc                                                               >  2:436259/1‑140 (MQ=255)
                                                                                               ttGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATg                                                          <  4:922856/145‑1 (MQ=255)
                                                                                               ttGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATg                                                          >  3:922856/1‑145 (MQ=255)
                                                                                                tGCGCCTATCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCAcc                                                                             >  4:92756/1‑125 (MQ=255)
                                                                                                    ccTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGtgctg                                                 <  1:887216/149‑1 (MQ=255)
                                                                                                     cTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGtgctgc                                                >  2:862568/1‑149 (MQ=255)
                                                                                                     cTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGtgctgc                                                <  1:1066796/149‑1 (MQ=255)
                                                                                                       gTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGTAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTc                                                                                <  3:906413/115‑1 (MQ=255)
                                                                                                       gTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGTAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTc                                                                                >  4:906413/1‑115 (MQ=255)
                                                                                                        tCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTgcgg                                              >  3:508444/1‑148 (MQ=255)
                                                                                                          gCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTgcggcg                                            >  3:72076/1‑148 (MQ=255)
                                                                                                           cTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCAAATGCCAGTGCTGCGGCGaa                                          <  2:404001/149‑1 (MQ=255)
                                                                                                           cTGACCTTAAACCCCGCGCATGCACAGCAAATAATCGCCACCTGCGATTTGCCTGCAAGCATTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGAAGTTCGACGCCATCATTGCCACCGCAAATGCGAGTGCCGAGGCCaa                                          <  2:382004/149‑1 (MQ=255)
                                                                                                            tGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAg                                         >  1:264628/1‑149 (MQ=255)
                                                                                                             gACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAg                                         >  4:878581/1‑148 (MQ=255)
                                                                                                             gACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAg                                         >  3:210733/1‑148 (MQ=255)
                                                                                                               ccTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGt                                       >  4:136655/1‑148 (MQ=255)
                                                                                                                cTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAg                                     >  1:1042483/1‑149 (MQ=255)
                                                                                                                 tttAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGt                                    <  1:705603/149‑1 (MQ=255)
                                                                                                                 tttAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGt                                    >  2:184780/1‑149 (MQ=255)
                                                                                                                       cccGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCaa                              <  2:783334/149‑1 (MQ=255)
                                                                                                                       cccGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCa                               <  4:928906/148‑1 (MQ=255)
                                                                                                                        ccGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGa                                                             <  4:546305/117‑1 (MQ=255)
                                                                                                                        ccGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGa                                                             >  3:546305/1‑117 (MQ=255)
                                                                                                                        ccGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGtg                                                    <  1:660003/126‑1 (MQ=255)
                                                                                                                        ccGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGtg                                                    >  2:660003/1‑126 (MQ=255)
                                                                                                                        ccGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCaa                              >  3:164815/1‑148 (MQ=255)
                                                                                                                        ccGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCaa                              >  3:632759/1‑148 (MQ=255)
                                                                                                                                        aaaaaTCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCAccc              >  4:311586/1‑148 (MQ=255)
                                                                                                                                         caaaTCGCCGCTGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCa                               <  3:865227/129‑1 (MQ=255)
                                                                                                                                          aaaTCGCCGCTGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGAt           >  1:33612/1‑149 (MQ=255)
                                                                                                                                          aaaTCGCCGCTGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGAt           >  1:943376/1‑149 (MQ=255)
                                                                                                                                                 cAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAAc    <  1:208611/149‑1 (MQ=255)
                                                                                                                                                   gcggcgATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAACTg  <  2:34531/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AATATCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAACTG  >  NC_000913/1836758‑1837053

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