Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 169,748 T→A 5.7% intergenic (+21/‑30) fhuA → / → fhuC ferrichrome outer membrane transporter/iron(3+)‑hydroxamate import ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913169,7480TA 5.7% 424.4 / 9.2 142intergenic (+21/‑30)fhuA/fhuCferrichrome outer membrane transporter/iron(3+)‑hydroxamate import ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (67/66);  new base A (3/6);  total (70/72)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 145 total aligned reads displayed.

GGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGTGA  >  NC_000913/169599‑169895
                                                                                                                                                     |                                                                                                                                                   
ggctCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAt                                                                                                                                                     >  2:466373/1‑150 (MQ=255)
   tCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTc                                                                                                                                                  <  1:1683787/150‑1 (MQ=255)
    ccAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTcc                                                                                                                                                 >  2:774577/1‑150 (MQ=255)
      aaCGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGATTACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGt                                                                                                                                               >  1:2058723/1‑150 (MQ=255)
        cGTGGCGCTGCATGTTAACAACCTGTTCGATCGGGAATACGTCGCCAGCTGCTTTAACACTAATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGc                                                                                                                                             <  2:280179/150‑1 (MQ=255)
         gTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGGAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGcc                                                                                                                                            <  1:91483/150‑1 (MQ=255)
           ggCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAAATCCGTGcccc                                                                                                                                          <  1:1892135/150‑2 (MQ=255)
           ggCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTAGGGGCACGGATTTCCGTGCCCa                                                                                                                                          >  1:1416800/1‑150 (MQ=255)
            gcgcTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCAt                                                                                                                                         >  1:625663/1‑150 (MQ=255)
            gcgcTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAAATCCGTGcccca                                                                                                                                         <  2:792718/150‑3 (MQ=255)
             cgcTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCAtt                                                                                                                                        >  1:220999/1‑150 (MQ=255)
               cTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTc                                                                                                                                      <  2:168168/150‑1 (MQ=255)
               cTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTc                                                                                                                                      <  1:496907/150‑1 (MQ=255)
               cTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGCACGGAAATCCGTGccccaaaag                                                                                                                                     <  1:1737400/150‑7 (MQ=255)
                tGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTca                                                                                                                                     <  1:773956/150‑1 (MQ=255)
                tGCATGTTAACAACCTGTTCGATCCTGAATCCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGGCAGGTCGTTGCAACCGCAACCATCCGAATCTAATTTCTCTTTTGGGGCACGGAAATCCGTGccccaaaag                                                                                                                                     <  2:146204/150‑7 (MQ=255)
                    tcttAACAACCTGTTCCATCGAGAATACGTCGCCAGCTGCTTTAACACTTATTGCTGCTGCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCATCCTTCCGTTTCTAATTTCTCTTTTGTGGCACGGATTTCCGTGCCCATTTCACAAg                                                                                                                                 <  2:1070671/148‑1 (MQ=255)
                    tGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAg                                                                                                                                 >  2:561301/1‑150 (MQ=255)
                     gTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGt                                                                                                                                >  1:1025018/1‑150 (MQ=255)
                      ttAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGtt                                                                                                                               >  1:1848656/1‑150 (MQ=255)
                       tAACAACCTGTTCGATCGCGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTTCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTg                                                                                                                              <  2:416783/150‑1 (MQ=255)
                        aacaacCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTgg                                                                                                                             >  2:1098965/1‑150 (MQ=255)
                         acaacCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGc                                                                                                                            >  2:1374739/1‑150 (MQ=255)
                         acaacCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGc                                                                                                                            <  1:1229668/150‑1 (MQ=255)
                          caacCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCt                                                                                                                           <  2:1249304/150‑1 (MQ=255)
                             ccTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTGTCTTTTGGGGCACGGAATTCCGTGCCCATTACACAAGTTGCCTGtt                                                                                                                        >  1:393215/1‑150 (MQ=255)
                              cTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTa                                                                                                                       >  1:542919/1‑150 (MQ=255)
                                gTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATg                                                                                                                     >  1:1625696/1‑150 (MQ=255)
                                  tCGATCGTGAATACGTCGCCAGCTGCTTTCACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCa                                                                                                                   <  2:2204202/150‑1 (MQ=255)
                                  tCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTATTTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCa                                                                                                                   >  1:1945118/1‑150 (MQ=255)
                                    gATCGTGACTACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTTCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAgg                                                                                                                 <  1:903809/150‑1 (MQ=255)
                                     aTCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGa                                                                                                                >  1:1548952/1‑150 (MQ=255)
                                      tCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGaa                                                                                                               <  1:2029736/150‑1 (MQ=255)
                                        gttaaTACGTCGCAAGCTGCTTTAACACTTTTGGCTGCTTCTGGGGCGAAGAACGTCAGGTCGTTTCAAGCGCAACCTTCCCTTTCTAATTTCTCTTTTGGGGCACGGTTTTCCGTGCCCATTTCAGAAGTTGGCTGTTATGCATTAATa                                                                                                             <  2:2073107/147‑1 (MQ=255)
                                         ggAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATac                                                                                                            <  2:77289/149‑1 (MQ=255)
                                           aaTACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCGGTTTGGCAGGAAGacac                                                                                                          >  1:880337/1‑150 (MQ=255)
                                             tACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGa                                                                                                        <  2:862835/150‑1 (MQ=255)
                                               cGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCACCCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAAt                                                                                                      >  2:1773028/1‑150 (MQ=255)
                                               cGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAAt                                                                                                      >  1:957284/1‑150 (MQ=255)
                                                gTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATc                                                                                                     <  2:820484/150‑1 (MQ=255)
                                                gTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATc                                                                                                     >  2:51602/1‑150 (MQ=255)
                                                   gCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCAtt                                                                                                  >  1:2057386/1‑150 (MQ=255)
                                                   gCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCAtt                                                                                                  >  1:40848/1‑150 (MQ=255)
                                                    ccAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTc                                                                                                 <  2:625663/150‑1 (MQ=255)
                                                      aGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCg                                                                                               <  2:1785653/150‑1 (MQ=255)
                                                          gctTTAACACTTATGGCTGCTTCTGTGGCGCAGAACGTCAGGTCGTTGCAACAGCAACATTAAGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCACATTTCACAAGTTGTCTGTTATGCAGGAATACACGAATCATTCCGATAc                                                                                           <  1:1435613/150‑1 (MQ=255)
                                                          gctTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGTGGCACGGATTTCCGTGCCCATTTCACAAGTAGGCTGTTATGCAGGAATACACGAATCATTCCGATAc                                                                                           >  1:1059844/1‑150 (MQ=255)
                                                           ctTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCACAGTTCACCAGTTGGCTGTTATGCAGGAATACACGAATAAATcccctaac                                                                                          >  2:1116360/1‑144 (MQ=255)
                                                             ttAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCAc                                                                                        >  2:1620452/1‑150 (MQ=255)
                                                               aaCACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACtt                                                                                      >  2:938475/1‑150 (MQ=255)
                                                               aaCACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACtt                                                                                      <  1:25368/150‑1 (MQ=255)
                                                               aaCACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCAATCCGAAACCACtt                                                                                      >  1:692799/1‑150 (MQ=255)
                                                                   cTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGTGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGc                                                                                  <  1:1582835/150‑1 (MQ=255)
                                                                   cTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTTGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGc                                                                                  <  2:1567640/150‑1 (MQ=255)
                                                                          ctgctTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGt                                                                           <  2:1815342/150‑1 (MQ=255)
                                                                             ctTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAat                                                                        >  2:373958/1‑150 (MQ=255)
                                                                             ctTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCAGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGATTCATTCCGATACCACTTTTGCACTGCGTAat                                                                        <  1:249230/150‑1 (MQ=255)
                                                                              ttCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAata                                                                       <  1:2027583/150‑1 (MQ=255)
                                                                               tCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAatat                                                                      >  1:1841746/1‑150 (MQ=255)
                                                                                cTGGGGCGCAGAACGTCAGCTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCAATTCACAAGTTGGCTGTTAAGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGAAataac                                                                     >  2:934106/1‑148 (MQ=255)
                                                                                  ggggCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATATCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATAtctc                                                                   >  2:339777/1‑150 (MQ=255)
                                                                                   gggCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTcc                                                                  <  1:1332244/150‑1 (MQ=255)
                                                                                      cgcAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCttt                                                               >  1:1802922/1‑150 (MQ=255)
                                                                                      cgcAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCttt                                                               >  1:140549/1‑150 (MQ=255)
                                                                                       gcAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTTCACTGCGTAATATATCCTTTc                                                              >  2:573233/1‑150 (MQ=255)
                                                                                       gcAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTTATATCTCCTTTc                                                              >  2:658523/1‑150 (MQ=255)
                                                                                       gcAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTc                                                              >  2:1451600/1‑150 (MQ=255)
                                                                                       gcAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATATCCTTTc                                                              >  2:1599536/1‑150 (MQ=255)
                                                                                        cAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCg                                                             <  1:711288/150‑1 (MQ=255)
                                                                                        cAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCg                                                             <  2:1841746/150‑1 (MQ=255)
                                                                                        cAGAACGTCAGGTCGTTGAAACCGCAACCTTCAGTTTCTAATTTCTCTTTTGGGGCTCGGATTTCCGTGCCCATTTCACAAGTTTGCTGATATGCAGGAATACACGAATCAATCCGATACCACTTTTGCACTGCGTTATATCTCCTTTCg                                                             >  2:507587/1‑150 (MQ=255)
                                                                                          taaCGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCgtg                                                           <  2:740152/149‑1 (MQ=255)
                                                                                            aCGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGtgtg                                                         >  1:2192864/1‑150 (MQ=255)
                                                                                            aCGTCAGGTCGTTGCAACCGCAACCGTCGGTTTCTAATTTCTCTTTTGGGGCCCGGATTTCCGTGACCATTTCACAAGTTCGCTGTTATGCAGGAATACACGAATCATTCCGATACCTCTTTTGCACTGCGTAATATCTCCTTTCGtgtg                                                         <  1:344742/150‑1 (MQ=255)
                                                                                                           aaccgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGCTACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtt                                          >  1:1206780/1‑150 (MQ=255)
                                                                                                            accgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTTACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCttt                                         <  2:2205957/150‑1 (MQ=255)
                                                                                                              cgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTg                                       <  2:1014304/150‑1 (MQ=255)
                                                                                                              cgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATAAACGAATCATTCCGATACAACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGTGCGCACGATTTTg                                       >  2:1178440/1‑150 (MQ=255)
                                                                                                               gcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGc                                      <  2:1180921/150‑1 (MQ=255)
                                                                                                               gcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGCGCCCGGGCACACgaaaggag                                      <  1:1165609/150‑8 (MQ=255)
                                                                                                                 aaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCAt                                    >  1:1107877/1‑150 (MQ=255)
                                                                                                                  accTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATc                                   >  2:1863260/1‑150 (MQ=255)
                                                                                                                   ccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTTCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATcc                                  >  1:1142415/1‑150 (MQ=255)
                                                                                                                   ccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATcc                                  <  1:2105213/150‑1 (MQ=255)
                                                                                                                          tgtCTATTTTCGCTTTTCGGGCACGGATATCCGTGCCCGTTTCTCAAGTTGGCTGTTATTCAGGACATCACGAATCAATCCGATACCTCTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCg                           <  2:359258/148‑1 (MQ=255)
                                                                                                                            tCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGtt                         <  2:1316825/150‑1 (MQ=255)
                                                                                                                               aaTTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAAc                      >  1:341387/1‑150 (MQ=255)
                                                                                                                                aTTTCTCTTTTGGTGCACGGATTTCCGTGAACATTTCACAAGTTTGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAAcc                     <  1:2097789/150‑1 (MQ=255)
                                                                                                                                aTTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGTATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAAcc                     <  2:2026434/150‑1 (MQ=255)
                                                                                                                                aTTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACAGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAAcc                     >  2:2075520/1‑150 (MQ=255)
                                                                                                                                   tctcTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCttt                  >  1:1541254/1‑150 (MQ=255)
                                                                                                                                    ctctTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGTAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTc                 >  2:37314/1‑150 (MQ=255)
                                                                                                                                      tgaaatgggCACGGAAATCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCt               >  1:22024/7‑150 (MQ=255)
                                                                                                                                      tgaaatgggCACGGAAATCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCt               >  1:780438/7‑150 (MQ=255)
                                                                                                                                           ggggCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTGATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCgg          <  1:1294960/150‑1 (MQ=255)
                                                                                                                                            gggCACGGATTTCCGTGCCCATTTCACTAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCTCTGTCGTTAACCTTTCCTGCCggg         <  2:863800/150‑1 (MQ=255)
                                                                                                                                            gggCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCggg         <  2:1229616/150‑1 (MQ=255)
                                                                                                                                             ggCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGa        <  2:1945118/150‑1 (MQ=255)
                                                                                                                                              gCACGGATTTCCGTGCCCCTTTCACAAGGTGGCTGTTATGCAGGAATACACGCATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCCCGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGaa       >  2:104490/1‑150 (MQ=255)
                                                                                                                                                  ggATTTCCGTGCCCATTTCACAAGTTGGCTCTTATGCAGGAATACACGAATCAATCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCCGGCAGGAAAggt   <  2:1425849/150‑3 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGTGA  >  NC_000913/169599‑169895

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: