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→GATsodCsuperoxide 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

GATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACC  >  NC_000913/1724252‑1724541
                                                                                                                                                |                                                                                                                                                  
gATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCcgc                                                                                                                                               >  4:46112/1‑149 (MQ=255)
 aTTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCcgcg                                                                                                                                              <  2:978585/149‑1 (MQ=255)
   ttCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCGCGGAATCcgcaggc                                                                                                                                           >  2:73008/1‑145 (MQ=255)
     cAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCGCGGAATCcgcaggcgg                                                                                                                                         <  1:73008/149‑7 (MQ=255)
       gACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTg                                                                                                                                        >  2:713864/1‑149 (MQ=255)
       gACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTg                                                                                                                                        >  2:82428/1‑149 (MQ=255)
       gACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTg                                                                                                                                        >  4:318468/1‑149 (MQ=255)
       gACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCAGCGCCGCGGAATCcgcaggcgg                                                                                                                                         >  2:6172/1‑143 (MQ=255)
       gACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCAGCGCCGCGGAATCcgcaggcgg                                                                                                                                         <  1:741657/149‑7 (MQ=255)
         cGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGc                                                                                                                                      <  2:904562/149‑1 (MQ=255)
            ggCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGTCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCGCGGAATCcgcaggcgg                                                                                                                                         <  1:984430/142‑7 (MQ=255)
            ggCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCGCGGAATCcgcaggcgg                                                                                                                                         >  2:984430/1‑136 (MQ=255)
              cgcgATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGc                                                                                                                                 <  3:318468/149‑1 (MQ=255)
               gcgATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGcc                                                                                                                                <  4:368612/149‑1 (MQ=255)
                  aTGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATc                                                                                                                             <  2:78962/149‑1 (MQ=255)
                    gACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCtt                                                                                                                           >  4:994718/1‑149 (MQ=255)
                       ggCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtgg                                                                                                                        <  2:305538/149‑1 (MQ=255)
                        gCATCGGTAGCTTTGTCGTCATTATTGACGAACAGTGCAGGCAGATCGTCTAAATGCCCGGCACCTTATGGCCCTTCATGTTTACCGGTAGTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGACATCTTtggt                                                                                                                       <  4:172500/149‑1 (MQ=255)
                         cATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtggtg                                                                                                                      >  3:59575/1‑149 (MQ=255)
                         cATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTtggtg                                                                                                                      >  3:157487/1‑149 (MQ=255)
                            cGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggct                                                                                                                   <  4:424526/149‑1 (MQ=255)
                            cGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggct                                                                                                                   <  4:59575/149‑1 (MQ=255)
                            cGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggct                                                                                                                   <  4:157487/149‑1 (MQ=255)
                             ggTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctg                                                                                                                  >  2:537707/1‑149 (MQ=255)
                                aGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctggct                                                                                                               <  4:140671/149‑1 (MQ=255)
                                 gCTTTGCCGTCATGATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctgggtg                                                                                                              <  1:46580/149‑1 (MQ=255)
                                   tttGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctggctggc                                                                                                            >  2:246096/1‑149 (MQ=255)
                                   tttGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGtggctgactggc                                                                                                            >  1:358323/1‑149 (MQ=255)
                                      gCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGc                                                                                                         <  2:921480/149‑1 (MQ=255)
                                        cGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTGCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGATGGCTGGCAGCtt                                                                                                       >  1:353453/1‑149 (MQ=255)
                                            attattGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCtt                                                                                                   <  4:832530/149‑1 (MQ=255)
                                             ttattGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCttt                                                                                                  >  1:378568/1‑149 (MQ=255)
                                                  gacgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCAt                                                                                             <  3:46112/149‑1 (MQ=255)
                                                   acgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCt                                                                                                    >  1:709784/1‑141 (MQ=255)
                                                   acgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCt                                                                                                    <  2:709784/141‑1 (MQ=255)
                                                   acgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCAt                                                                                             <  3:856710/148‑1 (MQ=255)
                                                   acgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCAt                                                                                             >  4:856710/1‑148 (MQ=255)
                                                   acgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATg                                                                                            >  3:100566/1‑149 (MQ=255)
                                                    cgacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGa                                                                                           <  4:100566/149‑1 (MQ=255)
                                                     gacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGc                                                                                                                                 >  2:694225/1‑110 (MQ=255)
                                                     gacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGc                                                                                                                                 <  1:694225/110‑1 (MQ=255)
                                                     gacCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGaa                                                                                          <  1:246096/149‑1 (MQ=255)
                                                        cAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCAGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAatat                                                                                       <  1:752517/149‑1 (MQ=255)
                                                                gcagATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATg                                                                                      <  2:340884/142‑1 (MQ=255)
                                                                gcagATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATg                                                                                      >  1:340884/1‑142 (MQ=255)
                                                                gcagATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCa                                                                               >  1:915366/1‑149 (MQ=255)
                                                                    aTCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCttt                                                                                                  <  4:917325/126‑1 (MQ=255)
                                                                    aTCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCttt                                                                                                  >  3:917325/1‑126 (MQ=255)
                                                                    aTCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGtt                                                                           >  3:140958/1‑149 (MQ=255)
                                                                         cTAGATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGACGCCATGTTCACCg                                                                      >  4:726541/1‑149 (MQ=255)
                                                                           aaaTGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCttt                                                                                                  <  1:815329/119‑1 (MQ=255)
                                                                           aaaTGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCttt                                                                                                  >  2:815329/1‑119 (MQ=255)
                                                                            aaTGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCgggg                                                                   <  3:150624/149‑1 (MQ=255)
                                                                               gCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGt                                                                >  2:1002860/1‑149 (MQ=255)
                                                                                 ccGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTaa                                                              >  1:789617/1‑149 (MQ=255)
                                                                                   ggCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATg                                                            >  1:696151/1‑149 (MQ=255)
                                                                                   ggCACCGTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATg                                                            >  4:283889/1‑149 (MQ=255)
                                                                                       ccTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCttt                                                        >  3:40435/1‑149 (MQ=255)
                                                                                        cTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTc                                                       <  1:82428/149‑1 (MQ=255)
                                                                                        cTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTc                                                       <  4:168660/149‑1 (MQ=255)
                                                                                         ttCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCa                                                      <  4:40435/149‑1 (MQ=255)
                                                                                           cTGGCGCTTCATGTTTAGCGGTATTTTGTGGATCAAGTTGCGCGCTTGCGGATTCCGCGGCGGTGGTTTTGACATCTTTGGTGGCTGACTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGa                                                    <  2:358323/149‑1 (MQ=255)
                                                                                             ggCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATc                                                  <  2:300656/149‑1 (MQ=255)
                                                                                             ggCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATc                                                  <  1:1002860/149‑1 (MQ=255)
                                                                                               cccTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCgg                                                >  2:85405/1‑149 (MQ=255)
                                                                                               cccTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCgg                                                <  4:624830/149‑1 (MQ=255)
                                                                                                ccTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCggg                                               <  2:353453/149‑1 (MQ=255)
                                                                                                      tGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGaaaa                                         <  1:713864/149‑1 (MQ=255)
                                                                                                       gTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAatat                                                                                       <  2:411641/102‑1 (MQ=255)
                                                                                                       gTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAatat                                                                                       >  1:411641/1‑102 (MQ=255)
                                                                                                       gTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAAc                                        >  2:629970/1‑149 (MQ=255)
                                                                                                                    ttGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCgggggg                                                                 <  3:3510/111‑1 (MQ=255)
                                                                                                                    ttGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCgggggg                                                                 >  4:3510/1‑111 (MQ=255)
                                                                                                                        ggATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTa                            <  2:638220/144‑1 (MQ=255)
                                                                                                                        ggATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTa                            >  1:638220/1‑144 (MQ=255)
                                                                                                                        ggATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGt                       >  4:997858/1‑149 (MQ=255)
                                                                                                                        ggATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGt                       <  1:537707/149‑1 (MQ=255)
                                                                                                                          aTCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGttt                     >  4:481789/1‑149 (MQ=255)
                                                                                                                           tCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTc                    >  4:508863/1‑149 (MQ=255)
                                                                                                                           tCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTc                    <  2:789617/149‑1 (MQ=255)
                                                                                                                           tCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTc                    <  2:786458/149‑1 (MQ=255)
                                                                                                                            cAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCa                   <  2:378568/149‑1 (MQ=255)
                                                                                                                                    ccGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGt           >  1:307634/1‑149 (MQ=255)
                                                                                                                                       ccTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGAc        <  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)
                                                                                                                                                |                                                                                                                                                  
GATTTCAGACGAGGCGCGATGACGGCATCGGTAGCTTTGCCGTCATTATTGACGACCAGTGCAGGCAGATCGCCTAAATGCCCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACC  >  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: