Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,031,455 C→A 6.9% L83M (CTG→ATG)  dadA‑I → D‑amino acid dehydrogenase small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,031,4550CA 6.9% 194.5 / 4.1 67L83M (CTG→ATG) dadA‑ID‑amino acid dehydrogenase small subunit
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base A (3/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

GGCCCCGCGATGGAAACAAGCTTCGCAAACGCAGGGCAGGTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGAATATAGCCGAGACTGCCTGGATGAACTCAGGGCTGAGACTGGTATA  >  NC_002947/5031308‑5031592
                                                                                                                                                   |                                                                                                                                         
ggCCCCGCGATGGAAACAAGCTTCGCAAACGCAGGGCAGGTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCt                                                                                                                                          >  4:273914/1‑149 (MQ=255)
  ccccGCGATGGAAACAAGCTTCGCAAACGCAGGGCAGGTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTg                                                                                                                                         >  5:250722/1‑148 (MQ=255)
   cccGCGATGGAAACAAGCTTCGCAAACGCAGGGCAGGTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGa                                                                                                                                       >  7:343742/1‑149 (MQ=255)
      gcgATGGAAACAAGCTTCGCAAACGCAGGGCAGGTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCa                                                                                                                                    <  5:336254/149‑1 (MQ=255)
          tGGAAACAAGCTTCGCAAACGCAGGGCAGGTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTa                                                                                                                                 >  5:201774/1‑148 (MQ=255)
            gAAACAAGCTTCGCAAACGCAGGGCAGGTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATc                                                                                                                               <  3:104726/148‑1 (MQ=255)
            gAAAAAATCTTTTCAAACGCATGGCATGTTTCTCCGGGTTATGCCTCACCTTGGTCAGCCCCTTGAGTGCGATTGAAGGCCATTAGATGGATGCTTCAGCAACACTCACCTTTAGCGATCAAGGCAACGTCTAGCCTGGATCAGTATc                                                                                                                               <  5:157701/148‑1 (MQ=255)
              aaCAAGCTTCGCAAACGCAGGGCAGGTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGtt                                                                                                                            >  5:300960/1‑149 (MQ=255)
              aaCAAGCTTCGCAAACGCAGGGCAGGTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGtt                                                                                                                            >  7:184199/1‑149 (MQ=255)
                 aaGCTTCGCAAACGCAGGGCAGGTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGga                                                                                                                         <  6:250722/149‑1 (MQ=255)
                 aaGCTTCGCAAACGCAGGGCAGGTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGga                                                                                                                         >  5:100670/1‑149 (MQ=255)
                      tCGCAAACGCAGGGCAGGATTCTCCGGGGTATGCCTCCCCTTGGTCAGACCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCCCCTTTAGCGATCAAGCCAACGTCTAGCCTGGCTCAGTATCGTTGGATGATc                                                                                                                    >  6:83592/1‑149 (MQ=255)
                       cGCAAACGCAGGGCAGGTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCAATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAATCAAACGTCTAGCCTGGATCAGTATCGTTCgatgata                                                                                                                    >  8:275212/1‑147 (MQ=255)
                        gCAAACGCAGGGCAGGTTTCTCCTGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCGAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCa                                                                                                                  <  5:18072/149‑1 (MQ=255)
                        gCAAACGCAGGGCAGGTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCa                                                                                                                  <  6:155047/149‑1 (MQ=255)
                              gcagggcaggTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGtt                                                                                                            <  7:162021/149‑1 (MQ=255)
                              gcagggcaggTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGtt                                                                                                            >  7:77743/1‑149 (MQ=255)
                                    caggTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTaa                                                                                                      <  8:310202/149‑1 (MQ=255)
                                     aggTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAAc                                                                                                     <  6:300960/149‑1 (MQ=255)
                                      ggTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACt                                                                                                    <  8:77743/149‑1 (MQ=255)
                                         ttCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGt                                                                                                  >  7:105816/1‑148 (MQ=255)
                                            tcCGTGGAATGCCTACCCTTGGACAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTAAGCAACACGCACCTATAGCGATCAACCCAACGACTACCCTGGCTCAGTATCGTGGCATCATCCACATGTTCCGTCACTGTACTc                                                                                              >  6:302808/1‑149 (MQ=255)
                                            tcCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATc                                                                                              >  8:117777/1‑149 (MQ=255)
                                              cGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCg                                                                                            >  2:164923/1‑149 (MQ=255)
                                                 ggTATGGCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCATCGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCa                                                                                         <  7:275212/149‑1 (MQ=255)
                                                     tGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGc                                                                                      <  3:273914/148‑1 (MQ=255)
                                                           cccTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGc                                                                                 <  4:373828/147‑1 (MQ=255)
                                                            ccTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGt                                                                              >  5:217495/1‑149 (MQ=255)
                                                                 gTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCATGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATg                                                                                  >  4:89502/1‑140 (MQ=255)
                                                                 gTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCATGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATg                                                                                  <  3:89502/140‑1 (MQ=255)
                                                                 gTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCATGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATg                                                                                  <  1:190382/140‑1 (MQ=255)
                                                                 gTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCATGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATg                                                                                  >  2:190382/1‑140 (MQ=255)
                                                                 gTCAGCCCCGGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTcaaca                                                                         <  8:333425/149‑1 (MQ=255)
                                                                   cAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACCCGCACCTTTAGCGCTCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGCTGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTcaaca                                                                         >  8:264687/1‑147 (MQ=255)
                                                                        ccTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACg                                                                  >  3:189714/1‑149 (MQ=255)
                                                                          tGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGa                                                                <  6:201774/149‑1 (MQ=255)
                                                                            tcgTGCCATTTCAGGCCAGTAGATGTATGCTTCATCAACACGCACCTTTCGCTATCAACCCAACGTCTAGCCTGGATCAGTATCTTTGGATGATCCAGATGTTGCGTAACTGTAATACGGCAAGCTATGCAGTCAAAAAAGAACGGATg                                                              <  4:252231/147‑1 (MQ=255)
                                                                                  cATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCg                                                         >  1:87398/1‑148 (MQ=255)
                                                                                   aTTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGaa                                                       >  8:200485/1‑149 (MQ=255)
                                                                                      gAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAAtt                                                     >  5:393177/1‑148 (MQ=255)
                                                                                       aaGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAAttt                                                    <  1:18443/148‑1 (MQ=255)
                                                                                        aGGCCATTATATGTATGCTTCAGAAACACGCACC‑TTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCg                                                 <  1:173811/149‑1 (MQ=255)
                                                                                        aGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTcc                                                  <  3:22830/149‑1 (MQ=255)
                                                                                        aGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTcc                                                  <  7:117777/149‑1 (MQ=255)
                                                                                          gCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGa                                                <  8:105816/149‑1 (MQ=255)
                                                                                          gCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCATGGATCAGTATCTTTGGATGATCCAGAAGTTGCGTAACATTAATAAGGCAATCTATGCAGTCAACAAACAACGGATCGTTCGAATTTcag                                                 >  8:174148/1‑146 (MQ=255)
                                                                                                aGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGAAtata                                           <  6:100670/148‑1 (MQ=255)
                                                                                                 gatggatgCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGAATATAGc                                         >  2:304507/1‑149 (MQ=255)
                                                                                                  atggatgCTTCAGCAACACGCACCTTTAGAGATCAAGACAACTGCTAGCCTGGATCAGTAGCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCACCAAACAACGGATGGGTCGAATTTCCGAATATAGcc                                        <  5:302808/149‑1 (MQ=255)
                                                                                                     gatgCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCTAATTTACGAAGATAgac                                        >  3:66387/1‑144 (MQ=255)
                                                                                                      atgCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATgg                                                             <  1:290024/124‑1 (MQ=255)
                                                                                                      atgCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATgg                                                             >  2:290024/1‑124 (MQ=255)
                                                                                                                    cGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGAATATAGCCGAGACTGCCTGGATGAAc                      <  7:200485/149‑1 (MQ=255)
                                                                                                                       aCCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGATCGAATTTCCGAATATAGCCGAGACTGCCTGGATGAACt                     >  5:329703/1‑147 (MQ=255)
                                                                                                                         cTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGAATATAGCCGAGACTGCCTGGATGAACTCAgg                 >  7:24692/1‑149 (MQ=255)
                                                                                                                         cTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGAATATAGCCGAGACTGCCTGGATGAACTCAgg                 <  6:393177/149‑1 (MQ=255)
                                                                                                                          tttAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGAATATAGACGAGACTGCCTGGAAGAACTCAgtg                >  4:236264/1‑147 (MQ=255)
                                                                                                                            tAGCGATCAAGCCAACGTCTAGCCTGTATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGAATATAGCCGAGACTGCCTGGATGAACTCAGGGc               <  5:83592/148‑1 (MQ=255)
                                                                                                                            tAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGAATATAGCCGAGACTGCCTGGATGAACTCAGGGc               <  8:184199/148‑1 (MQ=255)
                                                                                                                             aGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGAAGATAGCCGAGACTGCCTGGATGAACTCAGGGCTg             >  8:328629/1‑149 (MQ=255)
                                                                                                                              gctaccAAGCCAACGACTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTTCGTAACTGTATTCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGAATATAGCCGAGACTGCCTGGATGAACTCAGGGCTga            <  8:289227/144‑1 (MQ=255)
                                                                                                                              gCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGAATATAGCCGAGACTGCCTGGATGAACTCAGGGCTg             <  8:24692/148‑1 (MQ=255)
                                                                                                                              gCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGAATATAGCCGAGACTGCCTGGATGAACTCAGGGCTg             <  1:304507/148‑1 (MQ=255)
                                                                                                                                 aTCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGAATATAGCCGAGACTGCCTGGATGAACTCAGGGCTGAGAc         >  6:38493/1‑149 (MQ=255)
                                                                                                                                     aGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGAATATAGCCGAGACTGCCTGGATGAACTCAGGGCTGAGACTgg      >  6:4367/1‑148 (MQ=255)
                                                                                                                                     aGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGAATATAGCCGAGACTGCCTGGATGAACTCAGGGCTGAGACTGGt     <  1:188637/149‑1 (MQ=255)
                                                                                                                                        cAACCTATAGCCTGGATCAGTATAGTTGTATGATACAGATGTTGCGTAACTGTAATCCGTCAAGCTATGCAGACAACAAAGAACGGATGGTTCGAATTTCCGAATATAGCCGAGACTGCCTGGATGAACTCAGGGCTGAGACTTGtata  <  7:174148/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
GGCCCCGCGATGGAAACAAGCTTCGCAAACGCAGGGCAGGTTTCTCCGGGGTATGCCTCCCCTTGGTCAGCCCCTGGAGTGCCATTGAAGGCCATTAGATGGATGCTTCAGCAACACGCACCTTTAGCGATCAAGCCAACGTCTAGCCTGGATCAGTATCGTTGGATGATCCAGATGTTGCGTAACTGTAATCCGGCAAGCTATGCAGTCAACAAAGAACGGATGGTTCGAATTTCCGAATATAGCCGAGACTGCCTGGATGAACTCAGGGCTGAGACTGGTATA  >  NC_002947/5031308‑5031592

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