Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,977,942 G→T 6.3% L320M (CTG→ATG) ‡ lysA ← diaminopimelate decarboxylase, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,977,9420GT 6.3% 450.8 / 16.8 127L320M (CTG→ATG) ‡lysAdiaminopimelate decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base G (61/58);  new base T (4/4);  total (65/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

Only 100 of 127 total aligned reads displayed.

CCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAG  >  NC_000913/2977800‑2978087
                                                                                                                                              |                                                                                                                                                  
ccTGCCTTCAGTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCAc                                                                                                                                              >  3:669989/1‑148 (MQ=255)
      ttCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGatat                                                                                                                                        >  4:801517/1‑148 (MQ=255)
            tCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCa                                                                                                                                                >  3:1047214/1‑135 (MQ=255)
            tCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCa                                                                                                                                                <  4:1047214/135‑1 (MQ=255)
                 cAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGatat                                                                                                                                        <  2:775797/137‑1 (MQ=255)
                 cAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGatat                                                                                                                                        >  1:775797/1‑137 (MQ=255)
                             ttCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTAc                                                                                                                       <  3:354939/142‑1 (MQ=255)
                              tCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGTCGTCGGTTTACACCGTTGGCGCGTGTTACAGAGAACGACCATCAGCTGCCAGGGCACTGAGCTGGTGGTAACTAGCGTACATTGCCg                                                                                                                >  4:337580/1‑148 (MQ=255)
                              tCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATAtggt                                                                                                                                     >  2:65183/1‑127 (MQ=255)
                              tCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATAtggt                                                                                                                                     <  1:65183/127‑1 (MQ=255)
                               cAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATAtgg                                                                                                                                      <  3:1123643/125‑1 (MQ=255)
                               cAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATAtgg                                                                                                                                      >  4:1123643/1‑125 (MQ=255)
                                aaCATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCggg                                                                                                              <  2:1127381/148‑1 (MQ=255)
                                aaCATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCggg                                                                                                              >  1:1127381/1‑148 (MQ=255)
                                 aCATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATAtg                                                                                                                                       <  1:478318/122‑1 (MQ=255)
                                 aCATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATAtg                                                                                                                                       >  2:478318/1‑122 (MQ=255)
                                 aCATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGc                                                                                                                  <  2:729003/143‑1 (MQ=255)
                                    tttCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTa                                                                                                                                 >  1:355363/1‑125 (MQ=255)
                                    tttCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTg                                                                                                                   >  3:218985/1‑139 (MQ=255)
                                         cccTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGa                                                                                                     >  2:966511/1‑148 (MQ=255)
                                           cTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATc                                                                                                   >  4:687980/1‑148 (MQ=255)
                                             tCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGt                                                                                                 <  2:283513/148‑1 (MQ=255)
                                             tCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGt                                                                                                 <  1:296089/148‑1 (MQ=255)
                                             tCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGt                                                                                                 <  3:337580/148‑1 (MQ=255)
                                             tCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGt                                                                                                 <  1:56172/148‑1 (MQ=255)
                                              acaaatggggagccgccactttgtgctggttgatgccgggttcaacgatctgatgcgcccggcaatgtacggtagttaccaccatATCAGCTGCCAGGGCACTGata                                                                                                                                         >  4:534826‑M1/86‑107 (MQ=255)
                                               ctgctgGGTTAAGCAATCGCCAGCTTCACATAACGGTCCGGCGACGCCGGGTTCCACCGTTGGGGCGTGTACCAGAGAACGACCATCAGCGGCCAGGGCACTGAAATGCAGGACACTACCGTCCATTGCCGGCCGCATCAGATCGttt                                                                                               >  2:296089/1‑147 (MQ=255)
                                                 gctgGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGaa                                                                                             >  3:186587/1‑148 (MQ=255)
                                                 gctgGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGAGCGTTGaa                                                                                             <  3:801517/148‑1 (MQ=255)
                                                   tgGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATAtggt                                                                                                                                     <  4:495928/106‑1 (MQ=255)
                                                   tgGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATAtggt                                                                                                                                     >  3:495928/1‑106 (MQ=255)
                                                      gTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCgg                                                                                        <  4:669989/148‑1 (MQ=255)
                                                       gagccgccactttgtgctggttgatgccgggttcaacgatctgatgcgcccggcaatgtacggtagttaccaccatATCAGCTGCCAGGGCACTGat                                                                                                                                          >  4:159469‑M1/77‑97 (MQ=255)
                                                        aaaGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAAcc                                                                                           <  4:1026049/143‑1 (MQ=255)
                                                        aaaGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCa                                                                                      >  2:1055539/1‑148 (MQ=255)
                                                         aaGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCAt                                                                                     <  3:687980/148‑1 (MQ=255)
                                                            aCATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCaa                                                                                  >  3:953140/1‑148 (MQ=255)
                                                             cATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAAc                                                                                 <  3:1140111/148‑1 (MQ=255)
                                                              cactttgtgctggttgatgccgggttcaacgatctgatgcgcccggcaatgtacggtagttaccaccatATCAGCTGCCAGGGCACTGatat                                                                                                                                        <  2:59240‑M1/23‑1 (MQ=255)
                                                              cactttgtgctggttgatgccgggttcaacgatctgatgcgcccggcaatgtacggtagttaccaccatATCAGCTGCCAGGGCACTGatat                                                                                                                                        >  1:59240‑M1/70‑92 (MQ=255)
                                                              aTCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCa                                                                                      >  4:271605/1‑142 (MQ=255)
                                                               tCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGt                                                                                                 <  2:213930/130‑1 (MQ=255)
                                                                cGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATAtg                                                                                                                                       <  4:737563/91‑1 (MQ=255)
                                                                        tCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCg                                                                                                                <  1:1067201/106‑1 (MQ=255)
                                                                           caTAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGgcgc                                                                                                           >  2:111173/1‑108 (MQ=255)
                                                                             tAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCAt                                                                                     >  1:945094/1‑128 (MQ=255)
                                                                             tAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCAt                                                                                     <  2:945094/128‑1 (MQ=255)
                                                                              aaCGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGc                                                                                       >  3:49443/1‑125 (MQ=255)
                                                                              aaCGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGc                                                                                       <  4:49443/125‑1 (MQ=255)
                                                                              aaCGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTggcggc                                                                >  1:552491/1‑148 (MQ=255)
                                                                                cGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTc                                                              <  2:552491/148‑1 (MQ=255)
                                                                                cGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTc                                                              >  3:428737/1‑148 (MQ=255)
                                                                                  gTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTccc                                                            <  1:1055539/148‑1 (MQ=255)
                                                                                    ccGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCa                                                          <  4:428737/148‑1 (MQ=255)
                                                                                     cGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATAtgg                                                                                                                                      >  3:1058436/1‑71 (MQ=255)
                                                                                       gCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCgg                                                                                        <  3:719163/115‑1 (MQ=255)
                                                                                        cgacgaCGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCAtttg                                                      >  3:1110145/1‑148 (MQ=255)
                                                                                          acgacgGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACaa                                                                         <  2:108079/127‑1 (MQ=255)
                                                                                          acgacgGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACaa                                                                         >  1:108079/1‑127 (MQ=255)
                                                                                           cgacgGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCAtttgttt                                                   <  1:966511/148‑1 (MQ=255)
                                                                                            gacgGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCAtttgtttg                                                  <  3:1012097/148‑1 (MQ=255)
                                                                                            gacgGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCAtttgtttg                                                  <  4:186587/148‑1 (MQ=255)
                                                                                             acgGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTg                                                                                               <  1:695022/102‑1 (MQ=255)
                                                                                             acgGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCAtttgtttg                                                  <  3:519597/147‑1 (MQ=255)
                                                                                                  ttCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCATCACAAAGTGGCTGCTCCCCATTTGTTTGACGCTc                                            <  4:1110145/148‑1 (MQ=255)
                                                                                                  ttCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGc                                                                             >  3:336544/1‑115 (MQ=255)
                                                                                                  ttCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGc                                                                             <  4:336544/115‑1 (MQ=255)
                                                                                                   tCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTa                                                                                                                             >  4:1048595/1‑66 (MQ=255)
                                                                                                   tCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTa                                                                                                                             <  3:1048595/66‑1 (MQ=255)
                                                                                                      aCCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCa                                        >  4:622753/1‑148 (MQ=255)
                                                                                                      aCCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCa                                        >  2:179606/1‑148 (MQ=255)
                                                                                                       ccGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCg                                                                                                                          >  2:741901/1‑65 (MQ=255)
                                                                                                       ccGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCg                                                                                                                          <  1:741901/65‑1 (MQ=255)
                                                                                                        cGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATc                                                                                    >  3:106791/1‑102 (MQ=255)
                                                                                                        cGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATc                                                                                    <  4:106791/102‑1 (MQ=255)
                                                                                                             gcgcGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACAt                                                                                                                     >  3:214501/1‑64 (MQ=255)
                                                                                                             gcgcGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACAt                                                                                                                     <  4:214501/64‑1 (MQ=255)
                                                                                                              cgcgTGTTCCAGCGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGta                                <  1:179606/148‑1 (MQ=255)
                                                                                                                     tCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATa                         >  2:1125118/1‑148 (MQ=255)
                                                                                                                      ccAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATAc                        >  3:844547/1‑148 (MQ=255)
                                                                                                                       cAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGa                                                                                                     >  1:1030313/1‑70 (MQ=255)
                                                                                                                       cAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGa                                                                                                     <  2:1030313/70‑1 (MQ=255)
                                                                                                                        agagAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCAtttgtttg                                                  <  2:949756/120‑1 (MQ=255)
                                                                                                                         gagaACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCAt                                                                                                         >  2:409406/1‑64 (MQ=255)
                                                                                                                         gagaACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCAt                                                                                                         <  1:409406/64‑1 (MQ=255)
                                                                                                                          agaACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTccc                                                            >  1:138008/1‑108 (MQ=255)
                                                                                                                              accatATCAGTGCCCTGGCAGCTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAg                                                                       <  2:716447‑M1/88‑1 (MQ=255)
                                                                                                                              accatATCAGTGCCCTGGCAGCTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAg                                                                       >  1:716447‑M1/6‑93 (MQ=255)
                                                                                                                                catATCAGTGCCCTGGCAGCTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCaa                                                                                  >  1:822691‑M1/4‑80 (MQ=255)
                                                                                                                                catATCAGTGCCCTGGCAGCTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAGGTGGCGGCTCCCCAtt                                                        >  2:580151‑M1/4‑106 (MQ=255)
                                                                                                                                catATCAGTGCCCTGGCAGCTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCAtttg                                                      <  4:1067184‑M1/105‑1 (MQ=255)
                                                                                                                                aCCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTg                                                                     >  2:554428/1‑93 (MQ=255)
                                                                                                                                aCCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTg                                                                     <  1:554428/93‑1 (MQ=255)
                                                                                                                                     CAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTAc         >  3:622753‑M1/1‑148 (MQ=255)
                                                                                                                                     CAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTAc         >  4:401953‑M1/1‑148 (MQ=255)
                                                                                                                                      AGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTAcc        <  3:845638‑M1/148‑1 (MQ=255)
                                                                                                                                          gCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGa    <  4:845638/148‑1 (MQ=255)
                                                                                                                                           ccAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCAtt                                                        <  4:176734/95‑1 (MQ=255)
                                                                                                                                           ccAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCAtt                                                        >  3:176734/1‑95 (MQ=255)
                                                                                                                                            cAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAg  <  1:1125118/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCACATAACGGTCCGGCGACGACGGTTTCCACCGTTGGCGCGTGTTCCAGAGAACGACCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAG  >  NC_000913/2977800‑2978087

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: