Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 172,499 T→C 6.9% H346H (CAT→CAC fhuB → iron(3+)‑hydroxamate import ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913172,4990TC 6.9% 386.4 / 13.5 115H346H (CAT→CACfhuBiron(3+)‑hydroxamate import ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (52/55);  new base C (4/4);  total (56/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01

Only 100 of 117 total aligned reads displayed.

CTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATG  >  NC_000913/172353‑172646
                                                                                                                                                  |                                                                                                                                                   
cTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACAtgt                                                                                                                                                   <  1:460619/149‑1 (MQ=255)
  gACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGc                                                                                                                                                 <  4:1380580/149‑1 (MQ=255)
        cgTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGt                                                                                                                                           <  2:976277/149‑1 (MQ=255)
          tgtgGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGttt                                                                                                                                         <  4:116327/149‑1 (MQ=255)
            tggatggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGc                                                                                                                                       <  4:100708/149‑1 (MQ=255)
             ggatggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGcc                                                                                                                                      <  3:19207/149‑1 (MQ=255)
             ggatggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGcc                                                                                                                                      <  2:244000/149‑1 (MQ=255)
              gatggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGc                                                                                                                                                 <  2:63/137‑1 (MQ=255)
              gatggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGc                                                                                                                                                 >  1:63/1‑137 (MQ=255)
              gatggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGttt                                                                                                                                         <  2:1261309/145‑1 (MQ=255)
              gatggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGttt                                                                                                                                         >  1:1261309/1‑145 (MQ=255)
              gatggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGctc                                                                                                                                     <  2:777654/149‑3 (MQ=255)
                 ggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTcg                                                                                                                                  >  3:536521/1‑149 (MQ=255)
                    aGTGTCCACCGGTTCGGTCACTGCGATGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCgcgg                                                                                                                               <  4:1009757/149‑1 (MQ=255)
                       gtCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCgcgggcg                                                                                                                            >  1:4789/1‑149 (MQ=255)
                          cACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAGCACATGTTGGCg                                                                                                                                            >  3:713655/1‑130 (MQ=255)
                          cACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAGCACATGTTGGCGttc                                                                                                                                         <  1:1063918/133‑2 (MQ=255)
                          cACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAGCACATGTTGGCGttc                                                                                                                                         >  2:1063918/1‑132 (MQ=255)
                           aCCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGt                                                                                                                        >  3:132217/1‑149 (MQ=255)
                             cGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGtgc                                                                                                                      >  3:571937/1‑149 (MQ=255)
                              ggTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAGCACATGTTGGCGt                                                                                                                                           <  3:138210/127‑1 (MQ=255)
                              ggTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAGCACATGTTGGCGt                                                                                                                                           >  4:138210/1‑127 (MQ=255)
                               gTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGtgctg                                                                                                                    >  3:664479/1‑149 (MQ=255)
                                    gTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTg                                                                                                                                               <  3:334401/117‑1 (MQ=255)
                                    gTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTg                                                                                                                                               >  4:334401/1‑117 (MQ=255)
                                    gTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGtt                                                                                                               <  4:329192/149‑1 (MQ=255)
                                      cACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGa                                                                                                             <  4:664479/149‑1 (MQ=255)
                                      cACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGa                                                                                                             <  4:536521/149‑1 (MQ=255)
                                         tGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATgg                                                                                                          <  2:651594/149‑1 (MQ=255)
                                          gCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATgg                                                                                                          >  1:468289/1‑148 (MQ=255)
                                          gCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATgg                                                                                                          <  2:468289/148‑1 (MQ=255)
                                              tGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCtgtg                                                                                                     <  2:4789/149‑1 (MQ=255)
                                                aTCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTgtggt                                                                                                   >  3:565403/1‑149 (MQ=255)
                                                 tCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCtgtg                                                                                                     >  2:244802/1‑146 (MQ=255)
                                                 tCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCtgtg                                                                                                     <  1:244802/146‑1 (MQ=255)
                                                    gTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTgtggtggtg                                                                                               >  2:583410/1‑149 (MQ=255)
                                                      gcgcCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGAt                                                                                                            >  4:1134729/1‑134 (MQ=255)
                                                      gcgcCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGAt                                                                                                            <  3:1134729/134‑1 (MQ=255)
                                                            ctgctACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGctgctg                                                                                                                 <  4:796718/123‑1 (MQ=255)
                                                            ctgctACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGctgctg                                                                                                                 >  3:796718/1‑123 (MQ=255)
                                                                tACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGttt                                                                                   <  4:565403/149‑1 (MQ=255)
                                                                  cTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTgg                                                                                 <  4:132217/149‑1 (MQ=255)
                                                                    gtgGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTc                                                                               <  3:232734/149‑1 (MQ=255)
                                                                     tgGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCg                                                                              >  1:806131/1‑149 (MQ=255)
                                                                              ccGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTc                                                                               <  4:1093832/139‑1 (MQ=255)
                                                                              ccGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTc                                                                               >  3:1093832/1‑139 (MQ=255)
                                                                              ccGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTg                                                                            <  2:11562/142‑1 (MQ=255)
                                                                              ccGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTg                                                                            >  1:11562/1‑142 (MQ=255)
                                                                                gcgTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACg                                                                   >  1:735253/1‑149 (MQ=255)
                                                                                gcgTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACg                                                                   >  1:1170867/1‑149 (MQ=255)
                                                                                  gTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTg                                                                                                              >  1:660001/1‑104 (MQ=255)
                                                                                  gTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGc                                                                 <  1:583410/149‑1 (MQ=255)
                                                                                  gTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGc                                                                 >  2:1249485/1‑149 (MQ=255)
                                                                                        gcagcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCtgtg                                                                                                     >  4:417724/1‑107 (MQ=255)
                                                                                        gcagcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTgtggtg                                                                                                  <  4:802352/110‑1 (MQ=255)
                                                                                        gcagcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTgtggtg                                                                                                  >  3:802352/1‑110 (MQ=255)
                                                                                        gcagcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGttt                                                                                   >  4:188460/1‑125 (MQ=255)
                                                                                        gcagcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGttt                                                                                   <  3:188460/125‑1 (MQ=255)
                                                                                        gcagcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGa                                                                           >  3:1215215/1‑133 (MQ=255)
                                                                                        gcagcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGa                                                                           <  4:1215215/133‑1 (MQ=255)
                                                                                        gcagcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACg                                                           >  1:1100516/1‑149 (MQ=255)
                                                                                        gcagcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACg                                                           >  1:746254/1‑149 (MQ=255)
                                                                                          agcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGa                                                                           <  3:850724/131‑1 (MQ=255)
                                                                                          agcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGAGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGa                                                                           >  4:850724/1‑131 (MQ=255)
                                                                                           gcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAGCACATGTTGGCGt                                                                                                                                           <  4:625947/66‑1 (MQ=38)
                                                                                           gcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTc                                                                                       <  1:232702/118‑1 (MQ=255)
                                                                                           gcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTc                                                                                       >  2:232702/1‑118 (MQ=255)
                                                                                           gcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTgg                                                        <  3:1177284/149‑1 (MQ=255)
                                                                                              ttAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTgtggtgg                                                                                                 <  3:93587/105‑1 (MQ=255)
                                                                                              ttAGCGCGCCGGATATGAAGGCCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTgtggtgg                                                                                                 >  4:93587/1‑105 (MQ=255)
                                                                                                    cgcCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCggg                                               <  2:746254/149‑1 (MQ=255)
                                                                                                      ccGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGc                                             <  4:571937/149‑1 (MQ=255)
                                                                                                            atGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGt                                                                                                                <  4:396019/76‑1 (MQ=255)
                                                                                                            atGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGt                                                                                                                >  3:396019/1‑76 (MQ=255)
                                                                                                            atGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGttt                                                                                   <  3:48617/105‑1 (MQ=255)
                                                                                                            atGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGttt                                                                                   >  4:48617/1‑105 (MQ=255)
                                                                                                             tGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGttt                                                                                   >  2:1184805/1‑104 (MQ=255)
                                                                                                             tGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGttt                                                                                   <  1:1184805/104‑1 (MQ=255)
                                                                                                              gAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCg                                     >  1:1297372/1‑149 (MQ=255)
                                                                                                               aaGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTg                                                                                  <  1:92414/103‑1 (MQ=255)
                                                                                                               aaGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTg                                                                                  >  2:92414/1‑103 (MQ=255)
                                                                                                               aaGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGa                                    <  2:1170867/149‑1 (MQ=255)
                                                                                                                  gTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTc                                                                               >  4:1390235/1‑103 (MQ=255)
                                                                                                                   tCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTgg                                                        >  3:535615/1‑125 (MQ=255)
                                                                                                                   tCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAg                                   >  1:1245184/1‑146 (MQ=255)
                                                                                                                       cGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTaa                            <  2:1100516/149‑1 (MQ=255)
                                                                                                                        gATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAg                                   <  2:63903/141‑1 (MQ=255)
                                                                                                                        gATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAg                                   >  1:63903/1‑141 (MQ=255)
                                                                                                                        gATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAAt                           >  1:211330/1‑149 (MQ=255)
                                                                                                                         aTCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTg                                                               >  4:1120435/1‑112 (MQ=255)
                                                                                                                                    gCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTg               <  2:806131/149‑1 (MQ=255)
                                                                                                                                       gAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGcc            >  4:592701/1‑149 (MQ=255)
                                                                                                                                         aCGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGcc                        <  2:1414056/135‑1 (MQ=255)
                                                                                                                                         aCGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGcc                        >  1:1414056/1‑135 (MQ=255)
                                                                                                                                             cAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGtggtggtgg                                                                                              >  3:446091/1‑61 (MQ=255)
                                                                                                                                             cAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGtggtggtgg                                                                                              <  4:446091/61‑1 (MQ=255)
                                                                                                                                             cAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAAt      >  1:510910/1‑149 (MQ=255)
                                                                                                                                              aaCATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAAtt     >  3:333238/1‑149 (MQ=255)
                                                                                                                                               aCATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTa    >  4:732754/1‑149 (MQ=255)
                                                                                                                                                 aTGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATg  >  1:331539/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATG  >  NC_000913/172353‑172646

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