Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 165,279 G→A 5.9% V184I (GTC→ATC)  mrcB → fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913165,2790GA 5.9% 215.0 / 2.6 68V184I (GTC→ATC) mrcBfused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base G (35/29);  new base A (2/2);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

GTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTGAGATGATTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGTGGTTTCCCGGATTTGCTGGTGGATACTTTGCTGGC  >  NC_000913/165135‑165421
                                                                                                                                                |                                                                                                                                              
gtgtCGAAAATGACCCGTCCTGGCGATTTTACCGTGCAGGCCAACAGCATTGAGATGATTCGCCGTCCGTTTAATTTCCCGGACAGTAAAGAAGGAGAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCaa                                                                                                                                            >  1:202494/1‑149 (MQ=255)
gtgtCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTGAGATGATTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCaa                                                                                                                                            >  5:123635/1‑149 (MQ=255)
      aaaaTGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTGAGATGATTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGga                                                                                                                                      <  6:339280/149‑1 (MQ=255)
            aCCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTGAGATGATTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGaacaa                                                                                                                                <  1:472011/149‑1 (MQ=255)
                gTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTGAGATGATTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGt                                                                                                                            <  4:226427/149‑1 (MQ=255)
                   cTGGCGAATTTACCGTGCAGGCCAACAGCATTGAGATGATTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCa                                                                                                                          <  2:374924/148‑1 (MQ=255)
                     ggCGAATTTACCGTGCAGGCCAACAGCATTGAGATGATTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGtt                                                                                                                       >  4:389390/1‑149 (MQ=255)
                     ggCGAATTTACCGTGCAGGCCAACAGCATTGAGATGATTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCATCAATATGGAGAACAACCGTCAGtt                                                                                                                       >  8:362816/1‑149 (MQ=255)
                     ggCGAATTTACCGTGCAGGCCAACAGCATTGAGATGATTCGCCGACCGTTAGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCAACAGCCGACGAACATCAATATGGAGAACAACCGTCAGtt                                                                                                                       >  8:356415/1‑149 (MQ=255)
                          aTTTACCGTGCAGGCCAACAGCATTGAGATGATTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGtt                                                                                                                  <  5:457506/149‑1 (MQ=255)
                              aCCGTGCAGGCCAACAGCATTGAGATGATTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGtt                                                                                                                  <  4:353064/145‑1 (MQ=255)
                              aCCGTGCAGGCCAACAGCATTGAGATGATTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGtt                                                                                                                  >  3:353064/1‑145 (MQ=255)
                              aCCGTGCAGGCCAACAGCATTGAGATGATTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTttct                                                                                                               <  6:123635/148‑1 (MQ=255)
                               ccGTGCAGGCCAACAGCATTGAGATGATTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTttctt                                                                                                              <  7:437466/148‑1 (MQ=255)
                               ccGTGCAGGCCAACAGCATTGAGATGATTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTttctt                                                                                                              <  5:61698/148‑1 (MQ=255)
                                    cAGGCCAACAGCATTGAGATGATTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTc                                                                                                         <  5:60610/148‑1 (MQ=255)
                                                aTTGAGATGATTCGCCGTCCATTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCa                                                                                                                                             <  2:13231/100‑1 (MQ=255)
                                                aTTGAGATGATTCGCCGTCCATTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCa                                                                                                                                             >  1:13231/1‑100 (MQ=255)
                                                   gagaTGATTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTGGATCCGCGTCTGAt                                                                                         >  5:221683/1‑149 (MQ=255)
                                                      atgatTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCAc                                                                                      >  8:228107/1‑149 (MQ=255)
                                                      atgatTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCATCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCAc                                                                                      <  7:362816/149‑1 (MQ=255)
                                                      atgatTCGCCGTCCGTATGATTACCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGAACATCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCAc                                                                                      <  7:356415/149‑1 (MQ=255)
                                                        gatTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCgg                                                                                                                    <  5:365716/117‑1 (MQ=255)
                                                        gatTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCgg                                                                                                                    >  6:365716/1‑117 (MQ=255)
                                                         atTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCa                                                                                    <  3:398876/148‑1 (MQ=255)
                                                          ttCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGaa                                                                                                                                   >  6:119546/1‑100 (MQ=255)
                                                          ttCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGaa                                                                                                                                   <  5:119546/100‑1 (MQ=255)
                                                             gCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGAtc                                                                               >  8:97062/1‑149 (MQ=255)
                                                                gtccgtTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATctct                                                                            >  3:205416/1‑149 (MQ=255)
                                                                 tccgtTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCg                                                                                                           <  7:239010/117‑1 (MQ=255)
                                                                 tccgtTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCg                                                                                                           >  8:239010/1‑117 (MQ=255)
                                                                   cgtTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCg                                                                         >  2:219563/1‑149 (MQ=255)
                                                                    gtTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGc                                                                        >  5:86708/1‑149 (MQ=255)
                                                                     tttGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCa                                                                                                                          <  8:97538/98‑1 (MQ=255)
                                                                     tttGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCa                                                                                                                          >  7:97538/1‑98 (MQ=255)
                                                                     tttGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGAt                                                                                <  4:141143/140‑1 (MQ=255)
                                                                     tttGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGAt                                                                                >  3:141143/1‑140 (MQ=255)
                                                                         aTTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAAc                                                                   >  1:242906/1‑149 (MQ=255)
                                                                           ttCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACgg                                                                 >  4:139560/1‑149 (MQ=255)
                                                                             cccGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGt                                                                <  7:228107/148‑1 (MQ=255)
                                                                              ccGGACAGTAAAGAAGGCCAGGTGAGCGCGCTTCTGAACTTTGATGGCGCTCATCTGGCGCCGATCGTCAATATGGACAACAACCGTCAGTTAGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCAAGATCTCTTCGCCAAACGGt                                                                >  3:265064/1‑147 (MQ=255)
                                                                              ccGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACATCCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCtt                                                                           <  2:316779/136‑1 (MQ=255)
                                                                              ccGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCtt                                                                           >  1:316779/1‑136 (MQ=255)
                                                                              ccGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGa                                                              <  3:139560/149‑1 (MQ=255)
                                                                                    aGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGagcagc                                                         <  6:86708/148‑1 (MQ=255)
                                                                                    aGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCg                                                        >  6:153604/1‑149 (MQ=255)
                                                                                      tAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTc                                                      >  8:312962/1‑149 (MQ=255)
                                                                                        aagaagGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTg                                                    >  4:416278/1‑149 (MQ=255)
                                                                                           aagGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGttt                                                 <  7:312962/149‑1 (MQ=255)
                                                                                                      cgcgcgcgTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCa                                       <  3:389390/148‑1 (MQ=255)
                                                                                                       gcgcgcgTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGt                                     >  3:114387/1‑149 (MQ=255)
                                                                                                        cgcgcgTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGTg                                    <  4:205416/149‑1 (MQ=255)
                                                                                                        cgcgcgTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGTg                                    <  2:202494/149‑1 (MQ=255)
                                                                                                           gcgTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGTGGtt                                 >  5:11182/1‑149 (MQ=255)
                                                                                                                tGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGTGGTTTCCCg                            >  4:402503/1‑149 (MQ=255)
                                                                                                                  aCCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGTGGTTTccc                             <  1:219563/146‑1 (MQ=255)
                                                                                                                    cTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGGTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCTCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGTGGTTTCCCGGAtt                        >  1:391890/1‑149 (MQ=255)
                                                                                                                     tttGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGTGGTTTCCCGGAttt                       >  6:26478/1‑149 (MQ=255)
                                                                                                                                tcatcTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGTGGTTTCCCGGATTTGCTGGTGGATa            >  7:239163/1‑149 (MQ=255)
                                                                                                                                tcatcTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGTGGTTTCCCGGATTTGCTGGTGGATa            >  8:335863/1‑149 (MQ=255)
                                                                                                                                     tGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGt                                                                <  3:22874/92‑1 (MQ=255)
                                                                                                                                     tGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATCATCTCTTCGCAAAACGGt                                                                >  4:22874/1‑92 (MQ=255)
                                                                                                                                       gCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGTGGTTTCCCGGATTTGCt                    >  6:143022/1‑134 (MQ=255)
                                                                                                                                       gCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGTGGTTTCCCGGATTTGCt                    <  5:143022/134‑1 (MQ=255)
                                                                                                                                         gacgaTCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGTGGTTTCCCGGATTTGCTGGTGGATACTTTGCTgg   >  6:343318/1‑149 (MQ=255)
                                                                                                                                          acgaTCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGTGGTTTCCCGGATTTGCTGGTGGATACTTTGCTGGc  >  1:444794/1‑149 (MQ=255)
                                                                                                                                            gaTCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGTGGTTTCCCGGATTTg                      <  3:199966/127‑1 (MQ=255)
                                                                                                                                            gaTCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGTGGTTTCCCGGATTTg                      >  4:199966/1‑127 (MQ=255)
                                                                                                                                                |                                                                                                                                              
GTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTGAGATGATTCGCCGTCCGTTTGATTTCCCGGACAGTAAAGAAGGACAGGTGCGCGCGCGTCTGACCTTTGATGGCGATCATCTGGCGACGATCGTCAATATGGAGAACAACCGTCAGTTCGGTTTCTTCCGTCTTGATCCGCGTCTGATCACCATGATCTCTTCGCCAAACGGTGAGCAGCGTCTGTTTGTGCCGCGCAGTGGTTTCCCGGATTTGCTGGTGGATACTTTGCTGGC  >  NC_000913/165135‑165421

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