Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 172,503 C→T 6.9% L348L (CTG→TTG)  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,5030CT 6.9% 378.2 / 12.6 115L348L (CTG→TTG) fhuBiron(3+)‑hydroxamate import ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (52/55);  new base T (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 = 9.99e-01

Only 100 of 118 total aligned reads displayed.

GACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGC  >  NC_000913/172355‑172648
                                                                                                                                                    |                                                                                                                                                 
gACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGc                                                                                                                                                   <  4:1380580/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                                                                                                                                           >  1:1261309/1‑145 (MQ=255)
            gatggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGttt                                                                                                                                           <  2:1261309/145‑1 (MQ=255)
            gatggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGctc                                                                                                                                       <  2:777654/149‑3 (MQ=255)
               ggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTcg                                                                                                                                    >  3:963270/1‑149 (MQ=255)
               ggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTcg                                                                                                                                    >  4:232734/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)
                           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)
                                   tCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTg                                                                                                                >  4:1177284/1‑149 (MQ=255)
                                    cACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGa                                                                                                               <  4:536521/149‑1 (MQ=255)
                                    cACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGa                                                                                                               <  4:664479/149‑1 (MQ=255)
                                      cTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATg                                                                                                             >  2:329720/1‑149 (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                                                                                                       <  1:244802/146‑1 (MQ=255)
                                                  gTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTgtggtggtg                                                                                                 >  2:583410/1‑149 (MQ=255)
                                                    gcgcCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGAt                                                                                                              <  3:1134729/134‑1 (MQ=255)
                                                    gcgcCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGAt                                                                                                              >  4:1134729/1‑134 (MQ=255)
                                                          ctgctACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGctgctg                                                                                                                   <  4:796718/123‑1 (MQ=255)
                                                              tACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGttt                                                                                     <  4:565403/149‑1 (MQ=255)
                                                                  gtgGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTc                                                                                 <  3:232734/149‑1 (MQ=255)
                                                                   tgGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCg                                                                                >  1:806131/1‑149 (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:1170867/1‑149 (MQ=255)
                                                                              gcgTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACg                                                                     >  1:735253/1‑149 (MQ=255)
                                                                                gTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTg                                                                                                                <  2:660001/104‑1 (MQ=255)
                                                                                gTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTg                                                                                                                >  1:660001/1‑104 (MQ=255)
                                                                                gTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGc                                                                   >  2:1249485/1‑149 (MQ=255)
                                                                                gTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGc                                                                   <  1:583410/149‑1 (MQ=255)
                                                                                      gcagcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCtgtg                                                                                                       <  3:417724/107‑1 (MQ=255)
                                                                                      gcagcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCtgtg                                                                                                       >  4:417724/1‑107 (MQ=255)
                                                                                      gcagcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTgtggtg                                                                                                    >  3:802352/1‑110 (MQ=255)
                                                                                      gcagcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTgtggtg                                                                                                    <  4:802352/110‑1 (MQ=255)
                                                                                      gcagcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGttt                                                                                     <  3:188460/125‑1 (MQ=255)
                                                                                      gcagcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGttt                                                                                     >  4:188460/1‑125 (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)
                                                                                         gcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAGCACATGTTGGCGt                                                                                                                                             >  3:625947/1‑66 (MQ=37)
                                                                                         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                                                                                                                  >  3:396019/1‑76 (MQ=255)
                                                                                                          atGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGt                                                                                                                  <  4:396019/76‑1 (MQ=255)
                                                                                                          atGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGttt                                                                                     >  4:48617/1‑105 (MQ=255)
                                                                                                          atGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGttt                                                                                     <  3:48617/105‑1 (MQ=255)
                                                                                                           tGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGttt                                                                                     <  1:1184805/104‑1 (MQ=255)
                                                                                                           tGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGttt                                                                                     >  2:1184805/1‑104 (MQ=255)
                                                                                                            gAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCg                                       >  1:1297372/1‑149 (MQ=255)
                                                                                                             aaGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTg                                                                                    >  2:92414/1‑103 (MQ=255)
                                                                                                             aaGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGa                                      <  2:1170867/149‑1 (MQ=255)
                                                                                                                gTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTc                                                                                 >  4:1390235/1‑103 (MQ=255)
                                                                                                                gTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTc                                                                                 <  3:1390235/103‑1 (MQ=255)
                                                                                                                 tCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTgg                                                          >  3:535615/1‑125 (MQ=255)
                                                                                                                 tCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAg                                     >  1:1245184/1‑146 (MQ=255)
                                                                                                                     cGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTaa                              <  2:1100516/149‑1 (MQ=255)
                                                                                                                      gATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAg                                     >  1:63903/1‑141 (MQ=255)
                                                                                                                      gATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAAt                             >  1:211330/1‑149 (MQ=255)
                                                                                                                       aTCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTg                                                                 <  3:1120435/112‑1 (MQ=255)
                                                                                                                       aTCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTg                                                                 >  4:1120435/1‑112 (MQ=255)
                                                                                                                                  gCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTg                 <  2:806131/149‑1 (MQ=255)
                                                                                                                                     gAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGcc              >  4:592701/1‑149 (MQ=255)
                                                                                                                                       aCGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGcc                          >  1:1414056/1‑135 (MQ=255)
                                                                                                                                       aCGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGcc                          <  2:1414056/135‑1 (MQ=255)
                                                                                                                                           cAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGtggtggtgg                                                                                                <  4:446091/61‑1 (MQ=255)
                                                                                                                                           cAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGtggtggtgg                                                                                                >  3:446091/1‑61 (MQ=255)
                                                                                                                                           cAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAAt        >  1:510910/1‑149 (MQ=255)
                                                                                                                                             aCATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTa      >  4:732754/1‑149 (MQ=255)
                                                                                                                                               aTGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATg    >  1:331539/1‑149 (MQ=255)
                                                                                                                                                 gtgCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATggc  <  2:1297372/149‑1 (MQ=255)
                                                                                                                                                  tgCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGa                                      >  4:279770/1‑112 (MQ=255)
                                                                                                                                                  tgCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGa                                      <  3:279770/112‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGC  >  NC_000913/172355‑172648

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: