Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,607,367 C→G 5.5% A306G (GCC→GGC)  zntA → zinc, cobalt and lead efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,607,3670CG 5.5% 263.0 / 5.0 74A306G (GCC→GGC) zntAzinc, cobalt and lead efflux system
Reads supporting (aligned to +/- strand):  ref base C (36/33);  new base G (2/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

GGCGATGTGATTGAAGTCGCCGCAGGTGGGCGTTTGCCTGCCGACGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCGTCGCGCTCCCATTGAGCGGTTTATCGACCGTTTCAGCCGTATCTA  >  NC_000913/3607222‑3607511
                                                                                                                                                 |                                                                                                                                                
ggCGATGTGATTGAAGTCGCCGCAGGTGGGCGTTTGCCTGCCGACGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGccac                                                                                                                                               <  4:275317/149‑1 (MQ=255)
       tgATTGAAGTCGCCGCAGGTGGGCGTTTGCCTGCCGACGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTa                                                                                                                                        >  3:412148/1‑149 (MQ=255)
          ttGAAGTCGCCGCAGGTGGGCGTTTGCCTGCCGACGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGAc                                                                                                                                     >  4:134813/1‑149 (MQ=255)
             aaGTCGCCGCAGGTGGGCGTTTGCCTGCCGACGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGt                                                                                                                                  >  1:479189/1‑149 (MQ=255)
                  gccgcAGGTGGGCGTTTGCCTGCCGACGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTgtt                                                                                                                             >  1:367474/1‑147 (MQ=255)
                  gccgcAGGTGGGCGTTTGCCTGCCGACGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGt                                                                                                                             >  4:397970/1‑149 (MQ=255)
                   ccgcAGGTGGGCGTTTGCCTGCCGACGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTg                                                                                                                            >  4:296442/1‑149 (MQ=255)
                   ccgcAGGTGGGCGTTTGCCTGCCGACGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTg                                                                                                                            <  1:355036/149‑1 (MQ=255)
                      cAGGTGGGCGTTTGCCTGCCGACGGTAAACTGCTCTCACCGTTTGCCAGCTTTGATGAAAGCGCCCTGACCGCCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACg                                                                                                                         >  1:613134/1‑149 (MQ=255)
                       aGGTGGGCGTTTGCCTGCCGACGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGt                                                                                                                        <  1:571852/149‑1 (MQ=255)
                          tGGGCGTTTGCCTGCCGACGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTgg                                                                                                                     >  3:65909/1‑149 (MQ=255)
                          tGGGCGTTTGCCTGCCGACGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTgg                                                                                                                     >  2:263399/1‑149 (MQ=255)
                            ggCGTTTGCCTGCCGACGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGaa                                                                                                                   <  3:351448/149‑1 (MQ=255)
                               gTTTGCCTGCCGACGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTg                                                                                                                <  2:415588/149‑1 (MQ=255)
                                tttGCCTGCCGACGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGc                                                                                                               >  4:294222/1‑149 (MQ=255)
                                     ctgccGACGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCa                                                                                                          <  3:397970/149‑1 (MQ=255)
                                           aCGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCg                                                                                                    >  2:240878/1‑149 (MQ=255)
                                              gTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGa                                                                                                 <  4:412148/149‑1 (MQ=255)
                                                      ctctCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAgcgc                                                                                         >  3:366093/1‑149 (MQ=255)
                                                       tctcACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTc                                                                                                           >  1:611206/1‑130 (MQ=255)
                                                       tctcACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTc                                                                                                           <  2:611206/130‑1 (MQ=255)
                                                       tctcACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGcc                                                                                        >  3:546687/1‑149 (MQ=255)
                                                        ctcACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGGCACCAGcagg                                                                                                                                        <  2:294882/100‑4 (MQ=255)
                                                        ctcACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGGCACCAGcagg                                                                                                                                        <  4:539283/100‑4 (MQ=255)
                                                        ctcACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGGCACCAGcagg                                                                                                                                        >  1:294882/1‑97 (MQ=255)
                                                        ctcACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGGCACCAGcagg                                                                                                                                        >  3:539283/1‑97 (MQ=255)
                                                         tcACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCAt                                                                                      >  4:241978/1‑149 (MQ=255)
                                                          cACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCAtt                                                                                     >  1:75344/1‑149 (MQ=255)
                                                          cACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCAtt                                                                                     >  3:326767/1‑149 (MQ=255)
                                                           aCCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTg                                                                                                                               <  2:440157/106‑1 (MQ=255)
                                                           aCCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTg                                                                                                                               >  1:440157/1‑106 (MQ=255)
                                                             cGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGAc                                                                                  >  1:77068/1‑149 (MQ=255)
                                                                 tGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCt                                                                                                              >  1:161562/1‑117 (MQ=255)
                                                                 tGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCt                                                                                                              <  2:161562/117‑1 (MQ=255)
                                                                 tGCCAGTTGTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGa                                                                              <  3:136003/149‑1 (MQ=255)
                                                                   ccAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGAtt                                                                            >  1:262320/1‑149 (MQ=255)
                                                                    cAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTc                                                                           >  1:614509/1‑149 (MQ=255)
                                                                       ttttGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGa                                                                        >  3:527873/1‑149 (MQ=255)
                                                                       ttttGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGa                                                                        >  1:503617/1‑149 (MQ=255)
                                                                         ttGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGaaa                                                                      <  4:527995/149‑1 (MQ=255)
                                                                            atgaAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGa                                                                              >  3:475570/1‑138 (MQ=255)
                                                                            atgaAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGa                                                                              <  4:475570/138‑1 (MQ=255)
                                                                            atgaAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTg                                                                   <  4:366093/149‑1 (MQ=255)
                                                                             tgaAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGa                                                                  <  3:294222/149‑1 (MQ=255)
                                                                             tgaAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGa                                                                  <  3:134813/149‑1 (MQ=255)
                                                                              gaAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGAt                                                                 <  2:652335/149‑1 (MQ=255)
                                                                               aaaGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGaa                                                                       <  1:458124/142‑1 (MQ=255)
                                                                               aaaGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGaa                                                                       >  2:458124/1‑142 (MQ=255)
                                                                                   cgcCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCgaag                                                            <  4:326767/149‑1 (MQ=255)
                                                                                      ccTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGaagaag                                                         <  4:65909/149‑1 (MQ=255)
                                                                                           ccGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGaa                                                    <  2:503617/149‑1 (MQ=255)
                                                                                            cGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAag                                                   >  3:268000/1‑149 (MQ=255)
                                                                                               cGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGc                                                <  3:241978/149‑1 (MQ=255)
                                                                                                gAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCg                                               >  1:162574/1‑149 (MQ=255)
                                                                                                gAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCg                                               >  2:37004/1‑149 (MQ=255)
                                                                                                    ccATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGaaa                                                                      >  1:663585/1‑122 (MQ=255)
                                                                                                    ccATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGaaa                                                                      <  2:663585/122‑1 (MQ=255)
                                                                                                            gTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCGTCGCGCTCCCAt                                   >  4:635448/1‑149 (MQ=255)
                                                                                                                 gcgcgcAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCGTCGCGCTCCCATTGAGc                              <  2:479189/149‑1 (MQ=255)
                                                                                                                 gcgcgcAACGGGCGATAAAGTCCCTGCTGGTGACACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGATGAAGCCGAAGAGCGTCGCGCTCCCATTGAGc                              <  3:591756/149‑1 (MQ=255)
                                                                                                                     gcAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCGTCGCGCTCCCATTGAGCGGtt                          <  3:296442/149‑1 (MQ=255)
                                                                                                                       aaCGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCGTCGCGCTCCCATTGAGCGGTTTa                        <  3:156090/149‑1 (MQ=255)
                                                                                                                          gggCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCGTCGCGCTCCCATTGAGCGGTTTATCg                     >  3:326206/1‑149 (MQ=255)
                                                                                                                           ggCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCGTCGCGCTCCCATTGAGCGGTTTATCGa                    <  3:635448/149‑1 (MQ=255)
                                                                                                                            gCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCGTCGCGCTCCCATTGAGCGGTTTATCGAc                   <  2:367474/149‑1 (MQ=255)
                                                                                                                               aTAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGAt                                                                             <  4:684847/88‑1 (MQ=255)
                                                                                                                               aTAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGAt                                                                             >  3:684847/1‑88 (MQ=255)
                                                                                                                                tAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCGTCGCGCTCCCATTGAGCGGTTTATCGACCGtt               <  3:175512/149‑1 (MQ=255)
                                                                                                                                    gTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCGTCGCGCTCCCATTGAGCGGTTTATCGACCGTTTCAg           <  1:240878/149‑1 (MQ=255)
                                                                                                                                       ccTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCGTCGCGCTCCCATTGAGCGGTTTATCGACCGTTTCAGCCg        >  2:620425/1‑149 (MQ=255)
                                                                                                                                        ctgctgGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCGTCGCGCTCCCATTGAGCGGTTTATCGACCGTTTCAGCCGt       >  2:205433/1‑149 (MQ=255)
                                                                                                                                         tgctgGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCGTCGCGCTCCCATTGAGCGGTTTATCGACCGTTTCAGCCGTa      <  2:75344/149‑1 (MQ=255)
                                                                                                                                          gctgGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCGTCGCGCTCCCATTGAGCGGTTTATCGACCGTTTCAGCCGTAt     <  1:37004/149‑1 (MQ=255)
                                                                                                                                           ctgGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCGTCGCGCTCCCATTGAGCGGTTTATCGACCGTTTCAGCCGTATc    >  4:482419/1‑149 (MQ=255)
                                                                                                                                             ggTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCGTCGCGCTCCCATTGAGCGGTTTATCGACCGTTTCAGCCGTATCta  <  2:614509/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGCGATGTGATTGAAGTCGCCGCAGGTGGGCGTTTGCCTGCCGACGGTAAACTGCTCTCACCGTTTGCCAGTTTTGATGAAAGCGCCCTGACCGGCGAATCCATTCCGGTGGAGCGCGCAACGGGCGATAAAGTCCCTGCTGGTGCCACCAGCGTAGACCGTCTGGTGACGTTGGAAGTGCTGTCAGAACCGGGAGCCAGCGCCATTGACCGGATTCTGAAACTGATCGAAGAAGCCGAAGAGCGTCGCGCTCCCATTGAGCGGTTTATCGACCGTTTCAGCCGTATCTA  >  NC_000913/3607222‑3607511

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: