Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,386,276 A→G V638A (GTA→GCA)  aaeB ← p‑hydroxybenzoic acid efflux system component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,386,2760AG98.9% 332.6 / NA 88V638A (GTA→GCA) aaeBp‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (40/47);  minor base T (1/0);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCACGTTTCGCGGCGAGTTCTGCCGTGCATCGCTAAGTATAGCCATGAAAAAACCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTACCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGCAGATATCACATGGTCTGC  >  NC_000913/3386140‑3386406
                                                                                                                                        |                                                                                                                                  
cGTCACGTTTCGCGGCGAGTTCTGCCGTGCATCGCTAAGTATAGCCATGAAAAAACCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGcc                                                                                                                                  <  1:385954/139‑1 (MQ=255)
       tttCGCGGCGAGTTCTGCCGTGCATCGCTAAGTATAGCCATGAAAAAACCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCg                                                                                                                           <  2:1053107/139‑1 (MQ=255)
         tCGCGGCGAGTTCTGCCGTGCATCGCTAAGTATAGCCATGAAAAAACCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTTCCGGTTCCGTc                                                                                                                         >  2:1166247/1‑139 (MQ=255)
                 aGTTCTGCCGTGCATCGCTAAGTATAGCCATGAAAAAACCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCgg                                                                                                                 <  2:1477191/139‑1 (MQ=255)
                      tGCCGTGCATCGCTAAGTATAGCCATGAAAAAACCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGgcg                                                                                                            <  2:1574510/139‑1 (MQ=255)
                      tGCCGTGCATCGCTAAGTATAGCCATGAAAAAACCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGgcg                                                                                                            <  2:1214238/139‑1 (MQ=255)
                         cGTGCATCGCTAAGTATAGCCATGAAAAAACCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCtt                                                                                                         >  1:701892/1‑139 (MQ=255)
                             cATCGCTAAGTATAGCCATGAAAAAACCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTTCGTCACCTGCGGTGGCGCTTGCCa                                                                                                     <  1:1431222/139‑1 (MQ=255)
                                          aGCCATGAAAAAACCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCttt                                                                                        >  2:1491876/1‑139 (MQ=255)
                                            ccATGAAAAAACCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTc                                                                                                                         <  2:47474/104‑1 (MQ=255)
                                            ccATGAAAAAACCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTc                                                                                                                         >  1:47474/1‑104 (MQ=255)
                                             cATGAAAAAACCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTcc                                                                                     <  1:249976/139‑1 (MQ=255)
                                               tGAAAAAACCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTg                                                                                   >  2:512750/1‑139 (MQ=255)
                                                gAAAAAACCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGcc                                                                                                      >  2:672666/1‑119 (MQ=255)
                                                gAAAAAACCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGcc                                                                                                      <  1:672666/119‑1 (MQ=255)
                                                 aaaaaaCCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAATGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGt                                                                                 >  1:112116/1‑139 (MQ=255)
                                                 aaaaaaCCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGt                                                                                 <  2:1220135/139‑1 (MQ=255)
                                                      aCCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGAtt                                                                            >  2:726620/1‑139 (MQ=255)
                                                             cTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACACATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGtt                                                                     <  1:875588/139‑1 (MQ=255)
                                                                              gACTTAACTATCGGTCAATGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCa                                                    <  2:112116/139‑1 (MQ=255)
                                                                                    aCTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGgcg                                                                                                            >  1:1287815/1‑77 (MQ=255)
                                                                                    aCTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGgcg                                                                                                            <  2:1287815/77‑1 (MQ=255)
                                                                                       aTCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACcga                                           <  2:250444/139‑1 (MQ=255)
                                                                                       aTCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACcga                                           <  1:1491876/139‑1 (MQ=255)
                                                                                         cGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCa                                                    <  2:450922/128‑1 (MQ=255)
                                                                                         cGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCa                                                    >  1:450922/1‑128 (MQ=255)
                                                                                          ggTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGcc                                                                                                      >  1:310167/1‑77 (MQ=255)
                                                                                          ggTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGcc                                                                                                      <  2:310167/77‑1 (MQ=255)
                                                                                              aaCGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGt                                                                                                                               >  1:1337265/1‑48 (MQ=255)
                                                                                              aaCGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGt                                                                                                                               <  2:1337265/48‑1 (MQ=255)
                                                                                               aCGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGTTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGttt                                   >  1:252873/1‑139 (MQ=255)
                                                                                                cGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTcc                                                                >  2:1265374/1‑109 (MQ=255)
                                                                                                cGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTa                                  <  2:701892/139‑1 (MQ=255)
                                                                                                cGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCACTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTcc                                                                <  1:1265374/109‑1 (MQ=255)
                                                                                                  cATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACt                                                         <  2:428731/114‑1 (MQ=255)
                                                                                                  cATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACt                                                         >  1:428731/1‑114 (MQ=255)
                                                                                                     gTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGt                                                                                                                          >  1:164639/1‑46 (MQ=255)
                                                                                                     gTTGATACTTATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGt                                                                                                                          <  2:164639/46‑1 (MQ=255)
                                                                                                             ttATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACggg                     <  1:564616/139‑1 (MQ=255)
                                                                                                             ttATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACggg                     >  1:774621/1‑139 (MQ=255)
                                                                                                             ttATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACggg                     <  2:1249035/139‑1 (MQ=255)
                                                                                                              tATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACggg                     >  1:1324323/1‑138 (MQ=255)
                                                                                                              tATGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACggg                     <  2:1324323/138‑1 (MQ=255)
                                                                                                               aTGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACt                                                         >  1:1779443/1‑101 (MQ=255)
                                                                                                               aTGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACt                                                         <  2:1779443/101‑1 (MQ=255)
                                                                                                               aTGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGacgacg                                      <  1:1501532/120‑1 (MQ=255)
                                                                                                               aTGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGacgacg                                      >  2:1501532/1‑120 (MQ=255)
                                                                                                               aTGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGc                    >  1:484366/1‑138 (MQ=255)
                                                                                                               aTGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGc                    <  2:484366/138‑1 (MQ=255)
                                                                                                               aTGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGCa                   >  2:988921/1‑139 (MQ=255)
                                                                                                                tGGAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGCAg                  <  1:1517979/139‑1 (MQ=255)
                                                                                                                 ggAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACggg                     <  2:743291/135‑1 (MQ=255)
                                                                                                                 ggAGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACggg                     >  1:743291/1‑135 (MQ=255)
                                                                                                                 ggACCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGCAGa                 >  2:1213573/1‑139 (MQ=255)
                                                                                                                   aGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCa                                                                                                     <  2:131044/53‑1 (MQ=255)
                                                                                                                   aGCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCa                                                                                                     >  1:131044/1‑53 (MQ=255)
                                                                                                                    gCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAg                                                 <  2:1528722/104‑1 (MQ=255)
                                                                                                                    gCATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAg                                                 >  1:1528722/1‑104 (MQ=255)
                                                                                                                     cATCCCCGCCAGCCGATTTGCCGGTTCCGTCATCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACt                                                         <  2:1060781/95‑1 (MQ=255)
                                                                                                                     cATCCCCGCCAGCCGATTTGCCGGTTCCGTCAGCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACt                                                         >  1:1060781/1‑95 (MQ=255)
                                                                                                                     cATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGc                                                                                                >  2:1779487/1‑56 (MQ=255)
                                                                                                                     cATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGc                                                                                                <  1:1779487/56‑1 (MQ=255)
                                                                                                                     cATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTcc                                                                >  2:1786344/1‑88 (MQ=255)
                                                                                                                     cATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTcc                                                                <  1:1786344/88‑1 (MQ=255)
                                                                                                                     cATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACcgacg                                         >  2:1122874/1‑111 (MQ=255)
                                                                                                                     cATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACcgacg                                         <  1:1122874/111‑1 (MQ=255)
                                                                                                                     cATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACcgac                                          <  1:478747/110‑1 (MQ=255)
                                                                                                                     cATCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACcgac                                          >  2:478747/1‑110 (MQ=255)
                                                                                                                      aTCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACt                                                         >  2:1749081/1‑94 (MQ=255)
                                                                                                                      aTCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACt                                                         <  1:1749081/94‑1 (MQ=255)
                                                                                                                       tCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAg                                                                               >  2:1579610/1‑71 (MQ=255)
                                                                                                                       tCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAg                                                                               <  1:1579610/71‑1 (MQ=255)
                                                                                                                       tCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACggg                     <  2:1639603/129‑1 (MQ=255)
                                                                                                                       tCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACggg                     >  1:1639603/1‑129 (MQ=255)
                                                                                                                       tCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGCAGATATcac           <  1:726620/139‑1 (MQ=255)
                                                                                                                       tCCCCGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGCAGATATcac           >  2:1026342/1‑139 (MQ=255)
                                                                                                                        ccccGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAg                                                                                                    <  1:1138699/49‑1 (MQ=255)
                                                                                                                        ccccGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAg                                                                                                    >  2:1138699/1‑49 (MQ=255)
                                                                                                                        ccccGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGCAg                  >  1:430550/1‑131 (MQ=255)
                                                                                                                        ccccGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGCAg                  <  2:430550/131‑1 (MQ=255)
                                                                                                                        ccccGCCAGCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGCAGATATcaca          >  2:106070/1‑139 (MQ=255)
                                                                                                                                gCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGCAGATATcac           >  1:558404/1‑130 (MQ=255)
                                                                                                                                gCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGCAGATATcac           <  2:558404/130‑1 (MQ=255)
                                                                                                                                gCCGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGCAGATATCACATGGTctgc  <  1:1213573/139‑1 (MQ=255)
                                                                                                                                 ccGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCt                                                                                          >  1:1029216/1‑50 (MQ=255)
                                                                                                                                 ccGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCt                                                                                          <  2:1029216/50‑1 (MQ=255)
                                                                                                                                 ccGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACggg                     >  2:1175532/1‑119 (MQ=255)
                                                                                                                                 ccGATTTGCCGGTTCCGTCACCTGCGGTGGCGCTAGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACggg                     <  1:1175532/119‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                  
CGTCACGTTTCGCGGCGAGTTCTGCCGTGCATCGCTAAGTATAGCCATGAAAAAACCGACGCTTTTGGCGTCGGTTTTGACTTAACTATCGGTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTACCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGCAGATATCACATGGTCTGC  >  NC_000913/3386140‑3386406

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