Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,258,702 G→A 6.3% D321N (GAC→AAC)  dinF → oxidative stress resistance protein; putative MATE family efflux pump; UV and mitomycin C inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,258,7020GA 6.3% 207.1 / 4.5 63D321N (GAC→AAC) dinFoxidative stress resistance protein; putative MATE family efflux pump; UV and mitomycin C inducible protein
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base A (2/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

ACCGTACTTGGCGCGCGACTGGGGAGTGACATTATCGCTGTTAACGCGGTTCTGATGACGCTACTCACCTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAACACATCATTGCGTTACTGACGTCGTTAACCCAGATTCAGCA  >  NC_000913/4258561‑4258838
                                                                                                                                             |                                                                                                                                        
aCCGTACTTGGCGCGCGACTGGGGAGTGACATTATCGCTGTTAACGCGGTTCTGATGACGCTACTCACCTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTa                                                                                                                                    <  4:507271/148‑1 (MQ=255)
                     gggAGTGACATTATCGCTGTTAACGCGGTTCTGATGACGCTACTCACCTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGt                                                                                                                              <  4:856595/133‑1 (MQ=255)
                     gggAGTGACATTATCGCTGTTAACGCGGTTCTGATGACGCTACTCACCTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGt                                                                                                                              >  3:856595/1‑133 (MQ=255)
                               ttATCGCTGTTAACGCGGTTCTGATGACGCTACTCACCTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTg                                                                                                     >  1:35482/1‑148 (MQ=255)
                                 aTCGCTGTTAACGCGGTTCTGATGACGCTACTCACCTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTgcc                                                                                                   <  2:35482/148‑1 (MQ=255)
                                   cGCTGTTAACGCGGTTCTGATGACGCTACTCACCTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCa                                                                                                          <  4:739651/139‑1 (MQ=255)
                                   cGCTGTTAACGCGGTTCTGATGACGCTACTCACCTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCa                                                                                                          >  3:739651/1‑139 (MQ=255)
                                                      atgaCGCTACTCACCTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCAACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGt                                                                              >  3:260903/1‑148 (MQ=255)
                                                      atgaCGCTACTCACCTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCAACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGt                                                                              >  3:441181/1‑148 (MQ=255)
                                                       tgaCGCTACTCACCTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCg                                                                                                       <  1:1109868/122‑1 (MQ=255)
                                                       tgaCGCTACTCACCTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCg                                                                                                       >  2:1109868/1‑122 (MQ=255)
                                                         aCGCTACTCACCTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTAc                                                                           >  2:409028/1‑148 (MQ=255)
                                                                 cACCTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCg                                                                   >  2:532613/1‑148 (MQ=255)
                                                                 cACCTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCg                                                                   >  2:536893/1‑148 (MQ=255)
                                                                   ccTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCg                                                                                           <  3:726188/122‑1 (MQ=255)
                                                                   ccTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCg                                                                                           >  4:726188/1‑122 (MQ=255)
                                                                   ccTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCggt                                                                 >  1:452831/1‑148 (MQ=255)
                                                                   ccTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCAACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCggt                                                                 <  4:260903/148‑1 (MQ=255)
                                                                   ccTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCAACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCggt                                                                 <  4:441181/148‑1 (MQ=255)
                                                                     tttACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTg                                                                                                     >  3:433799/1‑110 (MQ=255)
                                                                     tttACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTg                                                                                                     <  4:433799/110‑1 (MQ=255)
                                                                            ccTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCa                                                                                                                                <  1:410066/76‑1 (MQ=255)
                                                                            ccTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCa                                                                                                                                >  2:410066/1‑76 (MQ=255)
                                                                             cTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTg                                                       <  2:452831/148‑1 (MQ=255)
                                                                              tATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGcgcgcg                                                                                                                                          >  1:782767/1‑64 (MQ=255)
                                                                              tATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGcgcgcg                                                                                                                                          <  2:782767/64‑1 (MQ=255)
                                                                                      ggATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCgg                                                                  <  2:992238/128‑1 (MQ=255)
                                                                                      ggATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCgg                                                                  >  1:992238/1‑128 (MQ=255)
                                                                                        aTGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGa                                            >  1:336575/1‑148 (MQ=255)
                                                                                           gTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAaca                                         <  2:336575/148‑1 (MQ=255)
                                                                                            ttttGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCgggg                                                                                        >  2:546766/1‑100 (MQ=255)
                                                                                            ttttGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCgggg                                                                                        <  1:546766/100‑1 (MQ=255)
                                                                                            ttttGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAacac                                        >  1:687531/1‑148 (MQ=255)
                                                                                             tttGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTgctggctg                                                <  4:928465/139‑1 (MQ=255)
                                                                                             tttGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAAcaca                                       >  1:806237/1‑148 (MQ=255)
                                                                                             tttGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGGTTTCGGTGGTTTATTTgctggctg                                                >  3:928465/1‑139 (MQ=255)
                                                                                                 ccTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAACAcatcat                                   <  1:409028/148‑1 (MQ=255)
                                                                                                 ccTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGATGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAACAcatcat                                   <  2:806237/148‑1 (MQ=255)
                                                                                                  cTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCgg                                                                  <  1:288616/116‑1 (MQ=255)
                                                                                                  cTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCgg                                                                  >  2:288616/1‑116 (MQ=255)
                                                                                                   tACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGAt                                                                                      >  1:1052569/1‑95 (MQ=255)
                                                                                                   tACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGAt                                                                                      <  2:1052569/95‑1 (MQ=255)
                                                                                                    aCGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGtttt                                                                     <  4:890113/111‑1 (MQ=255)
                                                                                                    aCGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGtttt                                                                     >  3:890113/1‑111 (MQ=255)
                                                                                                    aCGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGctggctgg                                               <  4:1002152/133‑1 (MQ=255)
                                                                                                    aCGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGctggctgg                                               >  3:1002152/1‑133 (MQ=255)
                                                                                                         gTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCggt                                                                 >  3:444044/1‑110 (MQ=255)
                                                                                                         gTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCggt                                                                 <  4:444044/110‑1 (MQ=255)
                                                                                                          ttGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTTGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCggt                                                                 <  2:918274/109‑1 (MQ=255)
                                                                                                          ttGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCggt                                                                 >  1:918274/1‑109 (MQ=255)
                                                                                                              aGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCg                                                                   <  3:837129/103‑1 (MQ=255)
                                                                                                              aGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCg                                                                   >  4:837129/1‑103 (MQ=255)
                                                                                                                    cTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAACACATCATTGCGTTACTGACGTCGTTa                <  1:536893/148‑1 (MQ=255)
                                                                                                                      ccGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAaca                                         <  4:979421/121‑1 (MQ=255)
                                                                                                                      ccGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAaca                                         >  3:979421/1‑121 (MQ=255)
                                                                                                                        ggTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAACACATCATTGCGTTACTGACGTCGTTAAccc            <  1:532613/148‑1 (MQ=255)
                                                                                                                           cAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCggtggt                                                              <  4:719936/95‑1 (MQ=255)
                                                                                                                           cAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAACACATCATTGCGTTACTGACGTCGTTAAcc             <  4:592201/144‑1 (MQ=255)
                                                                                                                           cAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAACACATCATTGCGTTACTGACGTCGTTAAcc             >  3:592201/1‑144 (MQ=255)
                                                                                                                           cAGGCATACGGTGCGCGCGACGGTAGCCAGTTACTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCggtggt                                                              >  3:719936/1‑95 (MQ=255)
                                                                                                                               cttACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCg                                                                                           <  4:1086630/60‑1 (MQ=255)
                                                                                                                               cATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCg                                                                                           >  3:1086630/1‑62 (MQ=255)
                                                                                                                               cATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAACACATCATTGCGTTACTGAcgtcgt                  >  2:111840/1‑135 (MQ=255)
                                                                                                                               cATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAACACATCATTGCGTTACTGAcgtcgt                  <  1:111840/135‑1 (MQ=255)
                                                                                                                                 tACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAACACATCATTGCGTTACTGACGTCGTTAACCCAGATTcagc   <  2:687531/148‑1 (MQ=255)
                                                                                                                                            cgACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAACACATCATTGCGTTACTGACGTCGTTAACCCAGATTcagca  <  4:3235/138‑1 (MQ=255)
                                                                                                                                            cgACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAACACATCATTGCGTTACTGACGTCGTTAACCCAGATTcagca  >  3:3235/1‑138 (MQ=255)
                                                                                                                                             |                                                                                                                                        
ACCGTACTTGGCGCGCGACTGGGGAGTGACATTATCGCTGTTAACGCGGTTCTGATGACGCTACTCACCTTTACCGCCTATGCGCTGGATGGTTTTGCCTACGCGGTTGAAGCGCACTCCGGTCAGGCATACGGTGCGCGCGACGGTAGCCAGTTGCTGGATGTCTGGCGGGCAGCGTGCCGCCAGTCGGGGATCGTAGCGTTACTGTTTTCGGTGGTTTATTTGCTGGCTGGGGAACACATCATTGCGTTACTGACGTCGTTAACCCAGATTCAGCA  >  NC_000913/4258561‑4258838

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