Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,724,395 T→A 6.8% E87D (GAA→GAT)  sodC ← superoxide dismutase, Cu, Zn, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,724,3950TA 6.8% 300.0 / 11.3 88E87D (GAA→GAT) sodCsuperoxide dismutase, Cu, Zn, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (38/44);  new base A (3/3);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

GATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACC  >  NC_000913/1724252‑1724541
                                                                                                                                                |                                                                                                                                                  
gATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCcgc                                                                                                                                               >  4:46112/1‑149 (MQ=255)
 aTTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCcgcg                                                                                                                                              <  2:978585/149‑1 (MQ=255)
   ttCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑C‑GCGGAATCcgcaggc                                                                                                                                           >  2:73008/1‑145 (MQ=255)
     cAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑C‑GCGGAATCcgcaggcgg                                                                                                                                         <  1:73008/149‑7 (MQ=255)
       gACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTg                                                                                                                                        >  2:713864/1‑149 (MQ=255)
       gACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTg                                                                                                                                        >  2:82428/1‑149 (MQ=255)
       gACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTg                                                                                                                                        >  4:318468/1‑149 (MQ=255)
       gACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCAGCGCCGCGGAATCcgcaggcgg                                                                                                                                         >  2:6172/1‑143 (MQ=255)
       gACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCAGCGCCGCGGAATCcgcaggcgg                                                                                                                                         <  1:741657/149‑7 (MQ=255)
         cGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGc                                                                                                                                      <  2:904562/149‑1 (MQ=255)
            ggCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGTCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑C‑GCGGAATCcgcaggcgg                                                                                                                                         <  1:984430/142‑7 (MQ=255)
            ggCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑C‑GCGGAATCcgcaggcgg                                                                                                                                         >  2:984430/1‑136 (MQ=255)
              cgcgATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGc                                                                                                                                 <  3:318468/149‑1 (MQ=255)
               gcgATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGcc                                                                                                                                <  4:368612/149‑1 (MQ=255)
                  aTGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATc                                                                                                                             <  2:78962/149‑1 (MQ=255)
                    gACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCtt                                                                                                                           >  4:994718/1‑149 (MQ=255)
                       ggCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtgg                                                                                                                        <  2:305538/149‑1 (MQ=255)
                        gCATCGGTAGCTTTGTCGTCATTATTGACGAACAGTGCAGGCAGATCGTCTAAATGCCCGGCACCTTATGGCCCTTCATGTTTACCGGTAGTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGACATCTTtggt                                                                                                                       <  4:172500/149‑1 (MQ=255)
                         cATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtggtg                                                                                                                      >  3:59575/1‑149 (MQ=255)
                         cATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtggtg                                                                                                                      >  3:157487/1‑149 (MQ=255)
                            cGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggct                                                                                                                   <  4:424526/149‑1 (MQ=255)
                            cGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggct                                                                                                                   <  4:59575/149‑1 (MQ=255)
                            cGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggct                                                                                                                   <  4:157487/149‑1 (MQ=255)
                             ggTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctg                                                                                                                  >  2:537707/1‑149 (MQ=255)
                                aGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctggct                                                                                                               <  4:140671/149‑1 (MQ=255)
                                 gCTTTGCCGTCATGATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctgggtg                                                                                                              <  1:46580/149‑1 (MQ=255)
                                   tttGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctggctggc                                                                                                            >  2:246096/1‑149 (MQ=255)
                                   tttGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctgactggc                                                                                                            >  1:358323/1‑149 (MQ=255)
                                      gCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGc                                                                                                         <  2:921480/149‑1 (MQ=255)
                                        cGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTGCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGATGGCTGGCAGCtt                                                                                                       >  1:353453/1‑149 (MQ=255)
                                            attattGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCtt                                                                                                   <  4:832530/149‑1 (MQ=255)
                                             ttattGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCttt                                                                                                  >  1:378568/1‑149 (MQ=255)
                                                  gacgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCAt                                                                                             <  3:46112/149‑1 (MQ=255)
                                                   acgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCt                                                                                                    >  1:709784/1‑141 (MQ=255)
                                                   acgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCt                                                                                                    <  2:709784/141‑1 (MQ=255)
                                                   acgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCAt                                                                                             <  3:856710/148‑1 (MQ=255)
                                                   acgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCAt                                                                                             >  4:856710/1‑148 (MQ=255)
                                                   acgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATg                                                                                            >  3:100566/1‑149 (MQ=255)
                                                    cgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGa                                                                                           <  4:100566/149‑1 (MQ=255)
                                                     gacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGc                                                                                                                                 >  2:694225/1‑110 (MQ=255)
                                                     gacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGc                                                                                                                                 <  1:694225/110‑1 (MQ=255)
                                                     gacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGaa                                                                                          <  1:246096/149‑1 (MQ=255)
                                                        cAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCAGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAatat                                                                                       <  1:752517/149‑1 (MQ=255)
                                                                gcagATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATg                                                                                      <  2:340884/142‑1 (MQ=255)
                                                                gcagATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATg                                                                                      >  1:340884/1‑142 (MQ=255)
                                                                gcagATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCa                                                                               >  1:915366/1‑149 (MQ=255)
                                                                    aTCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCttt                                                                                                  <  4:917325/126‑1 (MQ=255)
                                                                    aTCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCttt                                                                                                  >  3:917325/1‑126 (MQ=255)
                                                                    aTCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGtt                                                                           >  3:140958/1‑149 (MQ=255)
                                                                         cTAGATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGACGCCATGTTCACCg                                                                      >  4:726541/1‑149 (MQ=255)
                                                                           aaaTGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCttt                                                                                                  <  1:815329/119‑1 (MQ=255)
                                                                           aaaTGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCttt                                                                                                  >  2:815329/1‑119 (MQ=255)
                                                                            aaTGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCgggg                                                                   <  3:150624/149‑1 (MQ=255)
                                                                               gCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGt                                                                >  2:1002860/1‑149 (MQ=255)
                                                                                 ccGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTaa                                                              >  1:789617/1‑149 (MQ=255)
                                                                                   ggCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATg                                                            >  1:696151/1‑149 (MQ=255)
                                                                                   ggCACCGTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATg                                                            >  4:283889/1‑149 (MQ=255)
                                                                                       ccTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCttt                                                        >  3:40435/1‑149 (MQ=255)
                                                                                        cTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTc                                                       <  1:82428/149‑1 (MQ=255)
                                                                                        cTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTc                                                       <  4:168660/149‑1 (MQ=255)
                                                                                         ttCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCa                                                      <  4:40435/149‑1 (MQ=255)
                                                                                           cTGGCGCTTCATGTTTAGCGGTATTTTGTGGATCAAGTTGCGCGC‑TTGCGGATTCCGCGGCGGTGGTTTTGACATCTTTGGTGGCTGACTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGa                                                    <  2:358323/149‑1 (MQ=255)
                                                                                             ggCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATc                                                  <  2:300656/149‑1 (MQ=255)
                                                                                             ggCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATc                                                  <  1:1002860/149‑1 (MQ=255)
                                                                                               cccTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCgg                                                >  2:85405/1‑149 (MQ=255)
                                                                                               cccTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCgg                                                <  4:624830/149‑1 (MQ=255)
                                                                                                ccTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCggg                                               <  2:353453/149‑1 (MQ=255)
                                                                                                      tGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGaaaa                                         <  1:713864/149‑1 (MQ=255)
                                                                                                       gTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAatat                                                                                       <  2:411641/102‑1 (MQ=255)
                                                                                                       gTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAatat                                                                                       >  1:411641/1‑102 (MQ=255)
                                                                                                       gTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAAc                                        >  2:629970/1‑149 (MQ=255)
                                                                                                                    ttGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCgggggg                                                                 <  3:3510/111‑1 (MQ=255)
                                                                                                                    ttGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCgggggg                                                                 >  4:3510/1‑111 (MQ=255)
                                                                                                                        ggATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTa                            <  2:638220/144‑1 (MQ=255)
                                                                                                                        ggATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTa                            >  1:638220/1‑144 (MQ=255)
                                                                                                                        ggATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGt                       >  4:997858/1‑149 (MQ=255)
                                                                                                                        ggATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGt                       <  1:537707/149‑1 (MQ=255)
                                                                                                                          aTCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGttt                     >  4:481789/1‑149 (MQ=255)
                                                                                                                           tCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTc                    >  4:508863/1‑149 (MQ=255)
                                                                                                                           tCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTc                    <  2:789617/149‑1 (MQ=255)
                                                                                                                           tCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTc                    <  2:786458/149‑1 (MQ=255)
                                                                                                                            cAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCa                   <  2:378568/149‑1 (MQ=255)
                                                                                                                                    ccGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGt           >  1:307634/1‑149 (MQ=255)
                                                                                                                                       c‑cTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGAc        <  1:629970/149‑1 (MQ=255)
                                                                                                                                          tGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGc      <  2:915366/149‑1 (MQ=255)
                                                                                                                                              gATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTAcc  >  1:999391/1‑149 (MQ=255)
                                                                                                                                               aTTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCgggg                                                                   >  2:834243/1‑83 (MQ=255)
                                                                                                                                               aTTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCgggg                                                                   <  1:834243/83‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGC‑CTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACC  >  NC_000913/1724252‑1724541

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑