Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,962,276 G→T 7.1% G391V (GGC→GTC)  lysA‑II → diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,962,2760GT 7.1% 205.8 / 3.1 68G391V (GGC→GTC) lysA‑IIdiaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (31/32);  new base T (3/2);  total (34/34)
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.24e-01

CTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTCAGC  >  NC_002947/5962135‑5962422
                                                                                                                                             |                                                                                                                                                   
cTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAg                                                                                                                                              <  6:131081/149‑1 (MQ=255)
 tGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGCCGTCGACCAg                                                                                                                                              >  4:148316/1‑148 (MQ=255)
   gCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGGCCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCg                                                                                                                                           >  1:234763/1‑149 (MQ=255)
   gCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGCCGTCGACCAgga                                                                                                                                            <  4:88257/148‑2 (MQ=255)
   gCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGCCGTCGACCAgga                                                                                                                                            <  2:98186/148‑2 (MQ=255)
   gCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGc                                                                                                                                            <  3:287543/148‑1 (MQ=255)
   gCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGc                                                                                                                                            <  2:3616/148‑1 (MQ=255)
          ccTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGaa                                                                                                                                     <  1:265771/148‑1 (MQ=255)
           cTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAg                                                                                                                                    <  3:294346/148‑1 (MQ=255)
            tCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTa                                                                                                                                  >  4:276970/1‑149 (MQ=255)
                cAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTAcgcc                                                                                                                              >  2:207191/1‑149 (MQ=255)
                      ccGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTAcgccgccgc                                                                                                                         >  4:32049/1‑148 (MQ=255)
                      ccGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTAcgccgccgc                                                                                                                         >  2:187508/1‑148 (MQ=255)
                       cgcgTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTAcgccgccgc                                                                                                                         >  2:46250/1‑147 (MQ=255)
                       cgcgTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTAcgccgccgc                                                                                                                         >  7:71404/1‑147 (MQ=255)
                            ccTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCg                                                                                                                  >  6:9254/1‑149 (MQ=255)
                              tCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACAACCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTc                                                                                                                >  6:15175/1‑149 (MQ=255)
                                aaCCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTc                                                                                                                <  1:207191/147‑1 (MQ=255)
                                aaCCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGATCGAAGTACGCCGCCGCGAGACCGTCGc                                                                                                              <  1:187508/149‑1 (MQ=255)
                                      gCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGTCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGCACt                                                                                                        >  1:139658/1‑149 (MQ=255)
                                        cGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGTCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTg                                                                                                       >  5:81060/1‑148 (MQ=255)
                                        cGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGt                                                                                                      <  3:32049/149‑1 (MQ=255)
                                        cGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCAGCGAGACCGTCGCCGAACTg                                                                                                       >  3:137107/1‑148 (MQ=255)
                                          aaGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTAt                                                                                                    <  1:46250/149‑1 (MQ=255)
                                           aGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATg                                                                                                   >  6:118125/1‑149 (MQ=255)
                                              tGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTg                                                                                                <  2:200828/149‑1 (MQ=255)
                                                  ctgctgGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCACTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCg                                                                                             >  6:141247/1‑148 (MQ=255)
                                                   tgctgGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTg                                                                                                <  5:210602/144‑1 (MQ=255)
                                                   tgctgGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTg                                                                                                >  6:210602/1‑144 (MQ=255)
                                                    gctgGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGa                                                                                            >  5:195447/1‑147 (MQ=255)
                                                    gctgGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGa                                                                                            >  5:68306/1‑147 (MQ=255)
                                                     ctgGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGaa                                                                                           <  5:9254/147‑1 (MQ=255)
                                                                 tCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCg                                                                              >  8:86422/1‑147 (MQ=255)
                                                                 tCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGa                                                                            <  8:71404/149‑1 (MQ=255)
                                                                   cgcgGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGt                                                                          <  2:234763/149‑1 (MQ=255)
                                                                               tGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCAtgctg                                                              <  6:68306/149‑1 (MQ=255)
                                                                                        cATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGttt                                                      >  2:88984/1‑148 (MQ=255)
                                                                                          tGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTa                                                    <  1:45189/149‑1 (MQ=255)
                                                                                          tGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTAc                                                   >  3:194098/1‑149 (MQ=255)
                                                                                            aGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTAcc                                                  <  5:118125/148‑1 (MQ=255)
                                                                                               tCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTAcc                                                  >  2:165083/1‑145 (MQ=255)
                                                                                                 cAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGCCCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATg                                            >  8:114378/1‑149 (MQ=255)
                                                                                                 cAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATg                                            >  3:201898/1‑149 (MQ=255)
                                                                                                  aaCTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGc                                           <  4:201898/149‑1 (MQ=255)
                                                                                                     tACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTgg                                                                                               >  1:264845/1‑95 (MQ=255)
                                                                                                     tACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTgg                                                                                               <  2:264845/95‑1 (MQ=255)
                                                                                                       caacaCCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtg                                        <  5:15175/147‑1 (MQ=255)
                                                                                                        aacaCCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggc                                     >  6:332536/1‑149 (MQ=255)
                                                                                                         acacCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggct                                    >  4:35960/1‑149 (MQ=255)
                                                                                                           acCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggctgg                                  >  5:296114/1‑149 (MQ=255)
                                                                                                            cccGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggctggg                                 <  4:194098/149‑1 (MQ=255)
                                                                                                             ccGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggctgg                                  <  8:98458/147‑1 (MQ=255)
                                                                                                             ccGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggctgg                                  >  7:98458/1‑147 (MQ=255)
                                                                                                                  gCCGTTGCGCTTAAATCCTGGTCTACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGa                           <  4:148543/149‑1 (MQ=255)
                                                                                                                      ttGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTc                       <  3:276970/149‑1 (MQ=255)
                                                                                                                         cgcTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCAt                     <  4:137107/148‑1 (MQ=255)
                                                                                                                         cgcTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATg                    >  4:767/1‑149 (MQ=255)
                                                                                                                         cgcTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATg                    <  1:88984/149‑1 (MQ=255)
                                                                                                                             gAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTcc                >  1:85413/1‑149 (MQ=255)
                                                                                                                               aaTCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTg              >  5:242150/1‑149 (MQ=255)
                                                                                                                               aaTCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTg              <  6:242150/149‑1 (MQ=255)
                                                                                                                               aaTCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTg              <  5:332536/149‑1 (MQ=255)
                                                                                                                                 tCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAAc                             >  2:328952/1‑132 (MQ=255)
                                                                                                                                 tCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAAc                             <  1:328952/132‑1 (MQ=255)
                                                                                                                                     ggTCGCCGTCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTg        >  4:302251/1‑149 (MQ=255)
                                                                                                                                         gACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTca    <  3:35960/149‑1 (MQ=255)
                                                                                                                                          aCGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTcag   <  3:767/149‑1 (MQ=255)
                                                                                                                                           cGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTcagc  <  5:141247/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAG‑CCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTCAGC  >  NC_002947/5962135‑5962422

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: ‑